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FOREWORD 


The  work  described  in  this  report  was  accomplished  by  members  of  the 
Dunham  Laboratory,  Yale  University  under  subcontract  to  the  SUBIC 
program  (contract  NOnr  2512(00))  for  which  the  Office  of  Naval  Re¬ 
search  is  the  sponsor  and  General  Dynamics/Electric  Boat  the  prime 
contractor.  Cdr  F.  R.  Haselton,  Jr.,  USN,  is  Project  Officer  for 
ONR;  Dr.  A.  J.  van  Woerkom  is  Proje  t  Coordinator  for  Electric  Boat 
and  Chief  Scientist  of  the  Applied  Sciences  Department, 

The  SUBIC  program  encompasses  all  aspects  of  submarine  system  analysis. 
This  report  is  one  of  a  series  dealing  with  acoustic  signal  processing. 
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ABSTRACT 


This  report  describes  work  concerned  with  the  detection  of  a  single 
target  in  an  isotropic  noise  field.  In  each  case  studied,  the  data 
source  was  assumed  to  be  a  given  array  of  omnidirectional  hydrophones, 
but  many  different  schemes  for  processing  the  various  hydrophone  out¬ 
puts  were  considered.  The  problems  investigated  involved  comparison 
of  optimum  detectors  with  standard  detectors  and  comparison  of  detec¬ 
tors  operating  on  clipped  and  unclipped  hydrophone  data,  or,  more 
generally,  hydrophone  data  that  have  been  operated  upon  by  an  arbi¬ 
trary  non-linear  device.  The  results  of  these  comparisons  are  de¬ 
scribed  in  detail. 
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Hit?  following  is  a  summary  of  work  porfonnod  under  contract  53-00-10-0231 
between  Yale  University  and  the  Electric  Boat  Company  during  the  period 
21  January  1963  to  1  July  1963.  Detailed  results  and  supporting  arguments 
for  the  various  conclusions  are  contained  in  a  series  of  seven  progress 
reports  that  are  appended . 

With  minor  exceptions  the  work  reported  on  was  all  concerned  with  the 
detection  of  a  single  target  in  an  isotropic  noise  field,,  In  each  case  the 
data  source  was  assumed  to  be  a  given  array  of  omnidirectional  hydrophones, 
but  many  different  schemes  for  processing  the  various  hydrophone  outputs 
were  considered.  The  basic  assumptions  concerning  the  properties  of  noise 
and  target  signal  fell  into  one  of  two  general  categories; 

A.  Both  signal  and  noise  ere  independent  stationary  Gaussian  random 
variables;  the  signal  is  directional,  the  noise  isotropic.  In 
this  case  the  statistical  properties  of  signal  and  noise  are 
completely  specified  by  a  statement  of  total  power  and  spectral 
properties. 

B,  Signal  and  noise  have  one  component  satisfying  condition  A, 

In  addition  each  may  contain  one  or  more  periodic  components  in 
the  low  frequency  band  0  <  f  <  cps.  Furthermore  the  random 
component  of  the  signal  in  the  high  frequency  band  <  f  sf  Wg 
cps  may  be  modulated  in  amplitude  by  one  of  the  periodic  low 
frequency  components  of  the  signal. 

Assumption  A  clearly  constitutes  a  minimal  description  of  signal  and  noise 
and  therefore  leads  to  detection  schemes  of  only  moderate  efficiency.  On 
the  other  hand,  knowledge  of  power  and  spectral  properties  may  be  very  nearly 
all  that  is  available;  in  fact,  even  this  knowledge  may  be  quite  imperfect. 
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.The  best  detection  scheme  under  assumption  A  therefore  sets  realistic  bounds 
on  attainable  detectability  in  many  situations  of  practical  interest. 

Assumption  B  represents  an  attempt  to  obtain  more  refined  detection 
procedures  by  specifying  additional  signal  properties  that  do  occur  in 
important  practical  cases.  The  3,ignal  emitted  by  a  moving  ship  contains 
periodic  or  quasl-periodic  low  frequency  components  generated  by  the 
propeller  and  other  mechanical  components.  Amplitude  modulation  of  the 
high  frequency  range  at  the  propeller  frequency  is  also  a  commonly 
observed  phenomenon  in  such  situations.  Periodic  or  quasi-periodic  noise 
components  are  often  generated  by  moving  machinery  on  the  observing  ship0 
The  problems  investigated  may  also  be  divided  into  two  broad  classes: 

a)  Comparison  of  optimum  detectors  with  standard  detectors.  The 
term  "optimum  detector"  is  interpreted  as  a  likelihood  ratio 
detector  operating  on  the  outputs  of  an  array  of  hydrophones 
mechanically  or  electrically  steered  "on  target,"  The  "standard 
detector"  operates  on  the  output  of  the  same  hydrophone  array  by 
adding,  squaring  and  then  smoothing  by  means  of  a  low  pass 
filter.  A  block  diagram  of  a  standard  detector  is  shown  in 
Fig,  1  of  Progress  Report  No,  3c- 

b)  Comparison  of  detectors  operating  on  clipped  and  unclipped 
hydrophone  data,  or,  more  generally,  hydrophone  data  that  have 
been  operated  upon  by  an  arbitrary  nonlinear  device.  Where  no 
remarks  to  the  contrary  appear,  comparisons  are  always  based 

on  rms  or  mean  square  signal  to  noise  ratio  at  the  detector  output. 
The  input  signal  to  noise  ratio  at  each  hydrophone  is  assumed  to 


Results 


l)  Comparison  of  optimum  and  standard  detector  under  n,.",ur,|,i  . '  \ 
This  forms  the  subject  matter  of  Progress  Report  Ho,  3.  Using  r<r: 
given  by  Bryn^  the  figure  of  merit  (rms  signal  to  noise  ratio)  of  the 


Optimum  detector  for  reasonably  long  observation  times  can  be  written  in 


T  is  the  observation  time  in  seconds,  co^  and  are  the  lower  and  upper 
endpoints  of  the  processed  frequency  interval  in  rad/, sec*  S («)  and  N(w) 
are  the  input  signal  and  noise  spectra  at  each  hydrophone,  G(co)  is  the 
array  gain,  a  quantity  introduced  by  Bryn  and  defined  as  the  ratio  of  the 
contributions  of  signal  and  noise  to  the  average  detector  output  normalized 
with  respect  to  the  input  signal  tc  noise  ratio  at  each  hydrophone. 

Several  features  of  the  above  expression  are  significant.  The  factor 
Yt'  indicates  the  usual  dependence  of  statistical  fluctuation  on  sample 
size.  The  fact  that  the  integrand  is  non-negative  shews  that  the  detector 
should  utilize  all  frequencies  where  S(oo)  does  not  vanish  identically. 

In  cases  of  great  practical  interest  G(co)  varies  only  slowly  and  over  a 
moderate  total  range.  Hence  the  greatest  contribution  to  the  integral 
is  made  by  frequency  ranges  where  the  input  signal  to  noise  ratio  is 
relatively  large.  It  is  clear  that  changes  in  signal  or  noise  spectra  are 
important  only  to  the  extent  that  they  affect  the  integral.  Thus  knowledge 
of  the  detailed  structure  of  signal  and  noise  spectra  is  of  limited  value 


^F,  Bryn,  "Optimal  Signal  Processing  by  Three-Dimensional  Arrays 
Operating  on  Gaussian  Signals  and  Noise,"  J„  Acoust.0  Sac.  Am.,  vol.  3U, 
no.  3,  March  1962  pp.  289-297, 
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in  detection,  Even  the  knowledge  that  the  signal  c  ontains  a  very  narrow- 
band  component,  so  that  S(co)/N(co)  is  large  over  a  narrow  band,  contributes 
little  to  the  solution  of  the  detection  problem  unless  the  narrow-band 
component  contains  a  significant  fraction  of  the  total  signal  power  (or 
the  noise  power  level  happens  to  be  exceptionally  low  in  the  frequency 
range  of  the  narrow-band  signal). 

The  array  gain  G(co)  is  larger  for  the  optimum  detector  than  for  the 

standard  detector  using  the  same  set  of  hydrophones,  particularly  at  low 

frequencies.  This  effect  may  be  attributed  qualitatively  to  the  circumstance 

that  the  optimum  detector  tends  to  reduce  the  effect  of  noise  correlation 

between  different  hydrophones.  The  phenomenon  is  clearly  most  pronounced 

at  low  frequencies Thus  the  optimum  detector  is  superior  to  the  standard 

detector  in  two  respects:  l)  It  can  combat  noise  correlation  between 

hydrophones  and  2)  it  can  utilize  variations  in  input  signal  to  noise  ratio 

2 

over  the  processed  frequency  bands  Significant  differences  in  performance 
between  the  two  types  of  detectors  will  arise  only  when  one  or  both  of 
these  factors  are  important. 

If  signal  and  noise  spectra  are  both  flat  over  a  frequency  band 
0  <  co  <  co  and  vanish  for  co  >  to  ,  and  if  co  is  sufficiently  large  so 
that  the  noise  correlation  between  hydrophones  is  negligible,  neither  one 
of  the  above  factors  is  operative  and  the  performance  of  optimum  and 
standard  detectors  is  identical,  ^Report  No,  Eq,  (55)).  If  the  noise 

^At  high  frequencies,  where  there  is  no  noise  correlation  between 
hydrophones,  the  array  gain  simply  equals  the  number  of  elements  in  the 
array, 

2 

The  standard  detector  is  sensitive  primarily  to  total  power-.  Spectral 
properties  affect  its  design  only  to  the  extent  that  they  determine  what 
frequency  range  should  be  processed. 
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has  the  specified  properties  but  the  signal  deviates  from  the  flat 
spectrum,  the  optimum  detector  is  better  try  an  amount  depending  on  the 
extent  of  the  deviation  |  Eq.  (50 )  „ 

If  signal  and  noise  spectra  are  not  bandlirrdted,  the  figure  of  merit 
of  the  optimal  detector  can  have  a  peculiar  property.  Consider  for 
example  the  case  (quite  closely  approximated  in  situations  of  practical 
interest)  of  signal  and  noise  spectra  identical  in  shape.  Since  0(<o) 
approaches  a  constant  equal  to  the  number  of  hydrophones  as  co-^<»  , 
the  figure  of  merit  clearly  approaches  infinity  as  .  This  says 

that  detection  can  be  accomplished  to  an  arbitrary  degree  of  certainty 
in  any  finite  time.  The  statement  becomes  reasonable  as  soon  as  one 
recognizes  that  the  assumptions  imply  equal  input  signal  to  noise  ratio 
in  every  frequency  band.  Thus  the  optimum  detector  can  in  principle 
extract  as  much  useful  information  by  observing  a  100  cps  band  at 
extremely  high  frequencies,  where  signal  and  noise  power  are  both 
negligible,  as  by  observing  a  comparable  band  near  zero  frequency  where 
both  signal  and  noise  power  are  substantial. 

This  example  illustrates  a  common  danger  in  the  use  of  optimal 
procedures:  Assumptions  that  appear  at  first  glance  quite  innocuous 
may  have  exceedingly  far-reaching  consequences.  In  important  practical 
problems  the  signal  and  noise  spectra  are  in  fact  closely  similar  over 
the  frequency  range  where  most  of  their  pother  is  concentrated.  At  high 
frequencies,  however,  the  signal  spectrum  inevitably  falls  off  faster 
than  the  noise  spectrum,  if  only  because  of  white  noise  locally  generated 
in  the  circuitry,  and  this  effect  becomes  crucial  when  one  works  with 
optimal  techniques. 
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In  order  to  obtain  realistic  comparisons  one  must  therefore  restrict 
the  optimum  detector  to  the  frequency  range  over  which  the  postulated 
spectra  may  be  expected  to  describe  the  actual  state  of  affairs  with 
reasonable  accuracy.  Specific  calculations  were  carried  out  under  the 
assumption  that  signal  and  noise  spectra  both  fall  off  with  the  second 
power  of  frequency  above  1500  cps.  When  the  frequency  range  processed  by 
the  optimum  detector  extends  to  5000  cps  its  figure  of  merit  then  has  a 
value  2.6  times  that  of  the  standard  detector.  When  the  processed  frequency 
range  extends  to  10,000  cps  the  corresponding  figure  is  3*8«  (Report  No.  3? 
pages  28-29), 

In  order  to  avoid  the  somewhat  arbitrary  assignment  of  an  upper 
frequency  limit,  a  second  appr  oach  was  also  taken.  By  using  a  model  of  the 
transmission  characteristics  cf  sea  water  proposed  by  Eckart  (Report  No.  3> 
p,  29,  footnote  2)  one  can  obtain  a  modification  of  the  high  frequency 
characteristics  of  the  signal  spectrum  that  removes  the  basic  difficulty 
and  leads  to  a  finite  result,  even  when  an  infinite  frequency  range  is 
processed.  The  ratio  of  performance  indices  of  the  optimum  and  standard 
detector  is  now  range  dependent  Eq.  (94 )J  but  varies  only  from  a  value 
somewhat  less  than  2  to  about  8  for  the  variations  of  range  for  which  the 
analytical  assumptions  are  likely  to  be  valid.  Thus  the  two  lines  of 
analysis  lead  to  quite  comparable  results,  and  one  comes  to  the  conclusion 
that  improvements  of  this  order  of  magnitude  cannot  be  exceeded  unless 
signal  and  noise  spectra  differ  drastically  from  the  assumed  form. 

The  foregoing  comments  should  not  be  construed  to  infer  that 
improvements  by  factors  of  2  to  8  are  unimportant.  If  they  were  realizable 
by  even  moderately  simple  instrumentations,  they  would  be  quite  significant. 
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However,  in  interpreting  these  figures  it  is  necessary  to  keep  in  mind 
two  factors:  a)  The  larger  improvement  figures  correspond  to  conditions 
such  as  extremely  close  range  or  identical  shape  of  signal  and  noise 
spectra  over  a  very  wide  frequency  band,  conditions  that  are  either  of 
limited  interest  or  unlikely  to  be  satisfied  in  practice «  b)  The 
optimal  instrumentation  i3  likely  to  be  very  complicated.  In  practice, 
therefore,  one  is  likely  to  pay  a  very  high  price  in  complexity  of 
instrumentation  for  only  limited  gain  in  performance, 

2)  Comparison  of  optimum  and  standard  detector  undar  assumption  B„ 

The  effort  to  exploit  the  presence  in  the  target  signal  of  low 
frequency  periodic  functions  and  amplitude  modulation,  by  one  of  these, 
of  the  high  frequency  range  has  taken  two  distinct  directions. 

a)  A  study  was  made  of  the  optimum  scheme  for  detecting  a  Gaussian 
signal  amplitude  modulated  by  a  sinusoid  in  the  usual  isotropic 
Gaussian  noise  background.  This  is  very  nearly  equivalent  to  an 
optimum  detector  processing  only  the  frequency  range  <  f  <  Wg 
of  assumption  B.  In  practice  the  frequency  would  be  sufficiently 
low  so  that  the  power  of  the  random  signal  component  in  0  <  f  < 
would  be  a  small  fraction  of  the  tooal  signal  power 0  However,  the 
periodic  low  frequency  components  of  the  signal,  which  could  well 
be  important  to  the  detection  process,  are  ignored.  The  reason 
for  considering  a  detection  process  neglecting  such  important 
information  is  that  the  frequencies  of  the  periodic  signal 
components  are  not  known  a  priori  and  that  the  self  noise 
generated  near  the  receiver  is  likely  to  contain  very  strong 
periodic  components  in  the  same  frequency  range  so  that  the 
filtering  problem  would  be  very  formidable. 
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b)  An  attempt  was  made  to  utilize  the  low  frequency  periodic  signal 
components  as  well  as  the  high  frequency  amplitude  modulation 
without  requiring  either  knowledge  of  the  frequencies  of  the 
periodic  components  or  identification  and  removal  of  the  periodic 
self  noise#  The  analyzed  system  is  not  optimal,  but  it  ha3  at 
least  seme  of  the  features  that  one  would  expect  to  find  in  an 
optimal  system# 

The  approach  described  under  a)  is  covered  in  detail  by  Progress 
Report  No.  $0  The  following  specific  cases  are  considered: 

i)  Noise  and  unmodulated  signal  both  have  the  properties  of  white 
noise  limited  to  the  same  band.  The  receiving  array  consists  of 
a  single  hydrophone c 

ii)  Signal  and  noise  have  the  same  properties  as  in  i)e  There  are 
K  hydrophones  spaced  sufficiently  far  apart  so  that  the  noise 
correlation  between  different  hydrophones  is  negligible, 

iii)  The  noise  and  array  properties  are  the  same  as  in  ii)  but  the 
unmodulated  signal  spectrum  falls  off  witn  the  second  power  of 
frequency  above  a  certain  point. 

iv)  The  unmodulated  signal  has  the  properties  of  i)  but  the  noise 
spectrum  falls  off  with  the  second  poxjer  of  frequency  above  a 
certain  point.  The  receiving  array  consists  of  a  single 
hydrophone  » 

In  cases  i),  ii),  and  iv)  the  ratio  of  the  figure  of  merit  of  the 
optimum  detector  to  that  of  an  optimum  detector  operating  on  the 

1  2 

unmodulated  signal  in  the  same  noise  background  is  approximately  1  +  ^  b  , 
F Progress  Report  No.  5,  Eqs.  (30),  (39)  and  (5U)1  „  The  symbol  b  is  the 
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modulation  index,  a  number  certainly  less  than  unity  and  in  most  practical 

m 

cases  substantially  less  than  unity.  Thus  only  a  very  limited  improvement 
in  detectability  results  from  use  of  the  knowledge  that  the  -signal  is 
periodically  amplitude  modulated.  Ir.  case  iii )  the  improvement  ratio  i? 
even  smaller  ^Eq.  (UU)  .  More  complicated  cases  have  not  been  analyzed, 
but  there  is  no  reason  to  expect  significantly  different  conclusions. 

The  result  appears  even  more  decisive  if  one  considers  that  the 
analysis  assumed  purely  sinusoidal  modulation  of  known  amplitude,  frequency 
and  phase.  In  practice  the  phase  and  amplitude  would  certainly  be  unknown 
and  the  frequency  would  most  likely  be  known  only  to  lie  between  certain 
limits.  The  effect  of  these  uncertainties  would  be  a  further  degradation 
in  the  detectability  index.  Thus  one  is  forced  to  tne  conclusion  that  the 
existence  of  periodic  amplitude  modulation  of  the  random  target  signal  is, 
by  itself  at  least,  of  little  value  in  target  detection,* 

The  reason  for  this  phenomenon  becomes  clearer  if  one  considers  the 
instrumentation  required  for  optimum  detection.  For  case  i)  it  is  given 
by  the  simple  configuration  shown  in  Fig.  1  of  Progress  Report  No.  5> 

(p,  11 )0  The  circuit  evidently  consists  of  a  standard  power  detector 
followed  by  an  arrangement  for  coherent  detection  of  the  envelope.  Thus 
the  basic  detection  scheme  is  incoherent  and  the  slight  gain  over  a 
simple  power  detector  results  only  from  the  secondary  coherent  detection 
operation  on  the  envelope.  The  latter  would  obviously  be  degraded  by  lack 
of  knowledge  of  envelope  amplitude,  phase  or  frequency,, 

Progress  Report  No.  7  deals  with  the  second  approach  (b)  to  the 
problem  of  utilizing  the  added  information  supplied  by  assumption  Br 
It  postulates  the  specific  instrumentation  shown  in  Fig„  1  of  the  report. 
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The  outputs  of  the  various  array  elements  are  added  and  the  resulting 
signal  is  split  into  the  two  frequency  ranges  0  <  f  <  and  <■  f 
The  high  frequency  range,  with  its  sinusoidal  amplitude  modulation,  is 
processed  by  a  standard  power  detector.  The  output  of  this  detector 
contains  the  modulating  signal  in  addition  to  fluctuating  components. 

It  is  used  to  multiply  the  signal  carried  by  the  low  frequency  channel, 
which  contains  a  sinusoidal  component  at  the  modulation  frequency. 

Thus  the  high  frequency  channel  is  used  to  generate  a  sinusoidal  "reference" 
for  the  coherent  demodulation  of  the  sinusoidal  low  frequency  component 
of  the  signal.  The  attractive  feature  of  this  scheme  is  that  the  low 
frequency  sinusoid  and  the  "reference"  automatically  assume  the  same 
frequency.  It  is  not  necessary  to  have  a  priori  knowledge  about  this 
frequency,  and  the  presence  of  sinusoidal  noise  components  in  the  low 
frequency  range  presents  only  a  secondary  problem,  as  long  as  they  do  not 
coincide  too  closely  with  the  signal  sinusoid.  Only  the  assumption  that 
low  frequency  sinusoid  and  modulation  envelope  are  in  phase  may  be  somewhat 
artificial.  The  weakness  of  the  proposed  scheme  is  the  low  amplitude  of 
the  reference  sinusoid,  particularly  when  the  modulation  index  is  small, 
and  the  presence  in  the  "reference"  channel  of  large  non-sinusoidal 
fluctuations,  all  of  which  are  "noise"  to  the  proposed  instrumentation. 
Detailed  analysis  shows  that  this  defect  outweighs  the  advantages  of  more 
complete  utilization  of  assumption  B.  Under  typical  operating  conditions 
the  performance  index  of  the  proposed  instrumentation  is  substantially 
lower  than  that  of  a  simple  power  detector  operating  on  the  high  frequency 
band  alone  |"  Eq,  (2U)j .  One  must,  of  course,  keep  in  mind  that  the  proposed 
detector  is  not  optimal.  One  could,  for  instance,  use  it  to  supplement  a 


10 


simple  power  detector  ana  thus  presumably  obtain  a  performance  somewhat 
better  than  that  of  the  power  (standard)  detector*  However,  the  resultc 
if  the  report  raise  some  doubts  whether,  even  with  optimal  procedures, 
detectability  under  assumptions  B  i3  significantly  greater  than  under 
assumptions  A,  For  the  present  this  question  must  still  be  regarded  as 
unresolved. 

The  conclusions  from  l)  and  2)  concerning  possible  improvements  over 
the  standard  detection  scheme  (Fig,  1  of  Progress  Report  No,  3)  may  be 
summarized  as  follows: 

a)  If  signal  and  noise  are  characterized  only  by  total  power  and 
spectral  properties  and  the  spectral  properties  do  not  differ 
drastically,  the  performance  index  of  the  standard  detector  is 
not  greatly  inferior  to  that  of  the  optimum  detector.  One  must 
also  keep  in  mind  that  the  instrumentation  required  to  realize 
the  optimum  detector  is  apt  to  be  complicated.  Unless  simple 
modifications  of  the  standard  detector  can  be  found  that  make  it 
optimal,  the  attainable  improvement  would  not  appear  to  justify 
the  increased  complexity  of  equipment  in  most  practical  cases, 

b)  Knowledge  that  the  random  signal  is  periodically  amplitude 
modulated  is  of  little  value  in  detection, 

c)  It  appears  that  large  improvements  in  performance  relative  to  the 
standard  detector  can  be  made  only  under  one  of  the  following  two 
conditions: 

i)  An  effective  scheme,  different  from  the  one  discussed  in 
Progress  Report  No,  7>  can  be  found  for  detecting  low 


11 


frequency  periodic  signal  components  in  the  presence  of 
strong  periodic  components  of  self  noise c  This  requires 
continuous  monitoring  and  careful  processing  of  self  noise 
data* 

ii )  It  is  possible  to  give  a  more  detailed  characterization  of 
the  signal  than  is  implied  by  assumptions  A  or  B.  In 
particular  it  would  be  valuable  to  have  a  description  of 
the  typical  waveshape  of  the  signal,  with  as  few  random 
parameters  as  possible*  This  would  open  the  way  to  use  of 
genuine  coherent  detection  schemes. 

3)  Effect  of  clipping  and  other  nonlinear  operations  on  target 

detectability,,  Assumptions  A»  Same  noise  power  at  each  hydrophone, 

The  procedure  of  clipping  the  basic  hydrophone  data  before  further 
processing  has  certain  practical  advantages,  particularly  when  digital  data 
handling  techniques  are  used.  Theoretical  studies  were  therefore  undertaken 
to  determine  the  effect  of  data  clipping  on  target  detectability  and  to  see 
whether  nonlinear  operations  more  general  than  clipping  might  lead  to  more 
advantageous  results.  The  problem  was  approached  by  two  rather  different 
routes* 

a)  Progress  Report  No.  U  considers  the  standard  (power)  detector  as 
the  basic  instrumentation  but  allows  an  arbitrary,  odd  function, 
zero  memory  nonlinear  device  to  be  inserted  at  the  output  of  each 
hydrophone.  (See  Fig.  1  of  Progress  Report  No.  1*.  Note  that  the 
linear  processor  and  the  clipper  are  special  examples  of  such  a 
device , )  The  analysis  is  carried  out  under  the  assumption  that 
signal  and  noise  spectra  are  identical  in  shape  and  that  noise 
correlation  between  different  hydrophones  is  negligible.  The 
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result  is  stated  in  terms  of  the  usual  output  signal  to  noise  ratio'*' 
and  also  in  terms  of  a  second  figure  of  merit,  defined  as  the 
difference  between  on  target  and  off  target  average  output  divided 
by  off  target  average  output.  In  each  case  ^Eq3a  ($8)  and  (U0)| 
it  is  demonstrated  that  the  best  performance  is  achieved  with  a 
linear  processor.  However,  the  dependence  of  the  figure  of  merit, 
particularly  the  output  signal  to  noise  ratio,  on  the  specific 
nature  of  the  processor  is  not  at  all  critical  if  the  processor 
does  not  contain  dead  band.  The  signal  to  noise  ratio  for  the 
clipped  instrumentation  falls  below  the  optimum  by  a  factor  of 
only  0,89.  (Fig,  3), 

Progress  Report  No,  U  also  considers  modifications  in  the 
basic  detector  characteristic,,  The  processors  are  now  taken  to 
be  linear  but  the  conventional  square  law  device  is  replaced  by  a 
general  even  power  device*  It  is  shown  that  maximum  output  signal 
to  noise  ratio"*-  occurs  when  the  detector  is  a  square  law  device, 
but  once  again  the  index  of  performance  varies  only  slowly  with 
modifications  in  the  detector  characteristic a  The  second  figure 
of  merit,  the  normalized  difference  of  on  and  off  target  average 
output,  can  exhibit  a  very  different  behavior.  If  the  detector 
has  a  dead  band,  for  instance,  the  performance  index  grows 
monotonically  with  the  size  of  the  dead  band  and  approaches 
infinity  in  the  limit  [ Eq,  (112  )L  The  reason  for  this  peculiar 


Note  a  slight  difference  in  definitions  between  Progress  Report  No.  U 
and  the  reports  discussed  earlier:  In  Report  No.  U  the  average  signal 
component  at  the  output  is  defined  as  the  difference  between  average  output 
on  and  off  target.  In  the  reports  mentioned  earlier  it  is  the  increase  in 
average  output  which  occurs  when  a  signal  appears  that  was  previously  absent, 
the  array  being  steered  "on  target"  at  all  times 0 
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phenomenon  is,  of  course,  that  the  off  target  average  output 
tends  to  zero  more  rapidly  than  the  on  target  average  output. 

This  observation  is  not  entirely  of  academic  interest,  for 
substantial  reductions  in  off  target  average  output  can  be 
achieved  for  values  of  dead  band  that  cause  only  minor  degradation 
in  output  signal  to  noise  ratio  (Figures  8  and  9).  Thus  the  "off 
target  plateau"  of  the  directivity  pattern  can  be  reduced 
materially  without  significant  sacrifice  of  output  signal  to  noise 
ratio. 

b)  The  second  approach  to  the  problem  of  nonlinear  data  processing  was 
an  attempt  to  assess  the  inherent  "cost  of  clipping"  by  comparing 
the  output  signal  to  noise  ratio  of  optimum  detectors  operating 
on  clipped  and  unclipped  data0^  This  analysis  is  contained  in 
Progress  Report  No.  6.  Results  were  obtained  only  under  the 
following  special  assumptions:  Signal  and  noise  both  have  the 
properties  of  Gaussian  white  noise  band-limited  to  0  <  f  <  W  cps. 
The  hydrophone  outputs  are  sampled  at  regular  intervals  of  ^ 
seconds  and  there  is  no  correlation  between  noise  outputs  of 
different  hydrophones.  Under  these  conditions  the  rms  output 
signal  to  noise  ratio  of  the  detector  operating  on  clippeu  samples 


1  -  i  times  as  large  as  that  of  the  detector  operating  on 
W 

unclipped  samples,  where  M  is  the  number  of  hydrophones  in  the 
array.  It  is  probably  realistic  to  consider  this  as  the  actual 
cost  (in  terms  of  detectability)  of  clipping  sampled  data0  One 
must  keep  in  mind,  however,  that  part  of  the  information  loss 


"wien  the  detection  scheme  is  optimal  it  is  only  necessary  to  consider 
information  destroying  nonlinear  operation  (such  as  clipping )  on  the  basic 
data.  A  one  to  one  nonlinearity  could  always  be  removed  by  an  inverse 
operation  in  the  optimal  detector. 
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night  be  due  to  sampling  rather  than  to  clipping.  For  while  the 

assumed  sampling  rate  is  sufficient  to  avoid  information  loss  from 

sampling  in  the  unclipped  case,  the  spectrum  of  the  clipped  signal 

is  not  bandlimited,  so  that  no  finite  sampling  rate  can  guarantee 

complete  re c ons true t ability  of  the  continuous  signal  from  the 

samples.  It  is  interesting  to  observe  that  evaluation  of  Eq,  ($8) 

of  Progress  Report  No,  U  for  the  white  noise  spectra  postulated 

in  Progress  Report  No,  6  leads  to  a  value  very  close  to  —  for  the 

degradation  in  output  signal  to  noise  ratio.  This  at  least 

suggests  -  although  it  certainly  does  not  prove  -  that  the  figure 
2 

~  is  indeed  a  good  estimate  of  the  "cost  of  clipping"  for 
reasonably  large  arrays 0  It  further  suggests  that  the  optimum 
detector  for  processing  clipped  data  derived  from  spectra  with 
the  white  noise  properties  assumed  in  Report  No,  6  does  not  differ 
greatly  from  the  standard  (power)  detector, 

1  )  Effect  of  clipping  on  target  detectability.  Assumptions  A, 

Noise  power  may  vary  from  hydrophone  to  hydrophone. 

One  of  the  primary  reasons  for  clipping  hydrophone  data  prior  to 
further  processing  is  the  possibility  that  the  noise  power  may  vary  sharply 
from  hydrophone  to  hydrophone.  Due  to  malfunction  or  special  operating 
conditions  one  hydrophone,  or  a  group  of  hydrophones,  may  contribute  far 
more  noise  than  the  rest.  When  this  is  a  serious  possibility,  it  appears 
reasonable  to  clip  all  hydrophone  outputs,  incur  the  small  loss  of 
information  discussed  in  3)j  but  prevent  a  large  contribution  to  the 
output  noise  from  a  few  faulty  hydrophones.  This  type  of  consideration 
is  in  fact  one  of  the  primary  reasons  for  the  use  of  clipping  in 
practical  instrumentations. 
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The  theoretical  background  of  this  problem  is  explored  in  Progress 


Report  No.  2*  The  basic  instrumentation  investigated  was  that  of  the 
standard  detector  with  and  without  clippers*  j^Fig.  1  The  analysis 
assumed  signal  and  noise  spectra  of  identical  shape  and  the  absence  of 
noise  correlation  between  different  hydrophone  outputs* 

Output  signal  to  noise  ratios  for  the  clipped  and  unclipped 
instrumentations  were  compared  under  several  sets  of  assumptions 
concerning  noise  power  at  different  hydrophones. 

a)  The  probability  density  of  noise  power  for  the  i*  1 
hydrophone  is  uniform  on  a  logarithmic  scale  over 


and  vanishes  elsewhere. 


Eq„  (1U) 


b)  The  probability  density  of  IL  is  uniform  on  a  linear  scale  over 

<  Njj  and  vanishes  elsewhere.  Eq.  (26)  j 

c )  The  probability  density  of  varies  as  over  <N^  <  Njj 


and  vanishes  elsewhere, 


Eq„  (29  )J 

d)  The  noise  power  at  each  hydrophone  can  assume  only  the  values 

and  Njj  and  does  so  with  probabilities  p  and  (l  -  p)  respectively, 

[Eq.  (15)] 


In  each  instance  the  noise  powers  at  different  hydrophones  were  assumed  to 
be  statistically  independent# 

Comparisons  were  made  primarily  on  the  basis  of  a  figure  of  merit 
defined  as  the  difference  between  on  target  and  off  target  average  output 
divided  by  off  target  average  output  (as  in  Report  No.  U).  The  output 
signal  to  noise  ratio  was  also  considered  and  found  to  have  identical 
properties  except  for  a  multiplying  constant. 
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The  results  show  the  expected  trend  in  each  case0  When  there  is  little 
or  no  variation  in  noise  power,  the  unclipped  instrumentation  is  superior 
by  a  small  margin.  When  large  variations  in  noise  power  can  occur,  the 
clipped  instrumentation  becomes  superiors  The  exact  extent  of  this 
superiority  depends  rather  critically  on  the  probability  distribution  of 
the  noise®  In  cases  b)  and  c)  the  ratio  of  clipped  to  unclipped 
performance  indices  is  close  to  unity  for  physically  reasonable  values  of 
and  N^0  In  case  a)  and  particularly  in  case  d)  (when  approximately 
half  of  the  hydrophones  are  in  the  high  noise  state)  very  substantial 
improvements  in  performance  are  made  possible  by  clipping 0 

The  results  of  sections  3)  and  U)  may  be  summarized  as  follows: 

a)  The  clipping  of  hydrophone  data  never  causes  a  large  decrease 

in  signal  detectability.  It  may  result  in  substantial  improvement 
if  large  variations  in  noise  power  from  hydrophone  to  hydrophone 
are  likely, 

b)  Signal  detectability  is  not  at  all  critically  dependent  on  the 
specific  properties  of  nonlinear  processors  operating  on  the 
basic  hydrophone  data  unless  the  processors  have  a  significant 
amount  of  dead  band.  Neither  is  the  precise  nature  of  the 
nonlinearity  used  in  detection  critically  important.  Since 
some  properties  of  the  bearing  response  pattern  are  more 
critically  dependent  on  these  parameters,  it  appears  possible 
to  improve  the  bearing  response  pattern  without  seriously 
affecting  detectability. 
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Progress  Report  No.  1  has  not  been  mentioned  in  this  discussion. 

It  deals  with  the  correlation  between  noise  output  at  different  hydrophones 
and  provides  computational  results  in  this  area  that  serve  as  background 
for  several  of  the  other  reports.  In  particular  it  shows  under  what 
conditions  one  can  make  the  very  convenient  assumption  that  the  noise 
correlation  between  different  hydrophone  outputs  is  negligible. 
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SPACE-TIME  CORRELATION  IN  ISOTROPIC  NOISE  FIELDS 


I.  Introduction 

Several  investigations  have  been  made  in  order  to  determine  the 

cross- correlation  between  signals  produced  by  two  omnidirectional 
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transducers  located  in  a  field  of  many  scattered  noise  sources.  *  * 

Knowledge  of  the  noise  cross- correlation  is  often  necessary  in 
determining  the  performance  capabilities  of  arrays  of  such  transducers 
which  are  used  to  determine  the  direction  of  a  signal  source  located 
in  the  field  of  the  scattered  noise  sources. 

In  this  report,  an  attempt  is  made  to  expand  and  simplify  the 
interpretation  of  the  results  given  in  the  references  listed.  Three- 
dimensional  and  two-dimensional  isotropic  noise  fields  will  be 
considered.  The  field  is  considered  isotropic  in  that  the  power  received 
from  any  direction  is  constant. 

An  infinite  number  of  noise  sources  distributed  uniformly  on  the 
surface  of  a  sphere  of  infinite  radius  from  the  local  transducers  is 
used  to  simulate  the  three-dimensional  isotropic  field. 

For  the  two-dimensional  field,  an  infinite  number  of  noise  sources 
are  distributed  uniformly  along  a  circle  of  infinite  radius  in  the  plane 
of  signal  propagation. 

In  addition,  the  noise  produced  by  each  source  is  assumed  to  be 
stationary  and  statistically  independent  of  that  produced  by  every 
other  source,  and  each  noise  source  is  described  by  the  same  autocorrelation 
function. 


Al-1 


II.  Cross-Correlation  with  Time  Delay- 


Figure  1  describes  the  geometry  for  the  spherical  noise  field. 

A  z 
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It  has  been  shewn  that  the  normalized  cross- correlation  function  for  the 
delayed  output  of  transducer  1  with  the  output  of  transducer  2  is  given 
by  Equation  (l). 
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In  Equation  (l),  -  is  the  time  required  for  the  signal  to  travel  the 
0 

!  N 

distance  between  the  two  transducers,  and  p  (t)  is  the  normalized 
autocorrelation  function  for  each  noise  source.  A  chaige  of  variables 
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results  in  Equation  (3 ) 


Al-2 


Si* 


P<V  I  P  (^  )  dr 

3 


-X  +X 
3 


Equation  (3)  simply  represents  the  average  value  of  p^x*)  over  an 

t 

interval  2x  in  length,  centered  at  x  ■  x.  Figure  2  depicts  this 

3 

interpretation  graphically. 
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Fig.  2 

Graphical  Interpretation 
of  p(t  ,  t) 
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For  the  circular  noise  field,  the  cross-correlation  function  has 


been  shewn  to  beJ 
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With  the  variable  transformation  in  Equation  (2),  Equation  (li)  becomes 
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Equation  (5)  represents  a  weighted  average  of  p  (x)  over  an  interval 

i 

2x  centered  at  x  =  x.  The  weighting  function  is  shown  in  Figure  3. 
s 

It  can  be  seen  that  the  values  of  p^T*)  at  the  ends  of  the  interval 
are  weighted  very  heavily. 
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Note  also  that  the  average  value  of  the  weighting  function  expressed  in 
Equation  (6)  is  unity. 
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A  comparison  of  the  results  in  Equations  (3)  and  (5)  reveals  that  the 
evaluation  of  the  cross- correlation  function  for  the  three-dimensional 
noise  field  is  considerably  simpler  than  that  for  the  two-dimensional 
case  For  each  pair  of  transducers  in  an  electronically  steered  array, 
the  observation  is  also  made  that  the  maximum  relative  electrical  delay 
necessary  to  steer  the  directional  response  of  the  array  is  never  greater 
than  the  wave  propagation  time  between  the  pair  of  transducers  in  question. 
In  terms  of  the  symbols  used  previously,  the  following  is  always  true 
for  the  determination  of  bearing  response  patterns. 

T  <  T  (7) 
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An  examination  of  Figures  2  and  3  indicates  that  the  cross-correlation 
function  for  a  two-dimensional  field  is  less  than  that  for  a  three- 
dimensional  field  for  values  of  t  near  zero.  Identical  well-damped  auto¬ 
correlation  functions  are  assumed  in  both  cases.  In  the  two-dimensional 

t  i .  t  + 

case,  the  small  values  of  p  (t  J  near  x  ■  -  v  are  weighted  more 

s 

heavily  than  the  same  values  in  the  three-dimensional  case.  However, 
for  values  of  x  near  t  ,  p(x  ,  x)  for  the  two-dimensional  field  should  be 
greater  than  that  for  the  three-dimensional  field  since  the  large  value 

i .  t .  i 

of  p  (x  )  near  x  =0  receives  infinite  weighting  in  the  former  case. 

Evaluation  of  p(x  ,  x)  for  two  simple  forms  of  p  (x  )  follows, 
s 

III.  Evaluation  of  Cross-Correlation  Function  -  Case  A 


The  function  p(x  ,  x)  is  first  evaluated  for  the  simple  exponential 
s 

autocorrelation  function  given  in  Equation  (8) 

P  (x  )  =  e 


The  corresponding  spectral  density  is 
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The  integral  for  p(x  ,  x)  must  be  performed  in  two  parts  because  of 


the  term  x  appearing  in  the  expression  for  p  (x  ).  For  the  three¬ 


dimensional  noise  field: 
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If  t  °  t  ,  the  maximum  electrical  delay  necessary.  Equation  (10) 
s 


reduces  to 


p(v  Ts)  ■ 
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Furthermore,  if  cor  »  l  and  t  =  0, 
1  s  ’ 


p(T  o)  v  -L 


(12) 


CO  T 
O  S 


Equation  ( 12 )  also  places  an  upper  bound  on  p(x  ,  t)  for  any  t  and  t  . 

s  s 

If  the  separation  is  defined  in  terms  of  wavelengths  of  the  half¬ 
power  frequency  of  the  individual  noise  spectra,  the  relation  with  the 


parameter  co  t  is 
0  s 


The  upper  bound  for  p(x  ,x)  is 

s 


CO  T 

o  s 
2n 


(13) 


2n  /transducer  separation  in  wavelengths \ 
'  of  upper  half  power  frequency 


(lU 


For  the  two-climensional  noise  field,  Equation  (I4)  may  be  used  with 


several  approximations  to  evaluate  p(x  ,  x).  Equation  (l5)  is  obtained 


after  appropriate  substitutions  are  made. 

n 

1 

p(% ,  t)  -  i 


co  |x+x  cos  Q| 

6  °  3  d9 


0 


PQ  n 

-co  (xi'r  cos  9)  r  co  (x+x  cos  P) 
e°  s  dP+/e°  s  dP 


'9 


(15) 
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where 


(16) 


If  co  x  is  large,  the  integrand  has  negligible  value  except  in  the 
o  s 

vicinity  of  8  ■  Qq.  The  exponent  of  the  integrand  nay  be  approximated 
by  the  first  three  terms  in  the  Taylor  series  expansion  around  0  -  9qj 


co  (t  +  x  cos  9)  ^  0  -  co  x  sin  9(9  -  9  )  -  co  x  cos  9  l9. 

0  S  0  8  O  0  O  S  0  2 


•  0  -  co~\/x  2  -  t2  (0-9  )  +  co  X 
o\J  s  ~ 


(6  -  O2 


O  O  2 


(17) 


If  a  =  co  2  -  x2  and  b  *  coqt,  then  Equation  (15)  becomes 

P(Ts»  t)  "  I 


9°  a(9-9  i  )  -  £(0-©)2  f  -a(9-0  )  +  |(9-9)2 
e  0  c  0  d@  +  I  e  °  ^  de 


9 


(18) 


Befera  Equation  (18)  is  evaluated  for  the  general  case,  0  <  x  <  x  , 

9 

several  special  cases  will  be  evaluated.  For  x  =  0,  b  =  0,  a  K  “0TS>  and 


©o  =  tj.  Then  Equation  (18)  becomes 


r2  »t(0-|) 

P(^s,  0)  I  e  0  s  2  d9 


n 


CO  T 
0  S 


(19) 


For  T  =  T.a  =  01b  =  coT,  and  9  =  n.  The  seoond  integral  in 

S'  '  OS'  o 

Equation  (18)  disappears,  and  the  result  is 

n 


,  f  -l  co  x(9  -  a) 2 

c  de 

Jr. 


~\l  2n  co  x 
v  os 


(20) 
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The  results  in  Equations  (19)  and  (20)  are  valid  only  if  co  t  >  >  1. 

o  s 

For  values  of  t  near  zero,  such  that  the  quadratic  portion  of  the 
exponent  may  be  neglected,  the  result  in  Equation  (2l)  applies. 


Finally,  the  integrals  in  Equation  (18)  are  evalur.ted  by  completing 

i 


the  squares  in  the  exponents,  and  by  letting  9 

\2 


9  -  9  ,  and  9  =9-9. 

o'  o 


p(tr,  t)  =  l 


a2  f*eo  (^9*  -  a/Vb1)^  a^ 

.*/  -  d9;+e"^ 


o  (4%  -a/lfi?) 
- f—L 


2  1 
ii 


d9 


(22) 


Additional  variable  substitutions  found  in  Equations  (23)  and  (21+ )  and  the 
definition  of  tabulated  functions^*'’  in  Equations  (25)  and  (26)  allow  the 
final  result  to  be  written  in  Equation  (27). 


i2 


(23) 


y  * 


erf  (a) 


(2b) 

(25) 
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(26) 


0 

(  t2 

F(3)  =  I  e  '  dt 


f 


L  2  r  iV  » 

2b  ■'(F  F  I  — |  -  f|  -4^ -Jf(n  -  9  )]  >  ( 

IVzBI  Na  °l 

u  “V 

The  result  in  Equation  (10)  for  the  three-dimensional  noise  field  is 

plotted  in  Figure  U  for  two  values  of  ai_x  as  a  function  of  t/t  and 

s  s 

compared  to  similar  results  in  Equations  (19),  (20),  (21),  and  (27)  for 
the  two-dimensional  field.  Previous  quantitative  conclusions  ooncerning 
the  relative  size  of  p(x  ,  t)  for  the  two  types  of  noise  fields  are  seen 
to  be  substantiated  in  Figure  k. 

IV.  Evaluation  of  Cross-Correlation  Function  -  Case  B 

The  normalized  autocorrelation  function  found  in  Equation  (28)  is 


considered  next. 


pV)-?  I1  vr^l  e 


(-  0  +  V  <*.  2  -  T>o  |t  * 


i  1  +  ^ 

1  vFTi 


(-S-Vb2  -  i)“  t1 


If  <1,  Equation  (28)  nay  be  put  into  the  form 


1,  I,  1  “^olT 
p  (T  )  =  6 

V  1  - f 


COS  (Vl  -  j*  £0  T 


'I  +  COS  ^*^1  -  *  ^  )  (29) 


The  corresponding  spectral  density  for  p  (t  )  is  the  simple  narrow  band 
spectrum  which  follows. 


Al-9 


In  Equation  (30),  £0q  is  the  center  frequency,  end  2^0^  is  the  bandwidth 
between  half  power  points. 

In  order  to  find  p(r  ,  t)  for  the  three-dimensional  field.  Equation  (28) 

s 

is  substituted  into  Equation  (3).  One  step  in  the  integration  process 


yields 

p(V 


Evaluation  a  ■  the  limits  indicate^  in  Equation  (3l)  and  simplification  of 
terms  results  in  Equation  (32). 


Combination  of  the  terms  inside  the  brackets  of  Equation  (32)  yields 
Equation  (33) 5 
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—  ^60  X 
0  n 


p(t  ,  t)  -  -  (cosh2  *  co^t  -  sin2Vi  -<J2  w  t) 

Vl  -  '  Vs 


1/2 


x  sin 


'^1  -  ‘  2  co  t  —  tan-1  (ten  Vi  -  $>2' co  t  tanh  £  co  t) 

OS  0  0  •)  O 

(33) 


Special  cases  for  t  ■  0,  and  T  •  x  are  found  in  Equations  (3U)  end 

s 

(35)  respectively. 


-  rf  co  t  sin  -  z 
P(T„  0)  -  e  '  03 


CO  T 
0  S 


■\A-r 

-2  T  co  t  sin  2 

p<v  V  ■ 6  • 


CO  T 
O  S 


vi-T  Vs 


2Vl  -  r  T  »  T 

T  A  OS 


(3U) 


(35) 


Interesting  behavior  of  p(t  ,  t)  is  also  apparent,  if  f  is  small, 

s 


for  values  of  t  in  Equation  (36) 
s 


^  ^  x  nn 

T  83  T  +  *a*  p  ** 

a  2co  a 

0 


n  =  0,1, 2, 3,. . . .  (36) 
a  =Vl  -  if 2  -  1 


With  the  substitution  of  the  relation  in  Equation  (36),  Equation  (33) 
becomes 

-Jot  n*S 

nn  *  e  3  °  £  ,  ,  2,.  .  .2  ,V2 


p(x  +  JS.'-  ,  t) 

2aco  *  7 


aco  t  + 


oT  T  T 


( cosh^  co  t  -  sin  aco  t) 
nn  o  o  o 


x  sin 


aco 


t  +  -  tan”x(tan  aco  t  tanh  Vco  t) 

id  o  o 


-1 


(37) 


Examination  of  the  arctan  function  in  Equation  (37)  reveals  that 

nn 


tan” ‘'"(tan  ocoqt  tanh^oo^r)^ 


aco  t  when  aco  t  n 
o  o  2 


■  aco  t  when  ^co  t  >  1 
o  o 


(38) 
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I 


Thus 


p(T  + 


nn 


,  *0  * 


.  CO  T  - 

oT  e  7a 

—  (cosh  j  CO  T 


aoo  t  + 
o 


nn 


-  sin 


aro  t) 
o 


V2 


sin 


nn 


(39) 


Sufficient  computation  rrith  Equation  (32)  or  (33)  would  result  in  a 

contour  map  of  p(t  ,  t),  similar  to  that  in  Reference  2.  However, 
s 

Equation  (39)  sives  some  important  details  about  such  a  map,  without 
resorting  to  detailed  computation  procedures.  The  loci  described  by 
Equation  (36)  are  plotted  in  Figure  5.  The  axes  are  appropriately 
normalized. 


The  lines  in  Figure  £  for  which  n  is  even  give  approximate  loci  for 
which  p(t  ,  t)  is  zero.  The  lines  for  which  n  n  1,  9,  etc.  are  the 

approximate  loci  for  the  ildges  of  the  contour  plot  in  regions  where 
p(T  ,  t)  is  positive.  Similarly,  the  lines  for  which  n  ■  3*  7*  Uj  etc., 
describe  the  approximate  loci  for  the  ridges  for  negative  values  of 
P(t  ,  t).  The  relations  are  exact  where  cuo  t  = 

Finally,  from  Equation  (33)  it  can  be  seen  that  the  upper  bov  1  on 
the  magnitude  of  p(x  ,  t)  for  all  values  of  T  and  T  is 

3  S 


lp(T,  *Ol  5? 


-  o  u  T 
OS 


Vl  “  W  T 

v  0  OS 


cosh  <fCO  T  <  —  -y - 

•  0  •Vi-'2  0.x 


T  <  T 


0  S 


3 

(Uo) 


For  the  two-dimensional  noise  field,  substitution  of  Equation  (28) 

into  Equation  (h)  yields  an  integral  which  does  not  yield  to  the 

approximation  techniques  tried  in  Section  III.  Due  to  the  oscillating 

nature  of  p^x*),  the  integrand  can  have  appreciable  value  over  the  entire 

✓ 

interval  of  integration  for  small  * .  The  Taylor  series  approximation  in 
Equation  (l?)  for  the  exponent  is  not  useful  since  this  approximation 
becomes  inaccurate  for  values  of  I 0  —  ©  I  >  1.  Such  angular  differences 
are  well  within  the  interval  of  integration,  since  I  6  -  i  can  be  as 
large  as  n. 

Numerical  integration  seems  to  offer  the  most  fruitful  approach  to 

the  solution  of  p(x  ,  t)  for  narrow-band  two-dimensional  noise  fields. 

s 

This  will  not  be  attempted  here. 


V.  Evaluation  of  Cross-Correlation  Function  -  Case  C 

In  order  to  gain  insight  into  the  question  whether  the  results 
obtained  thus  far  depend  critically  on  the  somewhat  arbitrarily  postulated 
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Tom  oi'  p(v),  a  second  narrowband  node!  will  bo  considered.  For 
convenience  in  mathematical  manipulation  the  normalized  spectral  density 


of  each  elementary  noise  source  is  represented  by  the  equation 


S(co) 


I  co  I 


co  ~\l2n  o 
o  v 


(co-co  )‘ 


+  e 


(co+coq) 

~2^~ 


2-, 


co  >  >  a 


Oil) 


Equation  (Ul)  is  not  Wiener-Hopf  factorable  and  therefore  represents 
only  a  limiting  form  of  a  realistic  spectral  function.  For  the  purposes 
of  the  present  discussion  this  causes  no  difficulty. 

The  autocorrelation  function  corresponding  to  Equation  (III)  is 


2  2 
t  a 


p(t)  = 


cos 


CO  T  - 
0 


TO" 


a 

co 


sin 


co  T 
0 


(U2) 


With  co^  >  >  a  the  second  tern  is  small  compared  to  the  first  for  the 
significant  range  of  t. 

The  desired  cross- correlation  for  the  three-dimensional  field  is  most 
easily  evaluated  by  use  of  the  expression 


P(V  *0 


S(co) 


sin  cote 
arc 


dco 


U3)1 


Substitution  of  Equation  (hi)  into  Equation  (U3 )  leads  to  the  result 


p(T  t)  =  — 


2  CO  T 
0  s 


a  /  s2  a  , 

-  T  (ts-r)  -  T 

>  sin  co  (t  -  t)  +  e 

0  v  s 


sin  co  (t  +t) 
o  s  1 

Oil*) 


This  relation  is  equivalent  to  Equation  (32)  in  case  B.  The  special  cases 


t  «  0  and  t  =  t  yield 
s 
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p(Tg,  0)  “  e 


1  a 

x  -  CO  T 

2\  %  0  s 


sin  co  t 
o  s 

co  T 

o  s 


(W) 


p(ts<  V  "  6 


-2  -  CO  T  I  .  o 

co  os  sin  2  co  t 
0  o  s 

2  co  x 
o  s 


0i6) 


If  one  identifies  — -  with  and  considers  ^  <<1  Equations  (1*5)  and 

o 

(L*6)  are  closely  analogous  to  Equations  (3U)  and  (35)  respectively. 

An  upper  bound  equivalent  to  Equation  (1*0)  can  be  obtained  by 

replacing  the  sinusoidal  functions  by  unity.  Thus 

t2  9.  .o  1  o2  2, 


|p(V  T>|  5  ~ 


2co  T 
0  s 


~\~^7  Wo2(ts^>2  ■  \  — 2  “o2(ts^)2 

d  (ft  0  3  d  &  d  0  S 


+  6 


CO  T 

0  s 


(U7) 


Ch  the  other  hand  Equation  (1*0)  can  be  rewritten  in  the  form 


-  w  x 
o  s 


V  l  Vs 


cosh  H  co  t 
A  o 


2‘\/l  -  \  '  co  T 


0  S 


-  a  co  (t  ~*r)  -^co(t+t) 

]  0K  3  '  .  J  0X  S 

e  +  e 


(ua) 


X  <  T 


Again  identifying  ^  with  ~  and  letting  «  1,  one  observes  that  Equation 

o 

(1*8)  is  the  direct  analogue  of  Equation  (1*7). 

Finally  it  is  clear  that,  except  in  a  neighborhood  of  the  origin 

sufficiently  small  so  that  ct(t  +  t  )  £  0(l),  the  function  p(t  ,  t)  of 

s  s 

Equation  (1*1*)  is  characterized  by  ridges  and  valleys  parallel  to  the  line 

x  =  x  (for  x,  x  >0),  much  as  those  shown  in  Figure  5.  It  therefore 
s  s 
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appears  that  the  general  properties  of  p(x  ,  x)  are  relatively  independent 

s 

of  the  precise  nature  of  the  autocorrelation  function  assumed  for  the 
elementaiy  noise  sources. 

VT.  Conclusions 

The  results  found  in  Equations  (12)  and  ( UO )  are  valuable  in  that 

they  allow  a  quick  estimate  of  the  maximum  value  of  p(x  ,x)  for  a 

s 

given  separation  of  transducers.  The  results  for  the  wide-band  case(A) 

and  the  narrow-band  case  (B)  are  approximately  the  same  when  is 

interpreted  to  be  the  upper  half  power  frequency  in  the  first  case,  and 

the  center  frequency  in  the  latter  case.  These  answers  in  an  approximate 

sense  can  be  e:±ended  to  other  spectra,  wide  or  narrow,  if  co  is 

interpreted  to  be  the  upper  half  power  frequency  for  the  spectral 

density  of  each  tiny  noise  generator  in  the  field. 

Also  a  general  pattern  of  regularity  has  been  established  for 

p(x  ,  t)  for  the  narrow-band,  three-dimensional  noise  field.  Coupled 
s 

with  Equation  (UO),  Figure  5  can  provide  a  more  accurate  estimate  of 

p(T  ,  t).  Figure  5  is  also  of  value  in  determining  transducer  spacings 
s 

that  yield  little  or  no  cross-correlation. 

A  comparison  of  results  for  p(x  ,  x)  for  the  three-dimensional  noise 

s 

field  and  for  the  two-dimensional  noise  field  is  also  interesting.  Under 

the  assumption  that  the  spectra  of  the  individual  noise  generators  are 

given  by  Equation  (9)  in  both  cases.,  p(x  ,  x)  decreases  for  increasing  x 

and  constant  x  for  the  three-dimensional  case,  and  increases  for  the 
s  3 

same  conditions  in  the  two-dimensional  case  if  «cx^  is  large  enough. 
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SIGNAL  DETECTION  BY  ARRAYS  IN  NOISE  FIELDS 
WITH  LOCAL  VARIATIONS 

I.  Introduction 

A  considerable  amount  of  work  has  been  done  in  recent  years^'^5^ 

in  analyzing  systems  composed  of  local  arrays  of  transducers  used  .to  detect 

the  presence  of  a  plane  wave  signal  in  an  isotropic  noise  background. 

Such  systems  have  had  direct  application  in  undersea  signal  detection, 

and  of  particular  interest  in  this  field  is  the  DIMUS  system.  The  DIMUS 

2 

system  has  been  defined  to  consist  of  the  following  steps  in  signal 
processing : 

1.  Sampling  the  output  of  each  transducer  at  small,  equal  time 
intervals; 

2.  Providing  infinite  clipping  for  each  sampled  output; 

3.  Delaying  the  clipped  signals  electrically  by  integral  multiples 
of  the  sampling  interval; 

U.  Adding  the  delayed  signals; 

3.  Squaring  the  sum; 

6.  Filtering  the  squared  sum. 

3  2 

Analyses  carried  out  by  Faran  and  Hills,'  and  Rudnick  have  assumed 
that  the  background  noise  power  output  of  each  transducer  is  the  same. 
Furthermore,  the  average  noise  power  output  on  a  short-time  basis  was 
assumed  to  be  constant  for  each  transducer. 

Although  the  system  analyzed  by  Thomas  and  Williams^  was  not  a  DIMUS 
system,  infinite  clipping  of  signal  plus  noise  was  incorporated  in  the 
analysis.  Thomas  and  Williams  compared  the  performance  of  a  system  with 
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infinite  clippers  to  a  similar  one  without  infinite  clippers  for  the 
situation  in  which  the  short-terra  average  noise  power  output  of  each 
transducer  varies  with  time.  However,  this  analysis  is  limited  by  the 
assumption  that  the  short-term  noise  power  outputs  of  every  transducer 
are  described  by  the  same  time  function. 

Experimental  observations  have  shown  that  the  average  noise  power 
outputs  of  transducers  mounted  in  an  underwater  array  are  not  the  same 
from  transducer  to  transducer.  One  possible  explanation  lies  in  the 
fact  that  different  turbulent  flow  patterns  are  set  up  around  each 
transducer  when  the  array  is  in  motion.  Also,  the  self-noise  due  to 
the  internal  construction  and  possible  malfunctions  is  different  for  each 
transducer. 

The  purpose  of  this  report  is  to  incorporate  the  variation  in  noise 
power  output  from  transducer  to  transducer  in  the  analysis  of  the 
performance  of  the  array.  A  system  which  is  essentially  DIMUS  in  nature 
is  analyzed  and  the  performance  is  compared  to  a  system  which  is  identical 
in  ever  respect  except  that  the  infinite  clippers  are  absent.  For  each 
of  the  two  systems,  two  quantities  form  the  bases  of  comparison: 

1.  The  maximum  average  signal  output  of  the  system,  which  occurs 
when  the  array  "beam"  is  steered  directly  at  the  signal  source. 

2.  The  ratio  of  the  square  of  the  maximum  average  signal  output 
to  the  variance  of  the  output  of  the  low-pass  filter. 

For  purposes  of  comparison,  two  performance  indices  are  calculated  by 

forming  the  ratios  of  the  quantities  defined  above  for  the  system  with 

clippers  present  and  for  the  system  with  clippers  absent.  The  assumptions 

involved  in  the  analysis  are  listed  in  the  following  section. 
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II.  Assumptions  and  Definitions 


The  system  to  be  analyzed  is  shown  in  Figure  1. 


Transducers  Infinite  clippers 

(t)+n, (t)| 


o- 


s2(t j+n2(t) 


Sj-(t)+nM(t) 


Summer 


Square  Low  Pass 

Law  Filter 


Fig.  1  Essential  Configuration  of  a  DIIIUS  System 


1.  The  transducers  in  the  array  depicted  in  Figure  1  have  omni¬ 
directional  characteristics. 

2.  The  normalized  autc correlation  functions  of  both  the  signal  and 
the  noise  inputs  to  the  clippers  in  Figure  1  are  identical. 

3.  The  crosscorrelation  between  noise  inputs  for  different  channels 
is  zero.  In  an  actual  system,  the  cross correlation  may  be  made 
arbitrarily  small.  This  assumption  makes  the  following  analysis 
mathematically  tractable. 

b.  The  signal  and  the  noise  inputs  to  the  clippers  are  gaussian  and 
stationary.  Signal  and  noise  are  independent. 

3>.  The  signal  power  input  to  each  clipper  is  the  same,  but  the  noise 
power  input  is  allov.  ^d  to  vary  from  channel  to  channel. 

6.  It  is  assumed  that  the  short-term  averages  for  the  input  signal 
powers  and  input  noise  powers  do  not  vary  significantly  for  any 
particular  channel. 
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7.  The  effect  of  sampling  is  neglected.  Effectively,  the  sampling 
takes  place  at  a  high  enough  rate  that  the  higher  harmonic 


power  spectra  produced  by  sampling  do  not  overlap  the  fundamental 
power  spectral  density. 

8.  The  digital  time  delays  have  been  omitted,  since  the  analysis 

is  primarily  concerned  with  "on  target"  and  "off  target"  quantities. 

It  is  assumed  that  such  combinations  of  delays  are  available  to 
electrically  steer  the  array  directly  at  the  target,  and  to  steer 
the  array  far  enough  away  so  that  gain  effects  due  to  the  array 
configuration  are  negligible. 

III.  General  Form  of  Systen  Directivity  Patterns 

The  correlation  matrix  for  the  noise  components  is  given  by  Eq.  (l). 

It  is  assumed  that  n.  (t )  =  sTTtT  =  0. 

i  l 

1L  0  0-0 

0  N2  o  —  0 

0  0  N3  —  0 

0  0  0  —  nm 

M 

The  correlation  matrix  for  the  signal  components  of  the  array 
inputs  is 

p(^)  p(t12+  p(TlM+  T 

p(t2i+  tO  pC^)  .  — '  p(T2M+  T 

p(TI.a+  P(tM2+ 

(2) 
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In  Eqs.  (l)  and  (2),  the  are  the  noise  power  inputs  to  each  channel, 
and  S  is  the  common  signal  power  input  to  each  channel.  The  function 
p(t)  is  the  normalized  autocorrelation  function  for  the  signal  and  the 
noise.  The  delay  factors  in  Eq.  (2)  represent  the  combined  effects 
of  spatial  time  delays  introduced  by  the  plane  wave  signal  relative  to 
each  pair  of  transducers,  and  also  the  electrical  time  delays  which  are 
usually  introduced  in  all  steered  systems,  but  which  are  not  included 
in  Figure  1. 

The  average  output  of  the  squarer,  y,  yields  the  so-called  directivity 

pattern,  which  may  be  obtained  by  varying  the  electrical  time  delays  and 

keeping  the  physical  orientation  of  the  array  fixed,  or  try  keeping  the 

electrical  time  delays  fixed  and  varying  the  physical  orientation  of  the 

array  relative  to  the  plane  wave  signal.  For  the  system  in  Figure  1, 

2  3 

with  the  clippers  removed,  the  average  output  has  been  shown  *  to  be 
the  sum  of  all  the  matrix  elements  in  Eqs.  (l)  and  (2)  with  t  =  0. 

The  subscript  A  in  Eq.  (3 )  and  following  work  is  used  to  denote  the 
fact  that  the  clippers  are  removed. 


11  M 


^97 


Z  Z  kj +  s  p  (v 


i=l  j=l 
M 


M  M 


E  *i*s<*»[  £  p 

i-1  j-i+l 


(3) 


i=l 


The  angle  Q  is  defined  as  the  steered  beam  angle  aid 
function  of  the  electrical  time  delays  which  make  up 


is  a  complicated 
part  of  the  t^.. 
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Figure  2  shows  a  typical  directivity  pattern  resulting  from  Eq.  (3).  The 
signal  and  noise  are  assumed  to  have  low -pass  spectra.' 


Fig.  2  Typical  Diiectivity  Pattern 


At  some  angle,  9  ,  the  steered  beam  is  "on  target"  and  a  maximum 
output,  results.  When  the  beam  angle  is  greatly  different  from  the 

target  angle,  the  directivity  pattern  approaches  a  lower  asymptote, 
yAL,  very  closely,  because  the  p(t^_.)  approach  zero.  These  quantities 


E  Mi 


1 1 


+  SM 


With  the  infinite  clippers  present  in  the  system,  the  directivity 
pattern  has  much  the  same  shape  as  that  in  Fig.  2.  However,  in  order 
to  define  y  (9)  mathematically,  it  is  first  necessary  to  define  a 
normalized  correlation  function,  jj,.  (t),  for  signal  and  noise: 


Furthermore, 


p(T  +  t  )  i  /  j 


^j(T)  m\ 


Vs  ♦  v^7r?75  ^ 


(7) 


Pfr) 


i  =  d 


The  correlation  functions  for  the  clipper  responses,  with  the  assumption 
that  the  clipper  output  is  either  +1  or  -1,  have  been  shown^  to  be 


.  f 


t-L .  .(t)  =  v.  (t)  v.(t  +  t)  =  -  arcsin^  (i.  .(t) 


ij  i  j 


(8) 


Note  that 


f 


^<0>  -< 


sis  l 

-  arcsin  < - ra - tto  p(T«  J  f  i  /  0 

[(S  +  +  N  )1/2  1J  J  (9) 


i  =  0 


The  directivity  pattern  is  found  by  summing  Eq.  (9)  over  all  i  and  j: 

J*  M 

j®5  L  L  I  arcstoin.^o) 

i=l  ;j=l  L  J 


M  M 

C-1  C"1  2 

M  +  2  )  )  —  arcsin 

O  n 

i=l  j=i+l 


(S  +  Ni)ly/2(S  +  N,.)^2  ^ 


(10) 


The  "on  target"  and  "off  target"  values,  y0  and  y.  respectively,  are 


M  H 

r*1  r->  o 


yn  c  M  +  2  )  )  ^  arcsin  <, - - rrw 

°  ft  *  (S  +  N.)1/2(s  +  N.)1/2 

1=1  j=l+l  «-  l  3  J 


z  M 


(11) 


=  M 


(12) 
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IV.  Performance  Comparison  -  Ratio  of  Average  Signal  Outputs 

The  average  signal  output  is  defined  in  Eq.  (13 )  to  be  the  difference 
between  the  "on  target"  value  of  the  directivity  pattern  and  the  "off 
target"  value  for  both  systems.  Note  that  the  signal  source  is  still 
present  for  both  conditions. 


Ay  =  y0  -  yL 


(13) 


For  purposes  of  comparison,  the  directivity  patterns  for  both  systems  are 
normalized  with  respect  to  the  "off  target"  response  in  each  case.  The 
ratio  of  the  normalized  Ay  in  both  cases  then  becomes  meaningful  as  a 
performance  index. 

The  indicated  computation  nay  be  carried  out  for  any  particular 
system  for  which  the  are  known  by  measurement.  However,  in  order  to 
get  some  general  results,  it  is  convenient  to  assume  a  statistical 
distribution  for  the  NL  which  applies  to  all  transducers.  Of  course, 
physical  measurements  should  support  any  assumed  distribution.  Two  simple 
statistical  distributions  are  assumed  and  analyzed  in  detail,  one  of  which 
is  continuous,  the  other  is  discrete.  The  continuous  distribution  of  N.  • 
is  described  by  the  probability  density  functioh,  f(N. ). 


f(N±) 


(N.  In  Nj/NjT1 

y 

o 


N.  5 
i 


nh 


(1U) 


elsewhere 


The  discrete  distribution  is  two-valued  and  is  described  in  Eq.  (15). 


nh1 

r 

P 

1 

j*  with  probability  < 

> 

l  nl  J 

l 

(1  -  p) 

(13) 
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Statistical  averages  of  Ay  over  the  arc  used  in  computing  the  ratio 
of  the  Ay  for  both  systems. 

For  the  system  without  clippers,  from  Eqs.  (It),  (5)  and  (13) 


11 

co 

1 

(16) 

(yAi>  *  H«Ni>  *  3> 

(17) 

The  normalized  signal  output  found  in  Eq.  (18)  is  denoted  by  primes. 


.  _  0  <VA>  ,  .  S 

<AyA>  =  7— r  =  (M  -  1) - — 

<yAL>  S+<H;> 

For  the  system  with  clippers,  the  normalized  signal  output  from 
Eqs.  (11),  (12)  and  (13 )  is  f  ound  to  be 


(18) 


(Ay)  =  l , l  =  (h  -  l)  -  <(a 

<yT> 


H  =  (n  _ 


\ 


\arcsin^r^^~N“^j/ 


(19) 


The  signal  response  ratio,  which  is  an  index  of  relative  qystem 
performance  with  respect  to  average  signal  output,  is  defined  to  be  the 
ratio  of  Eq.  (19)  to  Eq.  (l8). 


Equation  (20)  remains  to  be  evaluated  for  the  two  cases  indicated  in  Eqs. 

(1U)  and  (13). 

For  the  continuous  distribution  in  Eq.  (1I4.) ,  which  describes  a 
uniform  distribution  of  noise  power  on  a  logarithmic  basis,  we  have 
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Furthermore,  the  evaluation  of  the  expectation  of  the  arcsin  function 
in  Eq.  (20)  is  simplified  if  the  restriction  is  made  that  >  a^>  1. 

Since  the  detection  of  strong  signals  presents  no  problem,  the  restriction 

.  7 

is  almost  always  satisfied,  Thus, 


(23) 


Substituting  the  results  in  Eqs.  (22)  and  (23)  into  Eq.  (20)  yields  R^: 


(210 


If  a^  >  a^  »  1,  and  a^  =  b,  Eq.  ( 21+)  further  simplifies  to 
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(25) 


(b1^  -  l)2(b  -  1) 
b(£n  b)^ 


The  results  given  for  in  Eqs.  ( 2U)  end  (2$)  are  evaluated  for 
several  values  of  and  plotted  as  functions  of  b  in  Figure  3.  The 
results  show  that  for  sufficiently  large  b,  the  performance  of  the  system 
with  clippers  is  superior  to  that  of  the  system  with  the  clippers  removed, 
although  the  results  are  certainly  dependent  to  a  great  degree  on  the 
assumed  distribution  of  noise  power  given  in  Eq.  (lli).  Other  assumed 
distributions  show  much  less  improvement. 

For  example,  if 

f(N.)  = 


NL  <  Ni  <  NH 


elsewhere 


(26) 


the  relation  for  becomes 


(1  +  aH)^2  -  (1  +  a^2 

ro 

,  .  aH  +  aL 

aH  ’  aL 

- 1 

CM 

H 

_ 1 

(27) 


When  a^  is  large  compared  to  unity,  Eq. 


b  +  1 

(b1^2  +  l)2 


(27)  becomes 

-  for  b  =  1 
n 

< 

-  for  b-?-c= 
n 


(28) 


Very  little  advantage  of  the  system  with  clippers  is  shown,  even  for 
large  b. 
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NOPMAL/ZEp  AVE&A&E 
SIGNAL  RESPONSE  PATIO 
AS  A  FUfiJCTlPhJ  OF  isJPV 


VARl  A  77  OhJ 


Similarly,  if 


"A 


-  < 


^  nl<ni<nh 


elsewhere 


(29) 


7 

the  relation  giving  is' 


\-i 


(1  ♦  aL)V=  a  *  aH)VV  x  J  [a  t  aH)l/2_  l][(l  *  l] 


Im 


’(1  + "  aH)^"2+  Xj[(l  +  aL)^2"  x] 


Vl 

3  - 

i 

aH-  aL 

1*—^  +  £n  "H 


Vl 


!i 

&T 


(30) 


Equation  (30)  reduces  to  Eq.  (3l)  with  considerable  manipulation  if  a^ 
is  large  compared  to  unity. 


(b3/2  -  l)2 
<b  -  l)3 


In  b 


9n 


In  b 


for  b  3  0 
for  large  b 


(31) 


Evaluation  of  Eq.  (31 )  for  b  3  1000  yields  =  1.95.  The  system  with 
clippers  again  shows  very  little  advantage  over  the  system  without 
clippers. 

Next,  Eq.  (20)  is  evaluated  for  the  discrete  distribution  found  in 
Eq.  (15).  The  average  noise  power  becomes 


(Nj)3  p  Nl  +  (1  -  p)nh 


(32) 
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Also, 

<(1  +  ai)"l/j(l  +  aJ)"l/2>  -  p2(l  +  aL)_i  +  2p(l  -  p)(l  +  a^’^d  +  a 


+  (1  -  pT(l  +  a) 


-1 


[p(l  +  aH)^-  +  (1  -  p)(l  +  a^2] 
(l  +  a^)(l  +  a^) 


(33) 


Combination  of  Eqs.  (20),  (32)  and  (33)  yields  the  relation  for 
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p(l  +  aH)^2  +  (1  -  p)(l  +  a^^2]  [p(1  +  aL)  +  (1  -  p)(l  +  aR) 


(1  +  aL)(l  +  aH) 


p(c^2  -  1)  +  1 


2  r 


p(c-1  -  l)  +  1 


1  +  a. 


where 


c  = 


H 


(310 


(33) 


1  +  aT 


Note  that  c  ^b  for  .  aT  > 

ll 

The  maximum  value  of 


>1. 

h 


occurs  for 


P  • 


2  c1/ 2  +  1 
3(c1/2  +  1) 


1/2  for  c  =  1 

< 

2/3  for  c  -*  CD 


(36) 


The  maximum  value  of  R^  is  plotted  as  a  function  of  c  in  Figure  It. 

Also  of  interest  are  the  limiting  forms  of  Eq.  ( 3 U )  for  which  either 
the  percentage  of  noisy  transducers' is  small,  or  the  percentage  of  quiet 
transducers  is  small.  For  a  small  percentage  of  noisy  transducers,  (l  -  p) 
is  small,  and  Eq.  (3I4.)  becomes 
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_ /2^g£  ± _ 

Noqmal.i Ave£A&e 
S/&/VAI  /zesrpofijss  42 Ary o 
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r  =  Li  =  5  +  >k 

I  +  al  s>+  nl 


I 

I 
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A2-15 


v-1 


1  +  (c  -  3  +  2c  ^^)(l  -  p) 


(37) 


For  a  small  percentage  of  quiet  transducers,  p  is  small,  and  Eq.  (3I4) 
becomes 


*1  ~  n 


1  +  (2c^^  -  3  +  c~^)p 


(38) 


V.  Performance  Comparison  -  Ratio  of  Output  Signal-to-Noise  Ratios 

The  output  signal-to-noise  ratio  for  both  systems  is  defined  in  Eq. 

2 

(39)>  where  (Ay)  is  the  average  signal  output,  and  (a  >  is  the  variance 
of  the  output  of  the  low  pass  filter,  shown  in  Fig.  1,  averaged  over 


the  N . . 
1 


(39) 


The  performance  index,  R^,  considered  in  the  following  analysis  is  the 
square  root  of  the  ratio  of  r,  for  the  system  with  clippers,  to  r.  for 
the  system  without  clippers.  It  can  be  seen  in  Eq.  (ho)  that  R^  is 
closely  related  to  R, ,  which  is  calculated  in  detail  in  Section  IV. 


/  W2 

It) 

W  K°za2'\ 

It/ )  ; 


From  Eq.  (IjO),  it  can  be  seen  that  the  major  part  of  the  analysis  in 
determining  R^  involves  the  calculation  of  the  average  variance  of  the 
output  of  the  low-pass  filter  for  both  systems.  Both  distributions  of 
input  noise  power  given  in  Eqs,  (llj)  and  (15)  must  also  be  considered. 
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Since  all  the  input  signal  and  noise  time  functions  are  assumed  to 
be  gaussian,  the  sum,  x(t),  shown  in  Figure  1  is  also  gaussian  for  the 
system  without  clippers.  Under  these  conditions,^ 


R  .(t)  =  R  2(0)  +  2  R  2(t) 

yA  x  X 


Oil) 


where  R^(t)  is  the  autocorrelation  function  for  y(t),  and  Rv(x)  is  the 

2 

autocorrelation  function  for  x(t).  In  Eq.  (Ijl),  R  (0)  simply  represents 
the  average  output  of  the  squarer,  which  has  been  given  the  symbol  y^. 
Thus  the  variational  component  of  y(t)  has  the  autocorrelation  function 
2  Rx2(t). 

The  spectral  density  of  the  function  y(t)  is 


S  A(»)  =  -  2  R  (t)  e"v4~r7 

yk  n  x 


dT 


(1*2) 


The  variance  of  the  output,  z(t),  of  the  low-pass  filter  is  found  in 

Eq.  (I43),  where  H  (jco)  is  the  frequency  response  function  of  the  filter, 
z 


zA 


S  (co) 

y 


Hz(jo) 


dw 


(1*3) 


In  order  to  make  the  value  of  r^  large,  the  bandwidth  of  the  low-pass 
filter  must  be  much  smaller  than  the  bandwidth  of  S  (co),  the  input 

y 

spectral  density.  Under  this  assumption,  we  have 


a_„2='  S  (0) 

v 


zA 


Hw(j») 


dw 


(kb) 
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1 


Furthermore,  if  H  (jw)  has  the  form 
z 


Hz(jco) 


(US) 


1  +  0 


.  CO 


“f 

then  Eqs  (li2),  (UU )  and  ( U5 )  combine  to  yield 


°zk2  =  2  “t?  I  RJ(r>  dT 


F  x 
0 


(U6) 


For  the  "on  target"  situation,  the  function  R  (t)  is  simply  the  sum 

X 

of  all  the  matrix  elements  in  Eqs.  (l)  and  (2),  with  r.  .  =  0. 

^"3 


/  M 


R  (t)  = 
x 


^  Ni  +  H2S  p(x) 


(U7) 


and 


'  i=l 


M 


R 


x(x)  -[£  N^lfcs 

i=l  n  j-1 


P2(^) 


M  M  M  \ 

y  y  n.  n  .  +  2  y1  ^  +  i^s2!  p2(v) 

‘i*l  j=l  1  3  i=l  '  (U8) 


Combining  Eqs.  (U6 )  and  (U8)  and  averaging  over  the  IT  we  have 


(o,A2)-  2»?  Jm^2)  +  (M2  -  M)  (u }/2  +  2  K35  (li,)  +  mV 


u 


f 


p  (t)  dT 

0  (U9) 


The  ratio  of  the  square  of  Eq.  (16)  to  Eq.  (1*9)  gives  the  signal- 
to-noise  ratio  for  the  unclipped  system. 
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The  first  denominator  factor  is  nearly  a  perfect  square  if  M  is  large 
so  that 


(5D 


The  result  in  Eq.  (£l)  is  plotted  in  Fig.  $  on  logarithmic  coordinates. 
For  values  of  S/ (N_.)  significantly  less  than  l/li,  the  value  of  r^  varies 


Fig.  5  Normalized  Output  Signal- to-Noise  Ratio  as  a 

Function  of  Input  Signal-to-Average  Noise  Ratio 


as  the  square  of  S/(lh),  and  the  first  two  terms  (noise  x  noise)  can 
be  used  to  accurately  determine 

This  result  is  of  quantitative  value  in  the  analysis  of  the  system 
with  clippers  in  estimating  the  region  of  validity  of  the  result  for 
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\°z  In  this  case  only  the  noise  3:  noise  terms  are  easily  evaluated. 

The  assumption  of  low  input  signal-to-noise  ratio  is  not  unduly 
restrictive,  since  only  low  input  signal-to-noise  ratios  cause  any 
dete  ction • p  roblems . 

In  the  analysis  of  the  system  in  Fig.  1  with  clippers  present,  it 
is  assumed  that  the  signal  components,  s_.  (t),  of  the  inputs  to  the  clippers 
are  negligible.  This  assumption  makes  the  mathematical  analysis  for 
the  evaluation  of  tractable. 

The  output  of  the  square  law  device  is  given  by 


M  M 

y(t)  =  y  y  v.(t)  v.(t) 
i-l  Pi 


(52) 


The  autocorrelation  function  for  y(t)  is 


R  (t)  *»  E 

y 


M  H 


M  M 


\ 


L  Z  vi(t)  v3(t)  £  Z  Tk(t  +  t)  vn(t  *  T) 

i=l  j=l  k=l  m=l 


(53) 


The  evaluation  of  Eq.  (53)  is  helped  by  breaking  down  the  terms  into 
different  categories  as  shown  in  Table  1. 


Subscript  relations 

No.  of  terms 

Evaluation  and  comments 

i  =  j  and  k  =  m 

M2 

2 

Each  term  is  unity  since  v^  =  1 

i  =  0  and  k  /  m 

M3  -  M2 

These  terms  are  all  zero  because 
each  is  of  the  form 

E  fv  2v,  v  1  =  E[v,v  1=0 
l  i  k  rnJ  L  km J 

k  /  m 

i  /  j  and  k  =  m 

i 

It3  -  M2 

i  /  j  and  k  /  m 

-  2V?  +  M2 

Only  terns  of  the  type  where 
i  =  lc,  0  =  n;  and  i  =  m.  j  =  k 
are  non-zero.  2(H<1  -  M;  in 
number. 

Table  1  Classification  of  Terms  in  Eq.  (53) 
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In  view  of  the  information  in  Table  1,  Eq.  (53)  reduces  to 


R  (t)  =  M  + 

y 


m  n 

o:  i  - 

i=l  j-i+1 


v.(t)  v.(t  +  t)  E  v.(t)  v.(t  +  t) 
i  i  3  3 


C5U) 


From  Eqs.  (7)  and  (8),  Eq.  (5^)  becomes 


R  (t)  =  M  + 

y 


-  M)  -4?  J”  arcsin  p(x) 

rr  L 


(55) 


2 

In  Eq,  (55),  M  represents  the  squared  average  output  of  the  squarer, 
and  the  second  terra  is  the  one  which  represents  the  time-varying 
component  of  y(t).  By  a  development  similar  to  that  in  Eqs.  (Ill) 
through  (Ii6),  it  may  be  sham  that 


CTZ2  "  “F  I  Ry,(T)  ** 


(56) 


where  R  (t)  is  the  second  term  in  Eq.  (55).  Then 

y 


o  ^  Wp  (11^  -  H)  I  I  arcsin  p(r) 

V 


dr 


<  <<'N.  ) 
i' 


(57) 


Since  the  clippers  effectively  normalize  their  inputs,  a  is  not  a 

z 

function  of  the  N.  and 
1 


(58) 


The  substitution  of  Eqs.  (Ii9),  (57),  and  (58)  into  Eq.  ( I4.O )  yields 
the  expression  for  Rp. 
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It  should  be  remembered  that  only  the  first  two  terms  in  Eq.  ( U9 )  are 
used  in  Eq.  (59)  and  also  the  S  in  Eq.  (1;0)  is  neglected.  Furthermore, 
if  M  is  reasonably  large,  the  /  may  be  neglected  in  Eq.  ($9).  Then 

R2  ~  K  1  h  (60) 

and  K  is  the  square  root  of  the  ratio  of  integrals  in  Eq.  (59). 

The  integral  terms  in  Eq.  (59)  remain  to  be  evaluated.  Three 
representative  normalized  autocorrelation  functions  are  chosen  in  Eqs. 
(61),  (62)  and  (63). 

-ou  |t| 

P1(t)  =  e  (6l) 

P2(t)  =  e  (1  -  «2ItI )  (62) 


<1 


The  function  p^(t)  represents  a  process  with  a  low-pass  spectral  density 
having  a  half -power  frequency  at  p2(t),  a  process  with  center 
frequency  co2  and  bandwidth  2  co^;  a  process  with  center  frequency 

co_,  and  bandwidth  2^co^. 

Table  2  gives  the  exact  value  of  the  integral 
three  autocorrelation  functions. 
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Px(^) 

P2(x) 

P3(t) 

L 

1 

1 

I 

Hj 

2co 

1 

l*»2 

CD 

Table  2  Values  for  J  p^(t)  drr 

0 


It  should  be  noticed  that  the  integral  yields  the  following  general 
result  in  each  case: 


P^(t)  dr  = 


2  (process  bandwidth) 


(6U) 


The  evaluation  of  the  integral  involving  the  arcsin  function  is  most 
easily  accomplished  by  means  of  a  series  approximation. 


arcsin 


p(t)  =  P  +  i  P3  +  ^  P5  +  TTo  P7  +  ••• 


112 


(65) 


Squaring  - 


to  J 


[arcsin  p(r)]  ■  p2  +  .333  P^  +  .178  p6  +  .111;  p6  +  ...  (66) 

CD 

Evaluation  of  the  integral  J  [sin-'*'  p(t)]  dr  using  the  first  four 


.  -1 


terms  of  the  series  approximation  of  [sin  J'  p(t)J  yields  the  values 
listed  in  Table  3.  The  integration  for  p^(^)  is  accomplished  by 
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replacing  the  even  powers  of  the  cosine  function  in  Eq.  (63)  by  the 
average  value  of  the  function.'4 


This  is  an  additional  approximation,  tut  the  result  should  be  accurate 
for  low  (p  .  It  can  be  seen  that  the  values  for  the  arcsin  integral  are 
remarkably  close,  regardless  of  the  nature  of  the  autocorrelation 
function. 

From  Tables  2  and  3 ,  the  factor  IC  is  evaluated: 

K  =  [  \  ^  =  .90  (68) 

\  *62  / 

Combining  Eq.  (68)  with  Eq.  (60),  we  have 

M  >  >1  (69) 

<N1)>>S 

An  obvious  consequence  of  Eq.  (69)  is  that  the  results  for  the  normalized 
signal  response  ratio,  IL ,  in  Figs.  3  and  4  nay  be  applied  directly  to 
the  determination  of  the  ratio  of  signal- to-noise  factors, 
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VI.  Conclusions 


The  results  in  Section  IV  for  the  average  signal  response  ratio, 

R^,  have  shown  that  the  system  -with  clippers  is  superior  to  the  system 
without  clippers  if  the  variation  in  transducer  noisiness  is  fairly 
large.  However,  in  the  case  of  the  continuous  distribution  for  the 
noise  powers,  the  value  of  R^  is  dependent  to  a  great  degree  on  the 
form  of  the  probability  density  function  that  the  1L  exhibit.  Very 
little  advantage  is  shown  for  most  monotonic  density  functions,  except 
the  one  which  varies  inversely  with  the  noise  power.  This  particular 
density  function  describes  the  situation  in  which  the  number  of  transducers 
having  noise  powers  above  the  geometric  mean  of  the  extreme  values  is 
equal  to  the  number  having  noise  powers  below  the  geometric  mean. 

Qualitatively,  this  result  agrees  with  the  result  obtained  for  the 
two-valued  discrete  distribution  of  noise  powers.  Maximum  benefit 
offered  by  the  system  with  clippers  is  obtained  when  approximately  half 
the  transducers  produce  the  lower  noise  power,  and  half  produce  the 
upper  noise  power. 

The  significant  result  of  Section  V  is  the  effective  equivalence  of 
the  two  performance  indices.  From  Eqs.  (60)  and  (69),  R2  is  just  a 
multiple  of  R^.  The  numerical  factors  entering  into  the  multiplier  are 
consequences  of  the  t;rpes  of  elements  used  in  processing  the  signals. 

The  factor  n/ 2  in  Eq.  (60)  clearly  arises  from  the  presence  of  the 
infinite  clippers,  and  the  factor  K  arises  from  both  the  clipping  and 
squaring  operations. 
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Finally,  the  result  for  Rg  is  seen  not  to  depend  greatly  on  the 
autocorrelation  function  for  the  signal  and  noise,  for  the  three 


different  functions 
exact  properties  of 
since  the  filtering 


chosen,  the  results  are  essentially  the  sane, 
the  low -pass  filter  are  not  important,  either, 
effect  is  exactly  the  same  for  both  systems. 


The 
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SUMMARY 


The  performance  of  an  optimal  detector  is  compared  with  that  of  a 
standard  detector;  each  uses  the  same  array  of  hydrophenes  in  an  effort 
to  detect  the  presence  of  a  directional  Gaussian  signal  in  an  isotropic 
Gaussian  noise  field.  Output  signal  to  noise  ratio  is  used  as  a  figure 
of  merit.  The  following  results  are  obtained: 

1)  If  signal  and  noise  both  have  flat  power  spectra  of  sufficient 
bandwidth  so  that  the  waveshapes  received  by  different 
hydrophones  are  uncorrelated,  then  the  optimum  detector  and 
standard  detector  are  identical  in  performance. 

2)  If  the  noise  spectrum  is  flat  as  in  l)  tut  the  signal  spectrum 
is  not,  the  optimal  detector  is  better  than  the  standard 
detector  by  an  amount  depending  on  the  differences  of  the 
spectral  properties  of  signal  and  noise. 

3)  If  the  signal  and  noise  spectra  are  identical  in  shape  the 
optimum  detector  is  theoretically  capable  of  detecting  a  signal 
with  any  desired  degree  of  certainty  in  an  arbitrarily  small 
length  of  time.  In  practice  the  noise  always  contains  more  high 
frequency  components  than  the  signal.  This  phenomenon  limits 
the  improvement  attainable  through  optimum  detection  to  a 
moderate  amount  in  most  practical  cases. 

I4)  Knowledge  of  small  differences  in  the  shapes  cf  signal  and 
noise  spectra  is  of  limited  value  in  detection. 

5)  Correlation  between  noise  disturbances  received  by  different 
hydrophones  appears  to  degrade  the  performance  of  the  standard 
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detector.  Cne  case  of  this  type  has  been  worked  cirt  in  detail. 

It  assumes  a  linear  array  of  five  hydrophones  spaced  at  two-foot 
intervals  and  signal  and  noise  spectra  falling  off  with  the 
second  power  of  frequency  above  1^00  cps. 

6)  An  attempt  is  made  to  estimate  the  probable  high  frequency 

behavior  of  the  signal  spectrum  in  order  to  set  realistic  bounds 
on  detectability  by  means  of  the  optimum  system.  Using  a  model 
for  signal  transmission  due  to  E chart,  a  range -dependent 
expression  is  obtained  for  the  ratio  of  performance  figures  of 
optimal  and  standard  detectors.  At  ranges  of  practical  interest 
the  ratio  assumes  values  not  larger  than  0. 

In  general,  if  signal  and  noise  are  characterized  only  by  total  power 
and  spectral  properties,  one  is  led  to  the  conclusion  that  the  optimal 
detector  exhibits  only  a  limited  advantage  over  the  standard  detector 
unless  the  spectral  properties  of  signal  and  noise  differ  drastically. 
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I.  Introduction 


This  report  deals  with  the  problem  of  optimal  detection  of  a  weak 
directional  signal  in  a  very  much  stronger  isotropic  noise  field  by  means 
of  a  given  array  of  hydrophones.  Both  signal  and  noise  are  assumed  to  be 
Gaussian  random  processes  with  known  power  spectra.  The  investigation 
has  two  primary  purposes: 

1)  To  set  bounds  on  the  detectability  of  signals  characterized  only 
by  their  total  power  and  spectral  properties  in  a  noise  background 
similarly  described.  If  a  signal  is  not  detectable  with  the  instrumentation 
considered  here,  then  it  is  not  detectable  at  all  with  the  given  array 
unless  additional  properties  of  signal  or  noise  are  known. 

2)  To  determine  how  far  conventional  data  processing  procedures  fall 
below  the  optimum.  As  in  many  optimization  problems,  the  best  instrumentation 
may  well  be  too  complex  for  practical  realization  and  serves  primarily  as  a 
standard  of  comparison  for  more  realistic  procedures. 

It  is  well  known  that  the  optimization  problem  described  above  is 
formally  solved  by  a  likelihood  ratio  detector.  The  physical  realization 
of  such  a  detector  has  been  discussed  by  F.  Bryn.  Bryn  goes  on  to 
compare  the  performance  of  the  optimal  detector  with  that  of  standard 
detectors,  primarily  on  the  basis  of  a  figure  of  merit  which  he  designates 
as  the  "array  gain,"  G^  =  G(con),  It  may  be  defined  as  follows 


G(w  ) 
n 


/Contribution  to  average  detector  output' 
of  signal  at  frequency  ccy, _ 

\Contributi on  to  average  detector  output 
\of  noise  at  frequency  «n 


/ Input  signal  to 
noise  ratio  at  each 1 
,  hydrophone  at 
\ frequency 


Finn  Bryn,  Optimal  Signal  Processing  of  Three  Dimensional  Arrays 
Operating  on  Gaussian  Signals  and  Noise.  J.  of  the  Acoustical  Society 
of  America,  Vol.  3h  No.  3,  March  1962,  pp.Tiw-297. 
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The  array  gain  is  therefore  a  steads'-  state  (bC)  performance  criterion 
depending  only  on  array  geometry.  For  the  purpose  of  the  present  discussion 
a  more  useful  performance  criterion  is  the  signal  to  noise  ratio  at  the 
output  of  the  detector.  Its  computation  will  therefore  be  the  first  task. 
II,  Nomenclature 

The  first  part  of  this  report  will  draw  heavily  on  Bryn»s  work  and 
his  nomenclature  will  therefore  be  used  as  far  as  practicable.  Thus  the 
signal  e^(t)  received  by  the  i^1  hydrophone  will  be  represented  over  the 
observation  interval  0  <  t  <  T  by  the  Fourier  series 


ei(t) 


\ 

/ 


A.(n) 


CCS 


to  t 

n 


Bi(n) 


sin  wt 1 
n  J 


n 


where  co^  =  and  the  summation  extends  over  all  frequencies  in  the  signal 
band.  Signal  and  noise  are  assumed  to  be  independent  Gaussian  processes 
with  power  spectra  S(w)  and  N(to)  respectively.  Then  if  T  is  large  compared 
to  the  correlation  time  of  the  signal  aid  noise  processes  cne  can  make  the 
following  statements: 

a)  A.(n)  and  B.  (n)  are  Gaussian  random  variables  with  zero  mean  and 
1  l 

variance  [S(wn)  +  N(wn)j  AoCl  j  where  Ato  =  ~ 


B)  A.(n)  is  statistically  independent  of  A.(m)  for  all  j  (including 

J 

j  *  i )  if  n  /  m 

c)  A  (n)  is  statistically  independent  of  B_^(n) 

However,  A_^(n)  is  not  in  general  independent  of  A^.(n).  For  ease  in  later 
manipulation  Bryn  introduces  the  symbolism 


x±(n)  =  < 


A.(n) 


B.(n) 

l 


1  <  i  <  M 
M+l  <  i  <  2M 


(1) 
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M  is  the  number  of  hydrophones  in  the  array.  The  symbol  (  )g  is  used 
to  indicate  the  statistical  average  of  the  bracketed  quantity  when  only- 
signal  is  present.  {  ^  designates  a  similar  average  when  only  noise  is 
present.  Still  following  Bryn  we  define 


<Xi(n)  "  S(a)n)  A<o  Pih(n) 

^  x^(n)  Xb(n))N  =  N(o>n)  At)  Wn) 


Thus  p^(n)  and  q^(n)  are  the  elements  of  the  normalized  correlation 
matrices  of  signal  and  noise  respectively: 


P(n)  = 


Q(n)  = 


P^n)  .  .  . 

pl,M(n) 

pl,M+l(n) 

PM,l(n)  ■  • 

•  pII,M(n) 

pM,M+l(n) 

PM+1,M+1^ 

• 

P2M,l(n)  •  ' 

•  p2M,M(n) 

P2M,M+1^ 

qn(n)  .  . 

ql,M+l^ 

qH^l(n)  .  . 

qM,M+l^ 

qM+l,l^ •  ■ 

•  %n,M(n) 

qM+l,M+l^ 

q2M,l(n)  •  • 

•  q2M,M^n) 

q2M,M+l^ 

PM+1,2M^ 


J2M,2MV 


ql,2M^ 
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(5) 


Finally  it  will  be  convenient  to  have1  symbols  for  the  square  matrices 
consisting  of  the  first  M  rows  and  columns  of  Eqs.  (It)  and  (5); 


and 


P^n)  = 


Pll(n)  •  •  •  plM(n) 


Pmi^h)  .  .  . 


Q^n)  = 


qn(n) 


qiM(n) 


(6) 


(7) 


The  inverse  of  Q(n),  and  hence  of  Q^(n),  will  be  assumed  to  exist.  Since 
Q(n)  is  a  correlation  matrix  this  represents  no  significant  restriction. 

III.  Optimum  Detection 

Bryn  gives  the  likelihood  ratio  LR,  the  output  of  the  optimum  detector, 


as 


LR  =  1  +  Z  r(c0n}  “  ^r?\ 

n  l 


(8) 


where 


w(»n) 


2M  2M 

•hi  I 


S(o3  )Aw 


2M  2M 


n 


Q.  .(n)  Q..(n) 


fe  ^  lQ(n)l  >  “ 

(?) 


is  the  ji  element  of  the  matrix  Q 


-1 


The  only  term  in  Eq.  (3)  dependent  on  the  signal  and  noise  actually  received 

by  the  various  hydrophones  is  W(co  ).  Hence  the  change  in  DC  output  due 

n 

to  the  presence  of  a  signal  is  simply 


A  (DC  output)  =  77  / 


n 

i 

'S(“n)_ 

"  2  L 

n 

“<“nl 

2  2M  2M  2M  2H 
^  ^ 


L  £ 


P.Jn) 

i=l  h=l  j=l  k=l  '^n'* 


/ 


V">  aki-.(n) 


lQ(n)l  jk 


p.,  (n)  (10) 


Using  the  fact  that  the  matrices  P  end  Q  are  symmetrical  this  expression 
can  be  written  in  the  matrix  form 

12 


1  V 

S(“n) 

2  L 

N (co  ) 

n 

L  n  . 

i  r 

rs(uo  )i 
n 

2  2n 

n 

N(co  ) 

L  n  J 

Tr  P(n)  Q-1(n)  P(n)  Q“x(n) 


-1/ 


Tr  <|  [P(n)  Q_1(n)l  \ 


(11) 


where  Tr)  jis  the  trace  of  the  indicated  matrix. 
Bryn  has  shown  that 


|Q(n)|  lQ(n)|  3k 


2M  2M  2M  2M  „  ,  x  .  ,  s 

7  v  y  y  p.  <„)  ^i(n)  ;*»(n) 

Li  i—>  Li  ‘L-i  ih  Iro  /■-. 

i=l  h=l  j=l  k=l 


or  in  matrix  form 


p.,  (n)  =  2 


M  M  A  (  \ 

kk^ni  ^ 


l2 


[p(n)  Q"1(n) 

x. 

2 

1 

>  =  2  Tr  . 

1 

1  L 

-1/ 


(12) 


(13) 
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Hence  Eq.  (ll)  becomes 


A (DC  output) 


Tr 


P1(n)  Q1""1(n) 


(14) 


Next  it  is  necessary  to  compute  the  rms  fluctuation  of  the  detector  output. 
This  quantity  will  be  designated  by  the  symbol  D( output).  If  the  input 
signal  to  noise  ratio  is  very  lew,  a  condition  already  assumed  in  the 
derivation  of  Eq.  (8),  essentially  all  of  the  output  fluctuation  can  be 
attributed  to  noise  and  the  signal  component  can  be  ignored.  Then 


D^(  output)  =< 


L 


W(w  ) 
n 


A 


(15) 


From  Eq.  (9) 


L 


S(w  )  Aw  5(w  )  Aw 

< _ n _ m _ 

N‘"(co  )(Aw)^  N^(w  )  Aw 
n  ra 


x 


2M 

L, 


Q.  .  (n)  Q.  1 .  i  (m)  Q. ,  (n)  Q.  (m) 
p^n)  p.,h,(m)  ldl  kh 


|Q(n)l  |Q(m)|  |Q(n)|  |Q(m)| 


(x  (n)  x,  (h)  x  (m)  x,t(m)) 

'  3  K  J  K  /  n 


.1.1,1.  1  .. 

i,j,h,k  =  1 


(16) 


Now  since  |xj(n),  x^.(n),  x^.(m),  x^m)  is  Gaussian  in  four  dimensions 
^Xj(n)  x^(n)  x^m)  x^m)^ 

(^(n)  Xk(n)};j(x.<m)  x^m))  ♦  (x.(n)  x^m))  (x^n)  x^mj)  *  (: x.(n )  x^m)^ 


V- 


"V 


V. 


term  (l) 


term  (2) 


V- 


'V 


term  (3 ) 


(17)' 


•^J.  H.  Laning  and  R.  H.  Battin,  Random  Processes  in  Automatic  Control, 
pp.  82-84.  McGraw-Hill  Book  Co. , 


The  substitution  of  Eq.  ( 17 )  into  Eq.  (16)  lea.ds  to  three  terns  that  will 
bo  treated  separately. 

\x.j(n)  x^n))  x^(m)^)  -  q^k(n)  qj^m^  N(“n)  A“  N(%j)  Aw  (18) 

Hence  term  (l)  is 


1  ,  S(co  )  Aco  S(co  )  A<o 

\  v1  v  n  m 


c  L  X 


N(co  )  Aco  N(co  )  Aco 
n  m  n  m 


2M 

V 


PA*(rn)  .JJ. 


Q^(n)  Q^m)  Qj^Cn)  QlJut(ra) 


W 


q.Jn)  q.i,j(m) 


/  ,  "ih .  I«(n)|  IQWI  |Q(n)l  IQMI 


L—  .  4 

i,j,h,k  =  1 


2M 


= 


1  V  S(“n5  A“ 

2  / 


Q.,  (n)  Q,,  (n) 


cih 


kh 


N2(co_  )  ( Aco)"  / |  ^  |Q  (n  )t  I Q  ( n  )| 


q,Jn) 


n 


i,j,lc,h  =  1 


(19) 


The  last  step  of  this  computation  follows  from  inspection  cf  Eq.  (9). 
Returning  to  Eq.  (17)  and  using  conditions  b)  and  c)  on  page  2 


qjji(n)  q^n)  N2(«n)(A<x>)2 

=  < 

0 


m  =  n 


m  f  n 
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(20) 


Thus  term  (2)  assumes  the  form 


2M 


1  V 
U  L. i 

n 


S(co  ) 
n 


-1 2 


N(co  ) 
n 


p.h(n)  p.Un)  ^  ^  ^  ,.Xn)  ,(n) 

lh  A  |Q(n)l  |Q(n)l  |Q(n)l  lQ(n)l  «  lck 


±j jjhjk'=  1 


’ll 


n 


rs(»„)i 

2 

Tr 

P(n)  Q-1(n)  Q(n)  Q_1(n)  P(n)  Q-1(n)  Q(n)  Q"1(n) 

N(co  ) 

L  n  J 

- 

n 


S(co  ) 
n 


Tr]  |  P(n)  Q“L(n)j2 


1  V  feki  Tr 

2  ^K> 


P^n)  Q1”1(n) 


The  last  step  results  from  Eq.  (13). 
Finally,  returning  once  more  to  Eq.  (17) 


x.(n)  x/m))  (x,  (n)  ximj>  = 
J  K  m  N  k  J  N 


qiki(n)  qkj»(n)  N2(u3n)(Aco)2 


0 


(21) 


m  =  n 


m  f  n 


(22) 
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Thus  term  (3)  reduces  to 


2M 


1  r 

H  L, 

n 


S(w  ) 

n 


N(a>  ) 
-  n  ■ 


\  '  Q±1(n)  QjjCn)  Q..(n)  Q.u(n) 

)  p  (n)  p.  (n)  -ii - a - — — -  —  q  .ji(n)  q.  in) 

/  ,  ‘  lQ(n)l  ]Q  (n  )|  |Q(n)l  lQ(n)l  J  J 


S(co  ) 
n 


u  U 

|n(w  ) 

n 

L  n 

iy\ 

,sun) 

2  L  1 

|  N  (co  ) 

n 

l  n 

Tr 


P(n)  Q  ^(n)  Q(n)  Q  ^(n)  P(n)  Q~^(n)  Q(n)  Q-1(n) 


.  2 


Tr 


P1(n)  Q1"1(n) 


(23) 


Substituting  Eqs.  (19),  (2l)  and  (23)  into  Eq.  (l5)  we  obtain 

2 


D  (output)  = 


/ 


n 


f  5(oo  ) 

r  in 

n 

Tr  P, (n)  Q  ~1(n) 

N(w  ) 
v  n 

[  1  1  J 

(21*) 


Equations  (ll*)  and  (2l*)  together  determine  the  ms  signal  to  noise  ratio  at 
the  output  of  the  detector 


A (DC  output) 
D( output ) 


'7^  Tr  D  (n)  «  -kn)] 

M»(®J  L  1  1  Jj 


(25) 


n  n 


Equation  (25)  may  be  rewritten  conveniently  in  terms  of  Bryn’s  "array  gain" 

(26) 


_M  _H 

G(“n)  =  L  )_,  Pih(n) 
i=l  j=l 


Qih(n) 


|  Q  (n  )| 


rhis  expression  is  similar  to,  but  not  identical  with,  an  expression 
described  by  Bryn  as  the  output  signal  to  noise  ratio  (SQ).  Since  Bryn 
gives  no  derivation  or  physical  definition  for  SQ  it  is  not  clear  whether 
he  is  in  fact  referring  to  a  closely  related  measure  of  performance. 
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or  in  matrix  form 


G(w  )  =  Tr  (n)  Q  ~J  (n)  (27) 

n  l  -  1 

Hence  Eq.  (25)  becomes 

A (DC  output) 

D( out out) 

The  frequency  index  n  ranges  over  the  entire  band  to  which  the  hydrophones 
are  sensitive. 

If  T  is  sufficiently  large  so  that  the  signal  and  noise  spectra  are 
essentially  constant  over  a  frequency  interval  of  Au>  =  rad/ sec. ,  Eq. 

(28)  can  be  expressed  in  integral  form 

m.  .y.tput)  ,  jT  hr-  -±nl  0(7)  to 
D(  output)  ^  1/  "n"  R 

(29) 

and  are  respectively  the  lower  and  upper  bounds  of  the  frequency  range 
being  processed  by  the  hydrophones. 

Two  special  cases  are  of  particular  interest. 

a)  If  signal  and  noise  spectra  are  identical  in  shape ,  a  situation 
closely  resembling  conditions  encountered  in  practice  if  one  ignores 
periodic  components  of  signal  and  noise 
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co  CO 

where  O  =  j  S(co)  dto  and  N  *  j  n(co)  dw  represent  the  total  signal 
0  0 

power  and  noise  power  respectively.  At  sufficiently  high  frequencies  the 
correlation  between  noise  disturbances  received  by  different  hydrophones 
tends  to  zero,  that  is 

( 

1  for  i  =  h 


qih(n)  -> 


0  for  i  /  h 


(31) 


Then  Q^(n)  reduces  to  the  identity  matrix  and  G(ccn)->  TrTF^(n)  .  Thus 
G(co)  approaches  a  positive  constant  so  that  Eq.  (30)  tends  to  infinity  if 
one  permits  to  grow  without  bound  In  other  words,  perfect  detection  is 
theoretically  possible  in  a  finite  time  under  the  stated  conditions.  This 
is  not  unreasonable,  for  the  ratio  of  signal  power  to  noise  power  is  the 
same  in  any  frequency  band  of  a  given  width.  Thus  the  optimum  detector 
can  be  thought  of  as  operating  simultaneously  on  an  infinite  number  of 
bands,  extracting  equally  significant  information  concerning  the  presence 
of  a  signal  from  each.  In  practice,  of  course,  the  noise  has  components 
that  do  not  fall  off  with  frequency  as  rapidly  as  the  signal  so  that  the 


input  signal,  to  noise  ratio 


S(co) 

NlcoJ 


tends  to  zero  for  sufficiently  large  u>. 


Noise  electrically  generated  in  the  hydrophones  or  associated  circuitry 
will  certainly  have  this  property  even  if  the  acoustical  noise  does  not. 
Determination  of  the  frequency  at  which  the  effective  input  signal  to 
noise  ratio  begins  to  decline  sharply  is  therefore  an  important  practical 
problem. 

b)  If  the  hydrophones  are  sufficiently  far  apart  so  that  Q^(n)  0 

for  all  i  /  h  throughout  the  frequency  range  processed  by  the  hydrophones, 


the  array  gain  is  simply  Tr 


P1(n) 


If  one  further  assumes  that  the  array 
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is  linear  and  mechanically  steered  on  target  or  that,  for  a  general  array, 
appropriate  delays  are  introduced  at  each  hydrophone  to  compensate  for 
relative  differences  in  distance  to  the  target,  then  the  elements  of  the 
matrix  are  all  unity.  It  follows  that  the  array  gain  is  simply  equal 
to  M,  the  number  of  hydrophones  in  the  array.  In  that  case  the  output 
signal  to  noise  ratio  is 


an 


A  (DC  output)  /jF  M  f  *lS(u) 

D(output)  v  2n  \j  J  N(u) 


d£0 


(32) 


The  general  expression  for  output  signal  to  noise  ratio  [Eq.  (29) J 

also  deserves  some  further  comment.  The  linear  dependence  of  detectability 

on  Vt  was,  of  course,  to  be  expected  on  general  statistical  grounds.  It 

is  interesting,  however,  that  knowledge  of  small  differences  in  the  shapes 

of  signal  and  noise  spectra  contributes  almost  nothing  to  the  detection 

process.  Differences  in  spectral  shape  become  important  only  if  they 

result  in  a  value  of  ^7 — -  very  much  larger  than  —  over  some  frequency 

N  (co)  N 

range.  A  case  of  this  type  would  be  cnc  in  which  a  finite  amount  of  signal 
power  is  concentrated  in  a  very  narrow  band  (i.e.,  the  signal  contains  a 
qua si-sinusoidal  component).  Because  the  integrand  is  squared  Eq.  (29) 
then  leads  to  a  large  value  of  output  signal  to  noise  ratio  even  if  the 
total  power  of  the  quasi- sinusoidal  component  is  quite  small.  One  must 
keep  in  mind,  however,  that  the  optimum  detector  now  depends  heavily  on 
knowledge  of  the  precise  frequency  range  in  which  the  quasi-sinusoidal 
signal  component  is  located.  In  most  practical  situations  it  would 
probably  be  unrealistic  tc  assume  such  knowledge.  An  interesting  problem 
that  has  not  been  solved  to  date  would  be  the  detennination  of  signal 


A3-1A 


detectability  in  the  sense  ol'  Eq.  ( 29 )  when  the  signal  contains  one  or 
more  quasisinusoidal  components  whose  frequencies  are  unknown,  or  at  least 
not  known  accurately. 

IV.  Standard  Detection 

The  performance  of  the  optimum  detector  wall  now  be  compared  with  that 
of  a  standard  detector  using  the  same  array  of  hydrophones.  The  standard 
detector  simply  adds  the  outputs  of  the  various  hydrophones,  squares  the 
sum  and  averages  the  resultant  signal  over  a  time  interval  T  (see  Fig.  l). 


Practical  instrumentations  frequently  include  'clippers  and  digital  time 
delays.  These  are  omitted  here  for  two  reasons: 

1)  It  is  desired  to  compare  the  likelihood  ratio  detector  with  the 
best  version  of  the  standard  detector.  Under  the  conditions  assumed  in 
this  computation  clipping  and  digital  delay  are  known  to  degrade 
performance. 

2)  The  degradation  in  performance  due  to  clipping  and  digital  delays 
has  been  discussed  elsewhere^  and  it  is  therefore  not  necessary  to 
introduce  this  added  computational  complexity  here. 

It  will  be  assumed  that  the  array  has  been  steered  "on  target,"  so 
that  the  signal  component  s(t)  of  each  hydrophone  output  is  the  same.  The 

^P.  Rudnick,  Small  Signal  Detection  in  the  Dimus  Array.  J ,  of  the 
Acoustical  Society  of  America,  Vol.  32  No.  7>  July  I960,  pp.  871-677. 
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th 


noise  output  of  the  i  hydrophone  will  be  designated  as  ru(t).  It  is 

assumed  to  have  aero  mean  and  to  be  independent  of  s(t).  Then  from  Fig.  1 

the  DC  component  of  y_  is 
/'  3 


'  M 

2 

,  v  M  M 

y3j  ‘  & 

y  |s(t)  +  n^(t)| 

’  ■  1,2  El  [»<*>] 2  [ + 1  Le 

ni(t) 

i-i '  1  i 

l  -  i=l  j=l 

Thus  the  increase  in  DC  output  due  to  the  signal  is 
A(DC  output)  =  M2  e|  s(t)j  I  =  ifi S 

l  ' 


(33) 


(3U) 


where  £  is  again  the  average  signal  power  at  each  hydrophone. 

Calculation  of  the  rms  output  fluctuation  is  perhaps  accomplished  most 
readily  by  integrating  the  continuous  portion  of  the  power  spectrum  W^(co) 
of  y^.  For  direct  comparison  with  the  optimal  detector  the  smoothing 
filter  will  be  defined  simply  as  a  device  which  averages  y  (t)  over  T 


seconds.  Its  weighting  function  is  therefore 

h(t)  = 


1 

T 


0  5  t  5  T 

elsewhere 


(35) 


By  Fourier  transforming  one  obtains  the  transfer  function 

T 

/- 

H(co)  =  l  j  dt  =  i-“  cos  +  1 

J 

0 


Hence 


|H(co)  |  2  =  1  -  2  cos  “T 

(coT)2 


/  .  oT 
sin 

ooT 

2 


(36) 


(37) 
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Therefore  the  variance  of  y^  is 

co  2 

r  <  .  » T  \ 

o  /  sin  T 

D  (y3)  “  /  w2^  ' '  coT  ' '  d0  (3^ 

J  T 

0 

where  is  the  power  spectrum  cf  y2> 

For  low  signal  to  noise  ratios  the  smoothing  time  must  be  long 
compared  to  the  inverse  bandwidth  of  w^(  w).  Under  these  assumptions  W^(  «) 
in  Eq.  (38)  can  be  considered  constant  at  W  (0)  =  lim  W9(w).  Thus 

r°[sin  °°t  \  2 

B2(y3)  -  W2(o)  /  -gj-M  d»  =  !»2(0)  (31 

jo'  T  ' 

WglO)  is  calculated  most  readily  from  the  corresponding  autocorrelation 
function  t  (t) 

CO 

W2(w)  =  ^  |  I2(t)  dr  (k( 


Now  ^(t)  =  Efy2(t)  y2(t  +  t) 


M  M  M  M  r 

r  r  r  J 


=  ^  E<1  [s^  +  +  +  T)  +  n^(t  +  t)] 

i=l  j=l  h=l  k=l  l 

\ 

[s(t  +  t)  +  nk(t  +  t)]  j 

M  M  M  M  r 

■  ^  )  y*  )_j  £|  s2(t  +  t) 

i=l  .1=1  h=l  k=l  1 


s(t)  s(t  +  T)[n.(t)  n^(t  +  t)  +  n^t)  nh(t  +  t)]  +  s2(t)  n^t  +  t)  nk(t  +  t) 

2  I 

+  s  (t  +  t)  n^t)  n  (t)  +  n^t)  n  (t)  n^( t  +  t)  nk(t  +  t)1 
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The  symbol  E|  j  is  used  to  designate  the  mathematical  expectation 
of  the  bracketed  quantity. 

The  first  term  of  Eq.  (i+l)  results  from  the  intermodulation  of  sigaal 
with  signal,  the  next  three  from  intermodulation  of  signal  with  noise,  and 
the  last  term  frem  internodulation  of  noise  with  noise.  When  the  signal 
to  noise  ratio  is  low  (the  situation  of  primary  interest)  the  last  tern  • 
dominates.  One  can  therefore  approximate 

M  H  M  M 

l  VI  . 

i=l  j=l  h=l  k=l 


n.  (t)  n.(t)  n,  (t  +  t) 
i  3  h 


\(t 


t) 


An  equivalent  approximation  was  made  by  Bryn  in  the  derivation  of  the 
likelihood  ratio  detector.  When  the  noise  is  Gaussian,  one  con  use  Eq.  (17) 
to  express  Eq,  (i|2)  in  the  form 


|  n.CtJnUtJn^Ct  +  x)nk(t  +  t) 

r 


<.  n.  .  . 

I  1  °  ) 

-  E|ni(t)n..(t)|>  e| n^( t  +  x)nk(t  +  t)|  +  E^Ctjn^t  +  -r )^E^ n^( t)nk( t  +  t)|> 

j  \ 


+  E<  n.  ( t)n,  (t 
|  1  k 


+  t)  l  E  n  .  (t)n,  (t  +  t) 


)  . 


(to) 


Introducing  the  nomenclature 


E<jn.  (t)  n.(t  +  t)1  =Nq.,(x) 


id 


(UU) 


Eq.  ( L;3  )  becomes 


Es' ni(-b)  n^(t)  n^t  +  r)  nk(t  +  t)1  =  f\[  jq^.CO)  qhk(0)  +  ciih(T) 

qjk(T) 


j 


+  %k(T)  Vt)]  (h5) 
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N  is  the  average  noise  power  at  each  hydrophone  (as  in  part  III)  and 


q^(T)  is  the  normalized  crosscorrelation  between  the  noise  components 
at  hydrophones  i  and  j.  Using  Eqs.  (39),  (UO),  (1|2)  and  (U3)  and 
recognizing  that  the  term  q.  .(0)  ^^(0)  contributes  only  EC  power,  one 
obtains 


M  M  M  M 


D 


2(y3)  =  --  '  ^  )  )  I 


T  c_,  L,  L,  I  'W'^jk' 

i=l  j-l  h=l  k=l  ^ 


U6) 


Hence  the  figure  of  merit  corresponding  to  Eq.  (29)  is 


A (DC  output) 

D(y3) 


M 


2  5 


r  ,,-2 


M  M  H H  « 


V  L  Ilf  qit.M  q)kM  dT 


i=l  j=l  h=l  k-1  _c= 


V2 


j£  -\/l  _ Ml 

j<J  r  2  r  M  M  M  I-i  « 


/  . 


V  v  \ 
/ ,  /  ,  / 


%h(T)  a-jk(T)  dT 


v* 


(hi) 


i=l  j=l  h=l  k=l  _co 

CO 
s 

Thus  the  figure  of  merit  depends  rather  critically  on  j  q^(T)  q^^)  dr. 


Two  cases  will  be  considered. 
l)  Uncorrelated  Noise 

If  the  hydrophones  are  separated  sufficiently  so  that  q  (t)  =  0  for 


i  /  h,  then 


CD 


J 

MCO 


qih(x)  qjk(T)  dx  =  > 


f  P2(-t) 


dr 


i  =  h,  j  -  k 


elsewhere 


(U8) 


^This* corresponds  to  the  assumptions  made  in  Frogress  Report  No.  2. 
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p(x)  =  q^^(x)  is  the  noiualizod  t.ut,csorrol.*\.,i  on  £-,\i\yU<-n  <r  +.K<-.  >i-,j  r,r,  a+, 
each  hydrophone.  When  Eq.  ([$)  is  satisfied 


A (DC  cutput) 

c(y3) 


-fi  V!  f- 


M 


w72 


(U9) 


p  tx)  dr 


-CO 


The  equivalent  result  for  the  likelihood  ratio  detector  is  Eq.  (32). 

A  typical  situation  in  which  Eq.  (U8)  is  at  least  approximately 

satisfied  would  be  the  following:  The  noise  has  a  power  spectrum  flat 

over  a  wide  band  0  <  o  <  o  and  vanishing  for  co  >  o  .  Then 

a  a 


sin  co  x 

p(t)  =  - ~  ,  w  large 

ox  a 


($0) 


Thus 


p  (x)  dx 


.  2 

sin  co  x 


—CO 


(<V>‘ 


,  n 
drr  -  — 
to 


(5i) 


Hence  Eq.  (U9)  assumes  the  form 


A(DC  output)  _  5 


D(y3) 
N 


M 


To 


2n 


M 


(52) 


Introducing  the  symbol  Nq  =  ~  ,  the  value  of  the  noise  spectral  density 


for  0  <  o  <  o 


A (DC  output) 

u(y3) 


S  -.IT 


N  V  2rt  "Vco 
o  a 


M 


(53) 


If  the  frequency  range  being  processed  by  the  optimum  detector  coincides 
with  0  S  co  <  co^  >  Eq.  (32)  can  be  written  in  the  following  form  |  for 


N(co)  =  N  ,  0  <  o  <  co 

v  '  o’  a 
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1 -  1 

-  1 

2 

A(DC  output)  S  _ 

£  Jr 

S(co) 

do 

D( output)  Nq  i 

no 

* 

..  S  J 

The  ratio  of  the  two  figures  of  merit  is 


(5U) 


performance  index  of  optimum  detector 
performance  index  of  standard  detector 


(55) 


The  reason  for  terminating  the  processing  range  cf  the  optimum 
detector  at  co  is  almost  certainly  that  virtually  all  of  the  signal  power 

Cl 

is  concentrated  in  the  range  0  <  to  <  co  .  Thus  it  is  reasonable  to 

°  a 

assume  that 


|  S(co)  do  =  5 

J 

0 


(56) 


If  one  expresses  the  spectral  function  by 


s 


S(co)  =  +  F(») 

a 


then  Eq.  (£6)  demands 


co 


F(co)  dco  =  0 


It  follows  that 


co 


co 


js(co)] 


-.2 


dco  = 


<~2 


+  2  F(co)  +  F^(co) 

co  2  “a 

a 


do 


0) 

r a  <-2 

—  +  I  F‘  (co)  dco  >  - — 

co  “  co. 

a  a 


(57) 
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with  equality  if  and  only  if  5(co)  -  ~  , 


5 


When  5(co)  =  —  the  ratio  of  the  performance  indices  ^  Eq.  1  is 


unity.  In  general,  therefore. 


performance  index  of  optimum  detector  ->  ^  ($Q) 

performance  index  of  standard  detector 

Thus  if  the  signal  as  well  as  the  noise  has  properties  of  white  noise 
bandlimited  to  0  <  co  <  co^  the  optimum  detector  is  identical  in  performance 
with  the  standard  detector.  If  the  signal  has  different  spectral  properties 
the  optimum  detector  makes  use  of  these  properties  to  improve  performance 
while  the  standard  detector  is  only  sensitive  to  total  signal  power. 

2)  Harkov  Noise 

If  the  cross correlation  between  hydrophones  docs  not  vanish,  additional 
terns  in  the  quadruple  sum  of  Eq.  (I47)  must  be  considered.  Assuming  that 
the  noise  at  each  hydrophone  consists  primarily  of  far  field  noise,  a 
reasonable  approximation  to  the  noise  autocorrelation  function  might  be 

-  00  ItI 

p(t)  =  e  0  (59) 


is  the  frequency  (in  radians  per  second)  at  which  the  noise  spectrum  has 

fallen  to  half  of  its  low  frequency  value.  Where  numerical  computations 

are  called  for  the  numerical  value  co  =  2n  x  1000  =  9lj20  radians/second 

0 

will  be  used. 

The  problem  of  computing  q  (t)  from  p(x)  has  been  discussed  in 
Progress  Report  No.  1.  With  a  spherically  isotropic  noise  field  the 
required  relation  is 

yTih 


4ih 


(t) 


2t., 

lh 


p(r)  dx 


T-T.. 

ih 


(60) 


where  x.,  =  .  d. .  is  the  distance  between  the  i^1  and  h+i'  hydrophouou 

lh  c  xh  ”  1 

(a  non-negative  quantity)  and  c  is  the  velocity  of  sound  in  water. 

Substituting  Eq.  (^9)  into  Eq.  (60)  one  obtains  by  straightforward 


computation 


qih(T)  1 


CO  X., 

o  lh 


-CO  X 

1  -  e  01  cosh  co  x 
o 


,  -CO  I  X I 

1  0  .  , 

-  e  sinh  co  x . 

_  o  lh 

CO  X. , 

o  lh 


*Tih 


|Tl  ^Tih 


q^(x)  is  evidently  non-negative  for  all  values  of  x  and  x^.  Thus  all 
terms  in  the  fourfold  sum  of  Eq.  ( I4.7 )  are  non-negative.  It  follows  that 
noise  coherence  from  hydrophone  to  hydrophone  (at  least  of  the  type 
assumed  here)  can  only  reduce  the  figure  of  merit  relative  to  the  value 
given  by  Eq.  (1+9).  ^ 

The  computation  of  j~  q_^(x)  q^(x)  dx  is  tedious  but  straightforward. 

Assuming  x  ^  >  x_^  and  taking  advantage  of  the  fact  that  q^(x)  is  an  even 

function  of  x  one  obtains 
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After  integration  and  some  algebraic  manipulation  this  becomes 
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(63) 


For  numerical  values  likely  to  occur  in  practice  this  equation  can  be 
simplified  considerably.  If  the  minimum  distance  between  hydrophones  is 
2  feet  and  one  uses  c  =  b7^0  ft/ sec,  Tih  >  =  k.  21  x  10"4  sec.  for 

i  /  h.  With  co^  =  9^20  rad/sec,  co  t  -  3.97  lor  i  /  h.  Hence  the 
approximation 

CO  T. 

sinh  w  t.  "=  cosh  w  t,  **  le  °  ^  ,  i/b  (61+) 

o  ih  o  ih  2  3 


introduces  negligible  error.  Using  Eq.  ( 6U )  and  discarding  terms  of 
-2co  Tju 

order  c  and  smaller,  one  obtains 
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■when  x  =  x  this  reduces  to 
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Similar  expressions  must  now  bo  obtained  for  the  cases  i  =  h,  j  =  k  and 
i  ■  h,  j  /  k. 
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q^(^)  q^(^)  dT  =  2  I  q1±(T)  qjk(T)  dx 
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Results  for  Linear  Array  with  Equally  Spaced  Elements 

All  of  the  terms  of  the  fourfold  sum  in  Eq,  (1*7)  are  now  available  and 
the  evaluation  of  the  output  signal  to  noise  ratio  becomes  in  principle  a 
straightforward  computational  problem.  Unfortunately  the  number  of  terms 
grows  very  rapidly  with  M  and  manual  computation  quickly  becomes  impractical 
unless  some  simplifications  arc  made.  A  configuration  leading  to 
particularly  simple  results  is  a  linear  array  with  equally  spaced  elements. 
Here 


T 


ih 


|i  -  h  |  t 


c 


(69) 


where  is  the  time  required  by  a  sound  wave  to  traverse  the  distance 
between  adjacent  hydrophones. 

If  one  uses  the  notation 
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then  for  the  array  satisfying  Eq.  (69) 
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The  further  change  of  variable 
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leads  to  the  form 


M 

V"1 

Ci»h,3,k 

J 

i,j,h,k  =  1 


M-l 

M2  C  +  2M  S  (M 
■  0,0  Cj 


M-l 


s-1 


s>  co,s  +  Z'L (M  -  e)  °e,c 
1=1 


M-l  M-l 


+  u  y  (M  -  t)0/L  -  s)  C 


y  S 


1=1  S=1 


M-l 

M2  C  +  U  M  V  (M 

0,0  <Lj 

s=l 


M-l  M-l 


s)  c  +  ii  -V 

O  ,  S  i _ I 


)  (M  -  £)(M  -  s)C 


lfs 


1=1  S=1 


(75) 


The  coefficients  e  and  C^,  g  are  given  by  Eqs.  (67),  (68),  and  (65) 


o,o  o,s 


respectively  (with  x =  sr  ,  t.,  =  lx  ). 

^  jk  o’  ih  o 

A  further  considerable  simplification  can  be  made  if  we  recognize  that 

for  the  numerical  values  introduced  on  page  22  (co^  =  9^20  rad./ sec. , 

-co  x 

hydrophone  spacing  x^  =  2  ft.),  cc^t^  =  3.97.  Hence  e  0  0  =  0. Ol89j 
so  that  the  erqponential  terms  in  the  definitions  of  the  coefficients  g 


can  be  omitted  with  only  very  minor  error.  Then 
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Hence  the  figure  of  merit  of  Eq.  ( U7 )  is 
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(80) 


1 — 1  lOJ 


Discussion 


a)  If  to  t  »>  1,  C.  "=  0  unless  s  =  Z  =  0.  In  this  cane  (noise 

OO  HjS 

components  at  different  hydrophones  uncorrelated)  Eq.  C 80 )  reduces  to  the 
simple  form 


This  expression  is  similar  in  form  to  Eq.  (52),  but  differs  by  a  numerical 
constant  because  the  noise  spectrum  has  a  different  shape  than  that  assumed 
in  the  earlier  computation. 

b)  If  =  3.97,  =  9U20  rad./ sec. ,  as  assumed  earlier,  and  M  =  5, 

the  figure  of  merit  becomes 


129  —It 
N 


(82) 


If  signal  and  noise  spectra  have  the  same  shape,  the  output  signal 
to  noise  ratio  for  the  optimum  detector  using  the  same  5~element  linear 
array  with  2-foot  hydrophone  spacing  is  from  Eq.  (30) : 


Bryn  has  plotted  the  array  gain  for  a  five-element  equally  spaced 
linear  array.  Adapting  his  plot  to  the  parameters  used  in  this  computation 
one  finds  that  G(co)  has  essentially  reached  its  maximum  value  of  5  at 
to  =  7500  rad./ sec.  Making  a  rough  approximation  to  the  curve  below  that 
frequency,  one  obtains  the  following  numerical  expression  for  output 
signal  to  noise  ratio  in  terms  of  the  frequency  top 
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T_  2. 

2n 


I/25  -  5.  55  x  103 


A  (DC  output) 
D( output ) 


If  w2  =  2n  x  5000  =  31*000  rad. /sec.,  this  reduces  to 


(8U) 


A  (DC  output)  ^  _£L  \|t 

D(  output )  ^ 


(8^) 


If  w  =  2n  x  10,000  =  62,800  rad,/ sec.  ,  the  corresponding  value  is 

A(DC  output)  ~  A_  fr 
D(  output)  N 


(86) 


Thus  the  rms  output  signal  to  noise  ratio  of  the  optimal  detector  is  better 
than  that  of  the  standard  detector  ty  a  factor  of  about  2.6  if  the  range 
of  frequencies  processed  by  the  optimum  detector  extends  to  5000  cps  and 
by  a  factor  of  about  3.8  if  the  processed  range  extends  to  10,000  cps.3" 
Which,  if  either,  of  these  figures  is  realistic  depends  on  the  frequency 
range  over  which  the  assumed  equality  in  shape  of  signal  and  noise  spectra 
corresponds  reasonably  well  to  the  actual  facts. 

c)  In  order  to  remove  the  somewhat  arbitrary  upper  bound  of  the 
frequency  range  processed  by  the  optimum  detector,  an  attempt  will  be 
made  to  gain  more  insight  into  the  probably  high  frequency  behavior  of  the 
signal  spectrum.  One  factor  that  should  undoubtedly  be  considered  is  the 
frequency  dependence  of  the  transmission  characteristics  of  water.  Using 
a  model  discussed  ty  Fckart  one  oan  write  the  following  expression  for 
the  signal  spectrum  received  at  each  hydrophone 

^IJo  band  limitation  was  assumed  in  the  standard  detector. 

2 

C,  Lckart,  Optimal  Rectifier  Systems  for  the  Detection  of  Steacfy 
Signals.  University  of  California,  Marine  Physics  Laboratory  of  the 
Scripps  Institution  of  Oceanography,  SIO  Reference  52-11,  U  March  1952, 
p.  13,  Eq.  (18). 
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Received  spectrum  *  — ^  e  S(«) 

r 

S(co)  is  the  signal  spectrum  at  a  point  close  to  the  target  and  r  is  the 
target  range  in  yards.  Assuming  S(co)  and  N(co)  identical  in  shape,  cne 
obtains  from  Eq.  (29) 
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[g(co)]  2  do  (88) 


C  =  J  S(co)  dw  is  the  signal  power  at  the  reference  point.  Evaluation  of  the 
0 

integral  is  cumbersome  because  of  frequency  dependence  of  G(«).  Hwever, 
since  G(co)  only  varies  from  a  minimum  of  3.55  near  »  =  0  to  a  maximum  of  5 
above  «  =  7500  rad. /sec.  ,  a  reasonable  approximation  can  be  obtained  by 
treating  it  as  a  constant.  Then 


A  (DC  output)  ^  $  <rzr  n  6.1*1  x  10^ 
D(  output)  $  r9^ 


For  G  =  JU  this  is 


A  (DC  output)  ~  2.56  x  109  r^- 

D(  output)  r'9^  ^ 


For  G  =>  5  the  value  is 


A  (DC  output)  -y  3 . 20  x  109 
D(  output)  -N 

Except  at  extreme  ranges  (above  10^  yards)  Eq.  (91)  should  be  th< 
better  approximation  and  will  be  used  in  further  comparisons. 
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For  the  standara  detector  the  shape  of  the  signal  spectrum  is 
unimportant;  only  the  total  signal  power  matters.  Hence  the  factor 

e  1X0  can  be  ignored  if  it  is  essentially  equal  to  unity  over  the 

effective  bandwidth  of  the  signal  spectrum.  This  is  true  if 


or 
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Then  the  figure  of  merit  for  the  standard  detector  becomes  ^from  Eq.  (82) | 
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The  ratio  of  Eqs.  (9l)  and  (93)  assumes  the  form 


rms  signal  to  noise  ratio  of  optimal  detector  _  2h.8  _  .  in5 

J""’  "  r  "  ■*  ‘  -  -  -  -  «r  n  i  iin-  m  t  IT7T  J  27  <^<C  JL  #  0^  X  1U 

rms  signal  to  noise  ratio  of  standard  detector  r' 

(9U) 

Thus  even  at  a  range  of  100  yards  the  possible  improvement  due  to  the  use 
of  optimal  detection  procedures  is  limited  to  a  factor  of  somewhat  less 
than  8.  (If  the  assumption  of  low  signal  to  noise  ratio  remains  valid.  ) 

V.  Concluding  Remarks 

The  report  has  compared  the  performance  of  a  likelihood  ratio  detector 
with  that  of  a  standard  detector  on  the  basis  of  output  signal  to  noise  ratio. 
Under  the  conditions  investigated,  the  improvement  attainable  through  use  of 
optimal  (likelihood  ratio)  techniques  was  in  most  cases  only  moderate.  It  may 
therefore  be  important  to  review  once  more  the  basic  assumptions  underlying 
the  analysis. 
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1)  Signal  and  noise  were  assumed  to  be  independent  Gaussian  random 
processes. 

2)  The  signal  to  noise  ratio  at  each  hydrophone  was  assumed  to  be 
small  compared  to  unity. 

3)  Signal  and  noise  processes  were  assumed  to  be  stationary  and  the 
signal  to  noise  ratio  at  each  hydrophene  was  assumed  to  be  the  same. 

4)  Numerical  calculations  were  concerned  primarily  with  cases  of 
relatively  broadband  signals  and  noises,  mostly  with  similar  or 
identical  power  spectra. 

5)  The  array  of  hydrophones  was  considered  fixed.  Numerical 
calculations  dealt  particularly  with  a  linear  array  of  equally 
spaced  elements. 

Assumption  l)  is  particularly  important.  It  implies  that  signal  and 
noise  can  be  characterized  only  by  their  total  power  and  spectral  properties. 
Of  these  two  features  the  standard  detector  uses  only  total  power,  hence  it 
is  in  general  inferior  to  the  optimum  detector.  However,  the  report 
indicates  that  (unless  the  spectra  of  signal  and  noise  are  drastically 
different)  the  benefit  to  be  derived  from  use  of  spectral  properties  is 
quite  limited.  Further  work  is  planned  to  investigate  the  improvements 
possible  with  optimal  detection  schemes  when  there  are  pronounced 
differences  in  signal  and  noise  spectra.  Signals  containing  narrowband 
components  at  known  or  unknown  frequencies  would  be  a  case  of  particular 
interest. 

The  analytical  treatment  of  non-Gaussian  random  processes  tends  to  be 
very  difficult,  and  an  attempt  to  solve  the  problem  under  such  assumptions 
does  not  appear  promising.  However,  if  information  concerning  waveshapes 
of  expected  signals  were  attainable,  great  improvements  in  detector 
performance  might  be  possible.  The  analytical  tools  for  handling  this  case 
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are  also  available.  An  effort  to  search  out  features  of  the  expected 
signals  related  to  waveshape  rather  than  spectral  properties  therefore 
appears  indicated. 

Assumptions  3)  and  5)  suggest  additional  possibilities  for  further 

investigations.  It  has  been  shown,  for  instance,  that  clipping  can  improve 

the  performance  of  a  power  detector  when  the  noise  power  varies  from 

1  2 

hydrophone  to  hydrophone,  or  as  a  function  of  time.  It  would  be 
interesting  to  investigate  whether  the  optimum  detector  also  performs  an 
operation  analogous  to  clipping  under  these  conditions.  If  hydrophone 
locations  are  regarded  as  an  additional  set  of  adjustable  parameters  the 
optimization  procedure  can,  in  principle  at  least,  be  carried  one  step 
further.  One  could  then  perhaps  inquire  into  the  best  location  for  a  fixed 
number  of  hydrophones  in  a  given  volume  and  the  performance  index  of  the 
associated  optimal  detector.  It  appears  unlikely  that  a  very  general 
solution  to  this  problem  would  be  computationally  feasible,  but  it  would  be 
of  interest  to  gain  even  limited  insight  into  the  dependence  of  system 
performance  on  hydrophone  location. 

Finally  the  theory  of  optimal  detection  may  provide  a  means  for 
assessing  realistic  costs  of  such  signal  processing  techniques  as  clipping. 

One  could,  for  instance,  compare  the  performance  of  an  optimal  detector 
operating  on  clipped  hydrophone  data  with  the  performance  of  an  optimal 
detector  operating  on  unclipped  data.  Whether  such  a  comparison  is 
computationally  manageable  for  cases  of  practical  interest  has  not  yet  been 
investigated. 

^T.  Usher,  Signal  Detection  by  Arrays  in  Noise  Fields  with  Local  Variations. 
Progress  Report  No.  2.  Electric  Boat  Co.  Research.  Yale  University,  March  I963. 
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J.  B.  Thomas,  T.  R.  Williams,  On  the  Detection  of  Signals  in  Ncn- 
Stationary  Noise  by  Product  Arrays.  J.  Acoust.  Soc.  of  America.  Vol.  3k,  No.U, 
pp.  U53-U62,  April  1959. 
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I .  Introduction 

The  signal  detection  characteristics, of  local  arrays  of  transducers 

12  3 

have  been  extensively  analyzed  *  *  for  the  situation  in  which: 

1.  The  output  of  each  transducer  is  subject  to  infinite  clipping 
or  purely  linear  processing; 

2.  The  processed  output  signal  from  each  transducer  is  summed; 

3.  The  summed  result  is  applied  to  a  square-law  detector; 

U.  The  squared  sum  is  filtered. 

However,  the  effect  of  a  general  nonlinear  processor  following  each 
transducer  relative  to  signal  detection  efficiency  has  not  been  determined. 
Neither  has  the  effect  of  a  general  non-square  law  detector  been 
examined.  It  is  the  purpose  of  this  report  to  investigate  both  effects. 

The  two  performance  indices  defined  in  Reference  3  are  used  to  provide 
a  numerical  indication  of  the  signal  detection  efficiency. 

The  signal  response  ratio,  R, ,  is  defined  to  be  the  ratio  of  the 
normalized  average  signal  output  of  the  system  with  either  the  nonlinear 
processor  or  the  non-square  law  detector  to  the  normalized  average  signal 
output  for  the  system  with  the  linear  processor  and  square-law  detector. 

Similarly,  the  index,  R2,  is  defined  to  be  the  ratio  of  the  square 
root  of  the  signal-to-noise  ratio  of  the  directivity  pattern  for  the 
system  with  either  the  nonlinear  processor  or  the  non-square  law  detector 
to  that  for  the  system  with  the  linear  processor  and  square-law  detector. 
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The  major  purpose  cf  the  analysis  is  to  determine  optimum  processor 
functions  and  optimum  detector  functions  for  maximizing  both  response 
indices.  Another  consequence  of  the  analysis  is  a  derivation  of  a  general, 
expression  giving  the  autocorrelation  function  for  the  output  of  a 
nonlinear  device  in  terms  of  the  input  autocorrelation  function.  The 
major  assumptions  of  the  analysis  are  listed  in  the  following  section. 
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II,  Assumptions  and  Definitions 


The  system  analyzed  in  the  following  sections  is  shown  in  Fig.  1. 


Transducers  Processors 


Fig.  1  Array  with  Nonlinear  Processors  and  Non-Square  Law  Detector 

Following  are  the  major  assumptions  inherent  in  the  analysis: 

1.  The  transducers  in  the  array  have  omnidirectional  characteristics, 

2.  The  autocorrelation  functions  of  both  the  signal  and  the  noise 
inputs  to  the  processors  are  identical.  The  signal  power  input 

and  the  noise  power  input  to  each  processor  do  not  vary  from  channel 
to  channel* 

3*  The  cross -correlation  between  noise  inputs  for  different  charnels 
*  is  zero. 

The  signal  and  the  noise  inputs  to  the  processors  are  Gaussian 
and  stationary.  Signal  and  noise  are  independent. 
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5.  The  processor  is  an  odd-function,  zero-memory  device, 

6,  The  detector  is  an  even-function  device. 
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III.  Correlation  Between  the  Output  Signals  of  Identical  Nonlinear  Devices 

Since  the  analysis  of  array  performance  in  following  sections  requires 
an  evaluation  of  the  cross-correlation  of  the  output  signals  from  two 
identical  nonlinear  devices,  or  the  evaluation  of  the  autocorrelation  for 
the  output  signal  of  a  nonlinear  device,  a  general  derivation  applicable 
to  both  situations  is  made  in  this  section. 

It  is  assumed  that  the  input  variables  are  Gaussian,  with  zero  mean, 
and  equal  variances.  Furthermore,  the  two-dimensional  probability  density 
function  is  given  by 


f(ere2)  — Tr=r exp  [' ei 


2  2 
**■  -  2ae 


2 no  Yl  -  a 


1=2] 


The  characteristic  function  of  f  in  'tenr*s  cf  the  complex 


variables  w^  and  w^  is 


00  00 

$(Vl,w2)  =  j  j  expj^e^  +  o2w2  f  (e^e^  de-jde. 


.£C  _00 


0,2  2  „  x 

exp  (w^  +  w2  +  aaw^J 


It  is  desired  to  find  the  statistical  average  of  g(e^)g(e2)  where 
g(e)  is  a  general  nonlinear  function  with  zero  memory.  The  development 
follows  the  transform  method  used  in  Chapter  13  of  Reference  U.  The 
correlation  function  g(e^ )g(e2 )  in  its  final  form  is  expressed  in  terms  of 
a  power  series  in  a,  the  normalized  correlation  between  e^  and  e2> 

It  is  convenient  to  define  the  function  g(e)  in  two  parts: 


g+(e) 


e  JO 
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e  <  0 


and 


g  (e) 


e  >  0 


0  <  0 


Then 


g(e)  -  g+(e)  +  g_(e ) 


The  Laplace  transforms  for  the  two  parts  of  g(e)  are 


G+(w)  =  I  g+(e) 


-we 


de 


-we 

G  (w )  =  j  g  (e )  e  de 


-OO 


The  inverse  relations  are  found  in  Eqs.  (8)  and  (?)« 


(U) 

(5) 

(6) 

(7) 


g+  (e) 


g  (e) 


»+(w) 


ew 

e  dw 


ew 

e  dw 


(8) 

(9) 


The  required  contours  of  integration  are  shown  in  Fig.  2.  The 
infinitesimal  displacement  A  in  each  case  is  required  for  convergence 
of  the  transforms.  It  is  assumed  that  g(e)  dees  not  increase  exponentially. 


Fig.  2  Contours  of  Integration  in  the  Complex  w-Plane 
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The  correlation  between  g(e^)  and  g(e2)  is 


00  00 


g(e1)  g(e2) 


g  (©-L )  g(e2)  f(e1,e2)  de^eg  (10 ) 


~P0  _0C 


Substitution  of  the  results  of  Eqs.  (5),  (8)  and  (9)  into  Eq0  (10)  gives 


oo  oo 


g(en  )  g(o  )  =  — *  I  I  |G  (w  )  e  1  1  dw  +  f  G  (w  )  e  1  1  dw 

1  2  (2"3)2J  J  Jo  1  D  J 


.jac  _co 


epW?  r  e  w 

G+(w2)  e  ^  dw2  +  J  G_(w2)  e  dw, 

D 


f(e1,e2)  de^eg 

(11) 


With  the  application  of  the  definition  in  Eq,  (2),  Eq„  (11)  becomes 


G+  (w]_ )  G+  (w2  )$  (w]_,w2  )  ^1^2 


+  J  J  G+  (w ]_ )  G_  (w2 )  ^  w2  )  dw^dw,, 

C  D 

^  ^G_(wl)  G+  (w2  )  ^(wi>w2  )  dw-jdifg 


D  C 


+  j  j  G_(w. i)  G_(w2)  ^(wi3w2')  ^1^2 


D  D 


(12) 


If  the  infinite  series  for  expj^a^a  w^w2  is  utilized  for  the  corresponding 
term  in  (|(wpW2),  the  double  integrals  in  Eq,  (12)  become  separable 
products  of  integrals. 


as 
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Since 


o°  2n  n  n  n 
a  a  w,  h„ 


r  2  1  V"  1 

exp[a  3  V2J "  L - ~ 

n-0  n‘ 


then, 


2  2 


00 


r — i  2n  n  _ 

g(e-  7  g(®2 ^  “  Y>  ,2in"  7 - 7 

-  2  £o  (2"J> 


" 

f  1 

1 

WjH  G+ (w1 )  c 

dwl  1 

.( 

) 

c 

(13) 


2  2 
a  w2 


w2n  G+  (w2)  e  ‘  dw2 


2  2 
a  w^ 

2  . 


2  2 
a  w2 


Win  G+(w]l)  e  *w1  w2n  G_(w2)  e 


dw„ 


D 


2  2 
a  w^ 

T 


2  2 
a  w2 


+  w^nG^(w1)  6  dw^  jv2n  G+(w2  )  6 


dw. 


D 


2  2 
a  w^ 


2  2 
a  w2 


dw. 


|w. i”  G_(w1)  e  dw1  |w2n  G_(w2)  e 
D  D 


(lU) 


Dae  to  the  product  form  of  Eq.  (Hi),  it  may  be  rewritten  as 

2 

y 

7 


r — r  2n  n  n 
V  Qf  a  1 


n-0 


2  2 
a  w 


2  2 
CT  W 


wn  G+(w)  e  dw  +  Jwn  G  (w)  e  dw 
LC  D 


(15) 


It  may  be  shown  that 


-we 


g+(e)  e  de  =■  w  G+(w)  -  g(0) 


(16) 


and 


-we 


g  (e  )  e  de  =  w  G  (w )  +  g  (0 ) 


(17) 
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_QO 


whore 


g  (e) 


d© 


(18) 


Furthermore,  the  relations  in  Eqs..,  (16)  and  (17 )  may  be  extended  for  the 
derivative  to  yield 
00 

f  r«1  “we  n  o  (n-m)  Lm-l] 

I  g^1  (e )  e  de  -  w  G+(w)  -  ^  w  g  (0) 

i  m«l 


and 


(19) 


■JOO 


r  ,  -we  (n-m )  [m-l] 

g^e)  e  de  =  wn  G  (w)  +  ^  w  g  (o) 

m=l 


(20) 


wnere 


g>rJ(al  - 

da 


n 


(21) 


When  the  results  of  Eqs.  (19)  and  (20)  are  substituted  into  Eq.  (15>), 
Eq.  (22 )  is  obtained. 


CO 


r — i  2n  n  .. 

«(vg(va  E  ^77-7 

n-0  ^  > 


2  2  » 
o  w  r 


-we 


|dw  e  ‘  I  g|n^(e)  e  de 


Lc 

2  2 


J 

0 


,  22  ,  220 

A  f  (n-m )  [m-l]  2*-  f  2-g-  C  .  , 

w  g  (0)  e  dw  +  dw  e  2  j 

mJl  C  D  i 


-we 

glilJ(e )  e  de 


^  [  (n-m )  [m-l) 

-  )  jw  g  (0)  e  dw 


m=l 


D 


(22) 
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Cm— l] 

It  can  be  Bhown  that  the  terms  involving  g  (0)  cancel  (See 
Appendix  A)  because  the  difference  of  the  integrals  over  contours  C  and  D 
may  be  made  as  small  as  desired  by  making  A  in  Fig,  2  arbitrarily  small. 
The  inverse  transform  of  the  single -dimensional  characteristic 
function,  yields  the  Gaussian  probability  density  function: 


2  2 
o  w 


1 

2nj 


2 

ie  2  e  dw  -  — e  2®~  «  f  (e  ) 


-we 


"V2n  a 


(23) 


C» 


The  contour  C*  is  the  imaginary  axis  of  the  Y7-planec  Since  C  and  D  may  be 
made  arbitrarily  close  to  the  imaginary  axis  of  the  w-plane,  Eq.  (22) 
becomes 


r — I  cn 

e(ei)  g(V  ■  Z 


2n  n 
a 


n=0 


,w 


g+“J(e)  f  (e  )  ds  +  I  g^(e)  f(e) 


de 


-O O 


(24) 


Finally, 


g(ex)  g(e2  )  -  ^  ~J  gtn)(e)  a° 


(25) 


n=0 


In  Eq,  (25)  the  bar  indicates  an  average  with  Gaussian  statistics 
for  e.  An  alternate  form  of  Eq,  (25)  is  possible  since  it  can  be  shown 
(See  Appendix  B)  that 


g[n\e)  =  a"n  Hn(^)g(e) 


(26) 


In  Eq,  (26)  Hnj~j)  is  the  Hermite  polynomial  of  order  n  with  ^  as  the 
argument. 

Then 

_ 2 

iT5TraV-  Z  ^iH„lf]g(e)  a"  <27) 

'  n=0 
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The  results  in  Eq.  (25)  and  (27)  can  be  used  in  determining  the  cross 
correlation  between  the  output  variables  of  identical  nonlinear  processors 
In  this  case,  a  is  the  normalized  cross-correlation  between  the  input 
variables  of  the  processors.  The  result  can  also  be  used  to  determine 
the  autocorrelation  function  of  the  output  of  a  nonlinear  processor.  In 
the  latter  case,  e^  and  e^  represent  the  input  variable  at  two  different 
times  and  a  is  the  normalized  autocorrelation  function  of  the  input 
variable , 
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IV.  Average  Signal  Output-Nonlinear  Processors  and  Square  Law  Detector 

The  symbols  and  definitions  in  Reference  3  will  be  followed  aa  closely 
as  possible.  For  the  system  with  the  nonlinear  processors,  the  directivity 
pattern  is 

M  M  _ _ 

y(Q)  -  Yj  8{8i(t)  +  + 

i-i  3-1 

For  the  input  signals  to  the  processors, 

a2  -  S  +  N  (29) 


The  normalized  cross-correlation  a  is 
[s.jJtJ  +  ni(t  (_s  j  (t )  +  nj(t)] 

fl  a  m  1  ■  .  i  --T— i  *  i  i  ■  ii  ii  i  i  ii  ii-  ■  -  -  n  ■' 

S  +  N 


S  +  N  p^Ti^ 


i  /  3 

i  =  3 


(30) 


In  Eqs.  (29)  and  (30),  p  is  the  normalized  autocorrelation  function 
for  the  signal,  S  is  the  total  signal  power,  and  N  is  the  total  noise  power. 
The  depend  on  the  spatial  time  delays  and  electrical  time  delays 

associated  with  the  typical  steered  array. 

The  results  in  Eqs.  (2$)  and  (27)  may  be  used  to  simplify  Eq.  (28). 

It  should  be  noted,  however,  that  when  a  3  1,  Eq.  (2$)  reduces  to  one  term, 
the  mean  square  value  of  g(e).  Furthermore,  the  processor  function  in  the 
array  has  been  assumed  to  have  odd  symmetry,  so  that  only  odd  order  terms 
in  Eqs.  (25)  and  (27)  are  present. 

Using  Eq.  (25),  we  have 


_ M 

y(e)  "Hg(e)2  +  2  ^ 

i-1 


M  oo 

Z  I 


j=i+l  n=0 


(2n+l ) 
S 

(2n  +  1)J 


C2n+1]  2  2n+l 

g  (e)  P  (t±.) 


(31) 
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The  "on  target”  response  occurs  when  all  ■  0,  and  the  "off  target” 
response  occurs  when  the  are  all  large  enough  so  that  the 
approach  zero.  These  responses  are  found  in  Eqs.  (32)  and  (33) 


respectively. 


70  ■  M  g.>)2  ♦  (M2  -  M)  V  - £»*»(.)* 

to  <2”  *  1>J 


7l  a  M  g(e) 


The  normalized  signs!  output  Ay  is 


4y’.  -  (M  .  1)  S=2 


y  s(2ntl>7^M2 

U  (2n  +  1)J 


For  the  system  with  linear  processors,  g(e)  =  Ke  and 


yA0  -  K2 (MN  +  SM2 ) 

yAL  -  K2M(N  +  S) 

AyA  -  (M  -  1)  — 

A  S  +  N 


The  signal  response  ratio,  R^,  has  previously  been  defined  to  be  the 
ratio  of  Ay  to  Ay^.  _ _ 

?\IsE±I 

T,  Ay  Is  Hr  N\  na0  (2n  +  1)J _ 
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An  alternate  form  of  is  obtained  by  using  Eq.  (26): 

(39) 

From  Eq.  (39)  it  is  evident  that  only  the  first  tern  of  the  series  is 
important  for  small  values  of  input  signal  to  noise  ratio.  Then 

for  |  <  <  1  (1*0  ) 

It  is  interesting  to  note  from  Eq0  (38)  that  for  large  values  of 
signal  to  noise  ratio 

S  - 

-  °  a  — )  1 

S  +  N 

2 

and  the  sum  in  the  numerator  cf  Eqa  (38)  approaches  g(e)  ,  Thus  the  signal 
response  ratio  approaches  unity.  Under  strong  signal  conditions,  then,  the 
average  signal  performance  of  the  system  with  nonlinear  processors  approaches 
that  of  the  system  with  linear  processors. 

Also  if  one  considers  the  non-negative  function 

(X  g(e)  +  e)2  >  0  (1*1) 

where  X  is  an  arbitrary  parameter,  we  have  by  expansion, 

l^g(e)2  +  2X  eg'(e)  +  o'  >  0  (1*2  ) 

Since  there  can  be  no  real,  unequal  roots  for  X, 

eg(e )‘  -  a  g(e)2  <  0  (1*3) 
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Applying  the  result  in  Eq„  (U3)  to  the  small  signal  approximation  for 
R^,  wo  have 

\  <  1  ^r  |  <  <  1  (tU) 

Since  all  terms  in  the  series  expression  for  R^  are  positive,  unity 
Js  an  upper  bound  for  R^  for  all  ratios  of  signal  to  noise.  Furthermore, 
the  equality  in  Eqs.  (Ill)  through  (U3)  holds  only  if  g(e)  is  a  linear  function 
of  e  (i.e.,  g(e)  =  ko ),  so  that  the  signal  response,  ratio  is  less  than  unity 
for  all  finite  input  signal-to-noise  ratios  for  all  nonlinear  processors. 

Thus  the  linear  processor  is  optimum  in  the  sense  of  producing  maximum 
average  output  signal  for  the  array. 
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V.  Signal -to -Noise  Ratio  at  Output  of  Array — Nonlinear  Processors  and 
Square  Law  Detector 

The  goal  of  this  section  is  the  development  of  an  expression  for  the 
performance  index,  Rg,  which  is  the  ratio  of  the  square  root  of  the  signal- 
to-noise  ratio  for  the  system  with  nonlinear  processors  to  that  for  the 
system  with  linear  processors.  Since  the  expressions  for  average  signal 
output  were  developed  in  the  previous  section,  most  of  the  analytical 
development  in  this  section  will  be  concerned  with  the  evaluation  of  the 
variance  of  the  output  of  the  low -pass  filter  in  Fig.  1. 

In  general 

Ry('T)  =  E[y(t)  y(t  +  t)]  -  y2  +  Ry'(^).  (U5) 

where  y  represents  the  average  array  output,  and  R  (x)  is  the  autocorrelation 
function  for  the  variational  component  of  y(t).  If 

HZ  ( 3®)  3  - - — r  0*6 ) 

1  +  3  — 

“F 

from  Eq.  (56),  Reference  3,  we  have 

00 

otj2  ~  oipj' Ry' (t)  dx  (dip  small)  (1+7 ) 

0 

For  the  analysis  of  array  performance  with  nonlinear  processors,  it 
will  be  assumed  that  the  input  signal-to-noise  ratio  is  sufficiently  small 
so  that  only  the  input  noise  contributes  to  noise  ( a ^  )  at  the  output. 

This  assumption  is  necessary  for  making  the  analysis  tractable,  and  also 
is  not  unduly  restrictive,  since  only  small  signal-to-noise  ratios  produce 
any  detection  problems. 
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Thus  the  autocorrelation  function  for  y(t)  becomes 


R  (t)  =  E  x^t)  x^t  +  t) 

y 


The  evaluation  of  Eq.  (1*8)  is  helped  by  separating  the  terms  into  the 
different  categories  as  shown  in  Table  1. 


Subscript 

relations 

Number 
of  terms 

Evaluation  and  commen'os 

i  -  j  and  k  =  l 

M2 

2 

i  f  k  yields  M  -  M  terms. 

E[g2(ni.)  33  g2^nk^  =  g2 

1  =  k  yields  M  terms. 

E  g2|r^(t)|  g2{ni(t  +  t)} 

i  -  j  and  k  ft 

i  /  j  and  k  =  l 

2  (M3  -  M2) 

Helds 

E  g2(ni)  gCr^)  g(n^)J  -  0 
because  n^  and  n^  are  uncorrelated 

and  g(n)  is  odd. 

i  /  j  and  k  /  l 

-  2M3  +  M2 

Only  terms  of  the  type  where  i  =  k, 

j  =  and  i  *  l ,  j  -  k  are  non-zero. 

2(M‘_  -  M)  in  number) 

Table  1  Classification  of  Terms  in  Eq,  (1*8) 
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From  Table  1,  Eq.  (U6 )  becomes 


R  (t) 

y 


(M2  -  M) 


g2(o) 


+  M  e  (t )  j  g2i  e  (t  +  t)  \+  2  (M*"  - 


m)  ipiy  +  t^2 
(U?) 


Both  the  second  and  third  terms  of  Eq„  (U9)  may  be  evaluated  using  the 

2 

result  in  Eq.  (26).  For  the  second  term,  since  g  (e )  is  an  even  function, 
•we  have 


For  the  third  term, 


(5D 


kk-lB 


Substituting  the  results  from  Eqs.  (50 )  and  ($l)  into  Eq.  (U9 )  and 
collecting  terms,  we  have 


R  (t)  = 
7 


/ 

2 

M  g 

+  M 

I 

1  „  2  2  1  „  2  U 

pT  H2g  P  +  HI  HUg  P  +  *  *  '* 


+  2 


(M2  -  M) 


rr—h  2  2  n— 2  75 — 2  h 

Hlg  P  +  37  Hlg  H3g  p  + 


(52) 


For  large  M,  the  terms  multiplying  M  in  Equ  (52 )  may  be  neglected  in 

2 

comparison  with  those  multiplying  2(M  -  M),  since  in  general  the  order 

2  ^ 

of  magnitude  of  internal  terms  multi plying  p  n  is  generally  the  same#  Also 

_  2 

the  term  (Mg2  )  in  Eq.  (52)  is  simply  the  y  term  in  Eq.  (U5 ) «.  From 

Eq.  (U7), 

2  „  /„„2 


00 

oo 

f- 

-  M) 

sighl 

p2dT  +  jj  Hjg2  1 

f  u 

p  dT  +  ... 

J 

C 

J 

)  c 

) 

(53) 


S  <<i 

N  M 

For  the  linear  system  for  small  signals,  from  Eq,  (U9),  Reference  3, 
we  have  „ 


ctza  =  2“f 


M2  N 2  (  p2dx 


S  <<i 

N  M 


<5U) 


The  definition  of  R^  yields 


y0  *  yL 


~ZA 

•  _ 

AO  "  J,AL  CZ 


(55) 


Substitution  of  the  results  in  Eqs.  (32  ),  (33),  (35),  (36),  (53)  and  (5U) 


for  small  results  in  the  following  expression. 
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In  Eq.  (56)  K,  is  defined  to  be 

OO 

A2n  3  J P2U(T)  dT  (57) 

0 


The  alternative  form  of  Eq.  (56)  utilizing  derivatives  of  g(e)  is 


An  inspection  of  Eq.  (58)  reveals  that  R ^  is  less  than  unity  for  all 
g(e)  for  which 

H2n*lll)g(e)  ■  °2ml  g2nt3)(e)  /  0  for  n  >  1  (59) 

The  only  function  excluded  from  the  inequality  in  Eq,  (59)  is  the  linear 
function  g(e)  ®  ke,  for  which 

H1(-“)g(e )  -  a  gtlj(e)  -  ok  (60) 

and 

H2n+l!  o)g  ^  =  a^2n+1^  i2n+1\e)  =  0  for  n  >  1  (6l) 

Thus  the  linear  processor  is  optimum  because  it  produces  a  signal-to-noise 
ratio  at  the  array  output  which  is  greater  than  the  signal-to-noise  ratio  for 
any  nonlinear  processor. 

In  the  following  section,  the  general  results  in  Sections  IV  and  V  are 
evaluated  for  several  types  of  nonlinear  processors. 
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VI. 


Evaluation  of  R^  and  R2  for  Specific  Nonlinear  Processors 

A  wide  class  of  nonlinear  processors  is  described  by  the  general  power 
law  relation  found  in  Eq.  (62).  Note  that  g(e)  is  an  odd  function  of  e. 

for  e  >  0 

for  e  =  0  (62) 

for  e  <  0 


g(e )  -<  0 

-Kiel” 


From  Reference  6,  page  20,  it  is  given  that 

00  2 

where  (k)j  is  the  general  factorial  function  (k  is  not  necessarily  an  integer )„ 
The  averages  of  the  functions  found  in  Eq,  (UO)  are  calculated  below. 
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(65) 


The  performance  index  R^  becomes 
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The  last  term  in  Eq.  (66)  is  evident  because  of  the  identity  (see  Reference 

6,  p.  18) 
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The  result  in  Eq.  (66)  is  plotted  in  Fig.  3e  It  can  be  seen  that  the 
maximum  value  of  R^  occurs  for  m  ■  1*  However,  the  peak  is  sufficiently 
broad  so  that  the  value  of  m  may  vary  between  zero  and  three  so  that  the 
value  of  R^  >  .6.  -Note  that  the  value  m  -  0  yields  the  function  which 
describes  the  infinite  clipper. 

Next,  the  performance  index  Rg  found  in. Eq,  (58 )  will  be  evaluated 
for  the  odd-function  power-law  processor.  The  averaged  derivatives  in 
Eq.  (08)  are 
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(70) 


Using 

one  obtains  from  Eq,  (57 ) 


— to  It  | 

p(x)  =  e  0 


(71) 


(72) 


Thus  Rg  becomes 


R2(m) 


1  +  -  l)2  +  -  1)2(">  -  3) 


180' 
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(73) 
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Fig.  3  Performance  Indices  for  Odd-Function  Power-Law  Processors 
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Equation  (73)  is  evaluated  and  plotted  as  a  function  of  m  in  Fig.  3.  It  may 
be  noted  that  the  series  terminates  when  m  is  an  odd  integer,  and  converges 
rapidly  when  m  is  reasonably  small. 

The  maximum  value  of  is  seen  to  occur  for  m  -  1,  as  expected.  The 
peak  of  the  function,  however,  is  again  very  broad. 

Another  processor  function  which  has  evoked  some  interest  is  the 
infinite  clipper  with  a  dead  zone  at  high  amplitude.  The  function  is  shown 
in  Fig.  U. 


Fig.  U  Infinite  Clipper  with  High  Amplitude  Dead  Zone 


In  order  to  evaluate  in  Eq.  (UO)  we  have 


e  g^e 

g(ef 


n 

In  Eq.  (75),  the  normal  probability  function  P(X)  is  defined  by1 


(7U) 

(75) 

(76) 
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Thus  the  index  is 


ya) 


P(a) 


For  small  a,  Eq,  (77 )  reduces  to 


a 


(77) 


a  <  <  1  (78) 


Equation  (77)  is  evaluated  and  plotted  in  Fig.  5o 

It  is  evident  from  an  examination  of  Fig,  5  and  Eq,  (78)  that  the 
infinite  clipper  with  high-amplitude  dead-band  is  an  extremely  poor 
processor  to  use  if  the  variance  of  input  noise  is  large  relative  to  th9 
level  at  which  the  dead  zone  appears.  For  all  values  of  a  the  device  is 
inferior  to  a  simple  infinite  clipper. 

An  expression  for  the  signal-to-noise  performance  index  may  be  written 
following  the  calculation  of  the  required  averaged  functions  found  in 


Eq.  (56), 


rl  e 
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(79) 


Also 


■Willie  ■ 2  >Wil§)  f(!l  «S)  ■ 2 
J0 


|H2n(0)  f(0)  -  H2n(a)  f (a) 
n  a  0^1^2^o«, 

In  Eq,  (80)  f(a)  is  the  Gaussian  density  function  of  the  normalized 


(80) 


variable  a.  Since 


Eq.  (80)  becomes 


H2n+1(f)  g(e)  *>  2  Ho)  -  M«)  -  2  B2n(a)  (82) 

The  functions  in  Eq.  (82)  are  tabulated  in  Reference  8.  The  performance 
index  R2  in  Eq0  (83)  is  obtained  by  using  Eq.  (56).  Note  that  higher  order 
terms  have  been  added  to  Eq.  (83). 


The  series  in  Eq.  (83)  does  not  converge  well  for  values  of  a  <  1. 

For  a  =  0,  the  series  is  definitely  divergent  so  that  R2(0)  =  0.  Results 
from  Eq.  (83)  are  plotted  in  Fig.  5« 

From  Fig.  5j  it  can  be  seen  that  the  performance  of  the  infinite 
clipper  with  high-amplitude  dead-hand  approaches  that  of  the  simple  infinite 
clipper  for  large  a  =*  e^/a,  as  expected.  Below  a  »  3 j  however,  the 
performance  indices  decrease  rather  sharply,  which  indicates  that  the  device 
has  little  utility  as  a  processor. 
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Fig.  5  Performance  Indices  for  Infinite  Clipper 
with  High  Amplitude  Dead  Band 
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VII.  Average  Signal  Output--— Linear  Processors  and  Uon-Scyaare  Law  Detector 

In  preceding  sections  it  has  been  determined  that  the  use  of  linear 
processors  with  a  square-law  detector  maximizes  both  performance  indices. 

In  this  section  the  performance  of  a  non-square  law  detector  with  linear 
processors  will  be  examined.  It  would  be  desirable  to  obtain  expressions 
for  the  performance  indices  of  arrays  utilizing  nonlinear  processors  and  a 
non-square  law  detector,  but  the  mathematical  expressions  become  extremely 
cumbersome,  For  this  reason,  the  following  analysis  is  limited  as  indicated. 
Since  the  use  of  linear  processors  with  square-law  detectors  was  optimum, 
perhaps  the  use  of  linear  processors  with  non-square  law  detectors  is  not 
unduly  restrictive. 

With  reference  to  Fig,  1,  the  detector  function  y(x)  is  assumed  to  be 
an  even  function  of  x.  The  input  x  is  Gaussian  because  the  transducer 
outputs  are  all  Gaussian,  and  the  processors  are  linear.  For  the  "on 
target"  case,  the  signal  components  from  each  channel  are  completely 
correlated  and  the  noise  components  are  assumed  to  be  uncorrelated.  For 
the  "on  target"  case  x(t)  becomes 

x(t)  =>  Xg(t )  ■  M  s(t)  +  M^^n(t)  (8U) 

where  s(t)  and  n(t)  are  representative  signal  and  noise  time  functions  at 
any  transducer  output.  For  the  "off  target"  case,  the  signal  components 
are  uncorrelated,  and  so  are  the  noise  components.  We  have 

x(t )  -  x^t)  -  M^sft)  +  K^nCt)  (85) 
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The  "on  target"  and  "off  target"  array  outputs  are  defined  by  Eqs.  (86) 


and  (87  ). 


7q  ■  y(xQj  -  y(Ms  +  M^n) 

=  “  y(W2s  +  Jf^n) 


(86) 

(87) 


Both  y(xg)  and  y(x^)  can  be  expanded  in  a  Taylor  series  around  the  argument 
M^n  and  the  averages  taken  as  indicated  in  Eqs,  (86)  and  (87).  These 


operations  yield 


+  S£bL  +  ... 
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Only  even  order  terms  appear  in  Eqs,,  (88)  and  (8?)  since  the  average  of 
odd  order  derivatives  of  even  functions  yield  zero. 

The  normalized  signal  output  is 
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For  the  square  law  detector, 


A  M(S  +  N)  S  +  N 


Thus  the  signal  performance  index  is 


It  is  evident  from  an  examination  of  Eq.  (9U)  that  only  the  first  term  of 

g 

each  series  expression  is  not  negligible  for  small  values  of  For  small 

g 

values  of  t?  we  have 


x  =  M^n 
a2  =  MN 
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(95) 


The  expression  in  Eq,  (95)  does  not  yield  an  optimum  detector  function 
y(x)«  A  demonstration  of  the  validity  of  this  statement  is  found  in 
Section  IX,  in  which  is  evaluated  for  different  detector  functions „ 
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VIII,  Signal-to- Noise  Ratio — Linear  Processors  and  Mon-Square  Law  Detector 


For  the  purpose  of  mathematical  simplicity,  the  signal-to-noise 

s 

performance  index  is  determined  for  small  values  of 

Since  x(t)  is  Gaussian  Eq.  (26)  may  be  applied  to  yield  an 
expression  for  the  autocorrelation  function  of  the  detector  output: 


Ry(T)  -  y(xr 
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(96) 


The  first  term  of  Eq.  (96)  is  simply  the  square  of  the  average  array  output, 
and  all  other  terms  represent  the  autocorrelation  of  the  variational 

s  * 

component  of  y(t).  Thus,  for  low  values  of  Ry(T)  becomes 
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s' (T)  .  Y  — i—  (Iffl i2k  p21c(t) 

7  tl  »» 


00 


00 

(MN)2k  y^2kXM1/2n)  f  p2k(^)  dx 
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(98) 


For  the  square-law  detector  (y  =  x  ),  either  Eq.  (5U),  in  which 
K  “  1,  or  Eq.  (98)  gives  o 


aZA2  “  2  cc>p  M2  N2  /  pc(t)  dT 
0 


(99) 


The  definition  of  R^  in  Eq.  (5$ )  together  with  the  first  two  terms  in 
Eqs.  (88)  and  (89)  and  also  Eqs.  (35),  (36),  (57),  (98),  and  (99)  yield 
the  expression  for  R^. 
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An  alternate  form  of  Eq,  (100 )  is 
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For  the  function  y  •  x  it  is  observed  that 

y[2k](x)  e  q  for  k  >  1 


-1/2 
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(101) 


and  thus 

For  all  other  functions 


y^  k\x)  -  0  for  k  >  1 


y^2lt\x)  t 


0  for  some  k  and  x 


and  thus 


y^2^(x)  /  0  for  some  k 


2 

Since  all  terms  in  Eqs.  (100)  and  (101)  are  positive,  y(x)  -  Kx  therefore 
produces  a  maximum  value  for  R^,  and  the  square-law  detector  is  optimum  with 
respect  to  output  signal -to-noise  ratio. 

In  the  following  section,  R^  and  Rg  are  evaluated  for  specific  detector 
functions. 
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IX,  Evaluation  of  R^  and  R2  for  Spetjfio  Detector  Functions 

A  wide  class  of  detector  functions  is  described  by  the  general  power 
law  relation  found  in  Eq.  (102).  Note  that  y(x)  is  an  even  function  of  x. 


y(x)  =»  K|x| 
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(102) 


T2T  . 

Equation  (63)  is  helpful  in  evaluating  y  (x)  and  y(x): 

2 


00 


HxJ  -  2  Kx] 


x  _ ,  m-1 

m  i_  /  2?  dx  -■yT  K  am  2  2  /2L=-H  » 


vs 


(103) 


-X 


yt2J(x)  -  2  /  K  m(m  -  1)  xJ 


0 


m 


-2  1  J? 


vs 


dx 


■vr 


m-1 


"  "  rr,  9  ^  I E  ,**  — )  ! 


K  a1™’'"  m  2 


1  (10U) 


From  Eq.  (95)  we  obtain  R^(m): 
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(105) 


As  can  be  seen  from  Eqc  (105)  the  signal  performance  index  increases  linearly 
with  m. 

The  higher  order  average  derivatives  are  computed  below. 
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for  the  low-pass  spectrum,  and  is  less  than  £  for  other  spectra. 
The  signal-to-noise  performance  index  from  Eq.  (100)  becomes 


HgOn)  - 


»  -  2)2(m  -  i,)2(n  -  6)2  *  . 


-1/2 


(108) 


Equation  (108)  is  evaluated  and  plotted  as  a  function  of  m  in  Fig0  6. 

The  optimum  detector  function  from  the  standpoint  of  signal-to-noise  ratio 
is,  of  course,  the  square-law  device,  but  the  maximum  is  quite  broad,  as 
can  be  seen  in  Fig.  6. 

Another  detector  function  is  described  in  Eq,  (109)  and  shown  in 
Fig.  7. 
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— x1  <  X  < 

elsewhere 
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The  following  equations  yield  the  results  for  the  averaged  function  and 
the  averages  for  the  higher  order  terms  found  in  Eq.  (101). 
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where  p(-^is  defined  as  in  Eq,  (76),  Also, 
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Fig.  6  Signal-to-Noise  Index  Rg  for  Various  Power  Law  Detectors 

y  -  KlxT 
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Fig.  8  Signal  Performance  Index  vs,  {3 
for  Detector  with  Dead  Band 
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If  we  let  p  =  —  ,  the  performance  index  R^  from  Eq.  (95)  for  small  ^ 


becomes 
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i_L 


-\[2n 


P  for  p  <  <  1 
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(112) 


|  P2  for  p  >  >  1 


The  approximation  for  p  >  >  1  is  obtained  by  using  the  asymptotic  expansion 
for  1  -  P(p)  found  in  Reference  7 >  P»  VIII,  Equation  (112)  is  plotted  in 
Fig.  8. 

Using  the  result  in  Eq,  (ill)  with  the  general  result  for  R£  in 
Eq.  (101),  we  have 
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Equation  (113 )  has  been  evaluated  with  the  help  of  the  tables  in  Reference  8 
for  six  terms  of  the  series.  For  large  values  of  p  the  series  converges 
slowly.  The  results  are  shown  in  Fig.  9« 

Figures  8  and  9  show  that  the  performance  of  the  detector  with  dead¬ 
band  compares  favorably  with  that  for  the  square  law  device  in  the  region 

p 

where  l<p<2  or  1<  x  /MN  <  U. 
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X,  Conclusions 


The  investigation  in  Sections  IV,  V,  and  VI  has  shown  that  the  selection 
of  linear  processors  from  a  general  field  of  odd-funci/Lon  processors  is 
optimum  from  the  standpoint  of  both  performance  indices.  It  is  assumed  that 
a  square  law  detector  is  used  in  conjunction  with  the  odd-function  processors. 
The  evaluation  of  the  performance  indices  in  Section  VI  for  specific 
processors  shows,  however,  that  the  array  performance  is  not  greatly 
degraded  for  certain  processor  functions  which  deviate  significantly  from 
the  linear  function*  For  the  odd  power  law  device, 

R^  >  *6  and  R£  >  .75  for  0  <  m  <  3 

However,  performance  was  degraded  greatly  with  the  infinite  clipper  with 
dead  band  when  the  rms  value  of  the  signal  plus  noise  was  equal  to  or 
greater  than  the  dead  band  level  of  the  processor*  The  examples  indicate 
that  significant  performance  degradation  occurs  only  when  the  processor  has 
dead  space  or  zero  output  for  a  significant  range  of  the  input  variable. 

Thus  the  use  of  the  infinite  clipper  (m  »  0)  as  an  array  processor 
does  not  cause  significant  degradation  in  performance,  particularly  in  view 
of  the  other  advantages  treated  in  part  in  Reference  3« 

The  investigation  in  Sections  VTI,  VIH,  and  IX,  which  treats  an  array 
with  linear  processors  and  a  general  even-function  detector,  has  shown  that 
the  square-law  detector  is  optimum  from  the  standpoint  of  output  signal-to- 
noise  ratio.  However,  other  detectors,  particularly  power  law  devices  with 
exponents  greater  than  two,  yield  greater  signal  detection  indices  than  the 
square  law  device.  This  is  a  distinct  advantage  since  suppression  of  the 
"off -target  plateau"  in  the  typical  directivity  pattern  is  a  desirable  result. 


A  4- 40 


As  in  the  investigation  with  the  processors,  the  optimality  of  the 
square  law  device  with  respect  to  the  signal-to-r.oise  index  is  not 
critical.  For  example 


^2  -  »7J>  for  0  <  m  <  6 

The  results  for  the  dead-hand  detector  are  interesting  in  that 
performance  comparable  to  or  somewhat  better  than  that  for  the  square  law 
detector  could  be  obtained  if  the  dead-band  level  of  the  detector  is 
adjusted  properly.  For  example 

=  2,5  and  R2  =  .93  if  P  -  2 

Further  work  might  include  investigation  of  arrays  with  multipliers 
as  detectors  since  suppression  of  the  "off-target  plateau"  in  the  directivity 
pattern  is  certainly  desirable  for  small-signal  detection.  The  inclusion 
of  non-zero  correlation  between  noise  inputs  for  different  channels  seems 
necessary,  since  noise  cross-correlation  degrades  performance  in  this  respect. 


A/f-41 


Appendix  A 


The  terms  summed  over  m  in  Eq.  (22),  will  be  considered  a  remainder  term. 


n 


«*)  ■  E 


m=l 


(n-m )  [m-l] 


2  2 
a  w 


(n-m)  [m-l] 


2  2 
o  w 


w 


g  (0)  e  dw  -  jw  g’  (0)  e 


dw 


(A-l) 


For  contour  C,  w  =  A  +  jv  ,  and  for  contour  D,  w  =  -  A  +  jv  where  -00  <  v  <  00 . 
Rewriting  Eq.  (A-l)  we  have 
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Figure  A-l  shows  a  closed  rectangular  contour  defined  by  A  and  R. 

The  integrals  in  Eq,  (A-2 )  are  taken  over  the  long  sides  of  the  rectangle . 


v 


Fig.  A-l  Closed  Rectangular  Contour 

Since  the  integral  over  the  entire  closed  contour  is  zero  by  Cauchy's  integral 
theorem,  Eq,  (A-2 )  may  be  rewritten  in  terms  of  the  integrals  over  the  short 
sides  of  the  closed  rectangular  contour. 
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r?  M 

F(A)  -  )  g  (0)  lim 

R-^oo 

m=l 


r- >  [m~l] 

)  g  (0)  lim 
R^oo 

m=l 


A+jR 


(n-m) 


w 


2  2 
a  w 

2  i  /  2' 

e  dw  +  /  w  e  dw 


-A-jR  2  2 

'  /  \  a  w 

(n-m) 


-A+jR 


A-jR 


f  (n-m) 

(n-m )' 

d'd  -I 

a  R 

“  9 

<kjR) 

-  (-jR) 

>  e  -  (2A ) 

as  long  as  A  is  finite, 

since 

-,n 


-R 


lim  R  e  =0 
R-»a> 


(A-3) 
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Appendix  B 


It  is  given  that 


g^nl(e )  s  |  g^nJ(e)  f  (e )  de 


—00 


where 


f(e)  =  — — —  e 


je 
2  a 


(B-l ) 


(B— 2 ) 


Integrating  Eq.  (B-l)  once  by  parts,  we  have 

,  00  00 


[n-1] 


g[nJ(e)  =  g  "(e)  f(e) 


[n“l]  _T  -1 1 

g  (e)  (e  )  -de 


—00  „00 


00 


[n-1] 


=  0  -  /  g  (e )  fl  (e )  de 


HI/ 


_0O 


(B-3) 


In  the  same  fashion,  the  right  hand  side  of  Eq,  (B-3 )  may  be  integrated  by 
parts  another  (n  -  l)  times  with  the  result  found  in  Eq,  (B-U). 


00 


g[rJ(e)  =  (~l)n  /  g(e )  fLnJ(e)  de 


[n], 


(B-U) 


-OO 


[n-m]  [m-l] 

The  product  terms  of  the  form  g  (e )  f  (e )  always  disappear  because  of 

Cm-lj 

the  exponential  nature  of  f  (e). 


By  definition: 


dx 


n 


.1) 

2  \  -  (-l)n  Hn(x)  e  2 


(B-5) 


where  H  (x)  is  the  Hermite  polynomial  of  order  n.  Applying  the  definition 


n 


in  Eq,  (B-5 )  to  the  probability  density  function  in  Eq<>  (B-2 ),  we  have* 
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ftnJ(e)  -  (-l)n  o"n  Hn(|)  f(e)  (B-6) 

Combining  Eqs„  (B-6)  and  (B-U)  we  have 

00 

g^(e)  “  j  g(e)  o"n  Kn(|)  f  (e  )  de 
*00 

gCn](e)  =  o"n  Hn(|)  g(e)  (B-7) 


Given  below  is  a  table  of  the  coefficients  of  the  Harmite  polynomials 
up  to  order  10.  Important  recursion  relations  are  also  given. 
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3 
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0 
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0 

— 
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0 

l 
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0 

-l 
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■1 

0 
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0 

3 

i 

0 
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0 

3 
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l 

0 

-10 

0 

15 

0 

5 

.  - 

I 
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H 

0 

-15 

0 

U5 

0 

-15 

6 

■ 

1 

0 

-21 

0 

105 

0 

-105 

0 

7 

- 

1 

■ 

0 

-28 

0 

210 

0 

-U20 

0 

10$ 

8 

-36 

0 

378 

0 

-1260 

0 

9A5 

0 

9 

1 

0 

-U5 

0 

630 

0 

-3150 

0 

U725 

0 

-9U5 

10 

Table  B-l  Coefficients  for  Hermite  Polynomials 
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dH  (x) 

H, ..  (x)  =  x  H  (x)  -  — - 


VI' 


n 


dx 


x  Hn(x)  -  n  Hn_1 (x ) 


From  Table  B-l,  for  example, 


H?(x)  -  x7  -  21  x5  +  105  x3  -  105  x 


(B-8) 

(B-9) 
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Summary 


The  detectability  of  signals  consisting  of  Gaussian  noise  modulated  by 
a  periodic  function  of  time  is  considered  in  this  report.  It  is  assumed 
that  the  modulation  has  a  precisely  known  form.  Detection  is  considered 
in  detail  for  four  different  examples.  These  are: 

1.  The  spectrum  of  both  the  background  noise  and  the  unmodulated 
signal  is  white  (uniform )j  the  receiving  array  consists  of  a 
single  hydrophone. 

2.  Same  as  (l)  except  that  the  receiving  array  consists  of  K 
hydrophones , 

3.  The  spectrum  of  the  background  noise  is  white,  that  of  the 
unmodulated  signal  drops  off  with  the  second  power  of  frequency 
above  a  certain  point.  The  receiving  array  consists  of  K 
hydrophones  . 

ii.  The  spectrum  of  the  background  noise  drops  off  with  the  second 

power  of  frequency  above  a  certain  point,  the  unmodulated  signal 

spectrum  is  white,  the  receiver  consists  of  a  single  hydrophone. 

In  all  of  these  cases  the  improvement  in  detectability  is  very  modest.  In 

cases  1,  2,  and  U  the  ratio  of  detectability  index  of  modulated  to 

1  2 

unmodulated  signal  is  1  +  ^  b  ,  where  b  is  the  modulation  index  which  is 
less  than  unity.  The  results  of  case  3  are  even  less  favorable  <,  For 
cases  1  and  2  the  structure  of  the  optimum  detector  has  been  obtained; 
it  turns  out  to  have  essentially  the  form  that  is  used  in  practice. 
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I .  Introduction 


The  signals  that  have  teen  considered  in  most  of  the  previous  studies 
of  underwater  detection  were  assumed  to  be  samples  of  stationary  gaussian 
noise. However,  in  practice  signals  having  a  certain  amount  of 
deterministic  structure  are  occasionally  observed.  Utilization  of  this 
deterministic  structure  by  the  detection  system  should,  in  general,  permit 
greater  detectability.  In  the  present  report  the  type  of  deterministic 
structure  that  is  considered  is  a  low-frequency  amplitude  modulation  of 
the  signal  power.  Such  a  modulation  can  be  observed  quite  clearly  in  the 
signal  produced  by  a  ship  with  a  rotating  propeller.  The  general 
appearance  of  such  a  signal  is  shown  in  Fig,  1,  The  signal  appears  to  be 


a  random  noise  whose  power  varies  at  the  propeller  frequency. 

We  assume  in  this  analysis  that  the  signal  has  a  gaussian  probability 
distribution  with  zero  mean,  and  that  the  variance  is  a  sinusoidal  function 
of  time  whose  amplitude,  frequency  and  phase  are  known  precisely.  In  practice 
the  modulation  is  probably  not  exactly  sinusoidal,  and  even  if  it  were, 
its  amplitude,  frequency  and  phase  would  be  unknown  a  priori  and  would  have 
to  be  estimated  from  the  signal.  Thus  the  situation  that- is  analyzed  is 
considerably  more  favorable  than  the  true  situation.  The  analytical  results 
can,  however,  be  considered  as  providing  an  upper  limit  on  the  attainable 
detectability  -  an  upper  limit  which  cannot  actually  be  reached  in  practice. 
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II.  Detection  of  Gaussian  Signals  in  Gaussian  Noise 

The  problem  of  detecting  gaussian  signals  in  gaussian  noise  has  been 
studied  extensively  in  the  liter  ature.'>,^>  The  most  complete  results  are 
probably  those  of  Middleton.^*  Middleton  takes  the  point  of  view  that  the 
received  signal  is  given  by 

x(t)  =  y(t )  +  n(t )  (1) 


where  y(t)  is  the  signal  at  the  "transmitter”  and  n(t)  is  the  noise. 
Although  the  received  signal  is  usually  observed  continuously  it  is 
convenient  to  assume  that  the  signal  is  sampled.  In  this  case  Eq,  (l) 
becomes 


x  =>  y  +  n 


(2) 


where 


x(t^) 

y(t^) 

n(t±) 

X  * 

x(t2). 

• 

>  I  = 

0 

,  n  - 

• 

• 

• 

0 

e 

_y(tn)_ 

0 

• 

Ln(tnl 

The  samples  do  not  have  to  be  taken  at  successive  time  intervals;  in 
fact  the  idea  of  signal  samples  extends  directly  to  the  simultaneous  observation 
of  several  different  signals,,  Thus,  if  signals  received  by  an  array  of 
K  hydrophones  are  to  be  considered,  the  signal  and  noise  vectors  of  Eq.  (2) 
would  be  of  order  n  “  mxK  ,  where  m  is  the  number  of  time  samples  taken  on 
each  of  the  K  hydrophones. 

In  general  the  samples  are  not  independent.  Thus,  if  £  is  assumed  to 
be  a  gaussian  variable,  then  the  n-fold  probability  density  function  of  jr  is 
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(3) 


p(y) 


(V2n‘f|P|1//2 


f  2  £ 


where  P  is  the  covariance  matrix  of  the  transmitted  signal  samples;  that  is, 
the  i j  element  of  P  is  given  by 


'  y(V  y(V 

|  P  |  is  the  determinant  of  P  and  y  is  the  transpose  of  y.  Similarly  the 
n-fold  probability  density  of  the  noise  is  given  by 


q(n) 


(V^)nlel1/2 


1  *  -1 

-  vf  n  Q  n 


(h) 


Hence  if  signal  and  noise  are  independent  gaussian  processes  the  probability 
density  function  of  the  received  signal  x,  given  that  there  is  a  transmitted 


signal,  is 


f(x/y) 


V 2n  |det(P  +  Q) 


P72 


-  \  x'  (p+Q)_12 

e 


<5) 


If  there  is  no  transmitted  signal,  the  probability  density  of  the  received 
signal  ia  simply  that  of  the  noise,  i„e,, 


f(x/0) 


_ 1 

*\Zf2n  (det  Q)1^2 


-  |  x  'q_-1x 
e 


(6) 


Hence  the  likelihood  ratio  is 


*x)  - 

~  f(x/0) 


det  Q 

det  (P  +  Q) 


(7) 
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In  the  threshold  detection  case  the  signal-to-noise  ratio  is  generally 
very  low,  so  that  the  elements  of  the  P  matrix  are  all  very  much  smaller 
than  the  corresponding  elements  of  the  Q  matrix.  In  this  case  Eq,  (7 )  can 
be  simplified  as  follows: 


The  coefficient 


det  Q 

1/2 

a 

l 

1 

“ 

>  1 

det  (P  +  Q) 

det 

[i  *  £  s'r 

1/2 


It  can  be  shown  (Middleton^*  Chapter  17  )  that  for  any  arbitrary  matrix  A 
and  arbitrary  scalar  variable  "X 


(8) 


dct 


I  +  ifA  J  -  ®xp  £  (-l)"*1  ^  te  An 


(9) 


m=l 


This  expansion  is  valid  so  long  as  for  the  largest  eigenvalue  of  k  the 
inequality  X^(,  <1  holds,  (The  "tr"  in  Eqe  (9)  means  "trace,"  the  sum 
of  the  diagonal  elements  of  the  matrix, ) 

Applying  Eq.  (9)  to  Eq,  (8)  we  set  &=  1  and  A  =  P  Q-^.  If  the 
signal-to-noise  ratio  is  small  all  eigenvalues  of  P  Q*"'1'  will  be  less  than 
unity,  so  that  Eq,  (9)  is  valid.  Hence  the  coefficient  becomes: 


det  Q 


det  (P  +  Q) 


1/2 


exp/ 


OO 


-i>>> 


m+1 


tr  (P  Q-1)" 


m 


m^l 


2 

-  |  tr  P  Q"1  +  ^  tr(P  Q"1)  -  ... 


(10) 
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The  exponent  in  Eq.  (7)  may  be  expanded  as  follows: 


-  |  x'  (P  +  Q )_1-Q_1  x  -  -  |  xV1  (I  +  P  Q"1)  -  I  x 

-  -  \  xV1  I  -  P  Q"1  +  (P  Q"1)2 . I  x 

=  I  Q-1*  -  I  X 'q_1p  Q_1P  Q_1X  +  ...  (11) 

For  small  signal-to-noise  ratio  the  higher-order  terms  become  negligible 
relative  to  the  first  term.  However,  instead  of  discarding  all  the  high- 
order  terms  completely,  it  is  desirable  to  retain  the  second  term 
approximately  by  replacing  it  by  its  average  value,  where  the  average  is 
taken  for  signal  absent 0  Thus 


|  <i'S_1P  Q_1P  Q-1^  -  1 1  £  “ij  •  I  Y,  £  *13  *1) 


\  tr  (PQ"1) 


where  ^  is  the  ij  element  of  Q-XP  Q”XP  Q”“*  a 


Combining  Eqs.  (10 ),  (11 ),  and  (12)  with  Eq.  (7)  yields  then  the 


result: 


£(x)  -  exp|  i  x  Q”1?  Q-Xx  -  |  tr  P  QJ  -  ^  tr  (P  Q-1)2  +  ...  j 


In  the  standard  likelihood-ratio  test  the  likelihood  ratio  is  compared 
to  a  threshold  Kj  if  the  threshold  is  exceeded  it  is  decided  that  a  signal 
is  present.  Hence  the  conditional  false  alarm  probability  a  is  the 
probability  that  £(x)  >  K  if  in  fact  no  signal  i3  present*-  the  conditional 
miss  probability  p  is  the  probability  that  ^(x)  <  K  if  a  signal  is,  in 
fact,  present.  It  is  usually  more  convenient  to  consider  the  log  £(x);  thus: 


00 


log  K 


pQ(u)  du 


p1(u)  du 


(lU) 


log  K 


-00 


where  u  =  log  t(x),  pQ(u)  and  p^(u)  are  the  probability  densities  of  the 
log  of  the  likelihood  ratio  for  signal  absent  and  for  signal  present, 
respectively. 

For  small  signal-to-noise  ratio,  po(u)  and  p^(u)  corresponding  to  the 
likelihood  ratio  of  Eqe  (13 )  will  be  approximately  gaussian  as  a  result  of 
the  central  limit  theorem.  In  this  case  p^(u)  and  p^(u)  are  completely 
determined  by  the  mean  and  variance  of  log  £(x).  From  Eq.  (13 ) 


r 

(log  4(xj)  =  |  <x'q_1P  Q_1X£  -  \  tr  P  Q"1  -  |  tr  (P  Q"1) 


\  tr  Q-1P  Q-1Q  -  |  tr  P  f 1  -  tr  (P  Q_1) 


-  -  l  tr  (  P  Q"1) 


(15) 


'r 

j^log  4(x) 


/  -^.og 


K 


N 

jj  ((x ' Q_1P  Q~V  )^>  -  |  tr  P  Q"1  +  ^  tr  (P  Q"1)  ^  <x  V^P  Q_1x> 


\  tr  P  Q"1  +  \  tr  (P  Q_1 )  |  -  |  77  tr  (P  (f1 ) 


2l 


|  tr  (P  Q"1)  +  ^ 


-2 


tr  (PQ"1) 


-  ^  tr  (P  Q"1 )  tr  (P  Q“x )  + 

f 

l  tr  (P  g1) 


tr  (P  Q"1) 


|  tr  P  Q_1  +  -r  tr(P  Q“X) 


|  tr  (P  Q_1 ) 


(16) 
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Similarly 


log  Z (x ) 


|tr 


(P  Q"1-) 


(17) 

(18) 


Thus  P0(u)  and  p^(u)  are  gaussian  distributions  having  the  same  variance, 

and  their  means  are  symmetrically  located  from  the  origin  by  a  distance 
1  —1  ^ 

^  tr  (P  Q  )  .  Substitution  into  Eq,  (lU)  then  yields 


where 


1 

2 


1  (n4)2'  1 


log  K 


]/tr  (Pa_1f| 


(19) 


It  is  clear  by  inspection  of  this  equation  (for  instance  with  K  =  1)  that 
a  and  p  both  decrease  as  tr  (P  Q-^")  increases;  thus  tr  (P  Q**^")  can  be 
considered  the  figure  of  merit  of  the  detection  system« 

It  is  also  interesting  to  note  that  the  "output  signal -to-noise  ratio" 
may  be  defined  by  r  y  4  y  _2 


<Jlog  ^(xj)  -  <Jlog  Z  (x  )^ 


N 


log  Z(x) 


Q- 


log  Z(x) 


|  tr  (PQ)2  (20) 


N 


1  —1 

Thus  tr  (P  Q“  )  can  be  identified  with  a  performance  index  via  a 
suitable  definition  of  signal -to-noise  ratio.  This  index  is  computed  for 
a  number  of  specific  examples  in  the  following  section. 


Ill,  Detection  of  Amplitude -Modulated  Noise  Signals 

An  amplitude-modulated  noise  signal  of  the  type  pictured  in  Fig,  1 
can  be  written  in  the  form: 

y(t)  -  f(t)  z(t)  (21) 


where  f(t)  is  a  deterministic  modulating  function  and  z(t)  is  a  stationary 
random  noise  function.  Then  the  general  element  of  the  P  matrix  i3 

-  <y(V  y(V> 

-  f(t±)  f(tj)  <z(t±)  f(t±)  f(tj)  Rz(t±  -  tj) 


(22) 


where  (t^  -  tj)  is  the  autocorrelation  function  of  the  stationary 
noise  z(t). 

Case  1:  Single  Hydrophone,  Uniform  Spectra  for  Signal  and  Noise 
If  the  noise  is  white  and  stationary,  the  Q  matrix  is  a  diagonal 
matrix  and  can  be  written  in  the  form 

Q  -  N  I 


(23) 


where  I  is  the  unit  matrix. 

Also,  if  the  signal  spectrum  is  uniform,  and  if  one  considers  only 
one  hydrophone,  then 

v*±  -  v =  s  (2U) 


ij 


where  S  is  the  signal  power  in  the  unmodulated  signal  and  6^  is  the 
Kroneker  delta;  6^  »  1,  6^  =  0,  i  /  j.  Then  the  P  matrix  is  diagonal, 
and  its  general  element  is 

Pit  -  S  f2(tt) 


(25) 
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The  matrix 


(p  2"1 ) 


is  then  clearly  also  diagonal, 


and  its  general  element  is 


7 


(26) 


so  that 


tr  (P  Q"1 ) 


S2 

7 


n 

Ls 

i=l 


fk(\) 


(27) 


To  perform  the  summation  a  particular  modulating  function  must  be  assumed. 

In  order  to  provide  a  fair  comparison  between  detection  of  modulated  and 
unmodulated  signals,  the  modulating  function  should  be  such  that  the  average 
noise  power  is  the  same  with  modulation  as  without.  Hence  let 


f2(t) 


1  +  b  cos  cot 


(28) 


The  modulating  frequency  is  assumed  to  be  relatively  low  so  that 
co(t^+^  *■  t^)  <  <  1.  Equation  (27)  becomes 


tr 


(P  Q"1)  - 


2  n 
i=l 


(l  +  b  cos 


N 


— 2 

AN 


(l  +  b  cos  cot)  dt 


where  A  is  the  time  interval  between  samples j  i.e.,  A  »  t^+^  “ 


*4 


(29) 


if  the  observation  interval  T  ■  nA  lasts  for  a  large  number  of  cycles  of 
the  modulating  function,  so  that  cJT  >  >  1, 
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If  there  were  no  modulation,  b  =  0.  Thus  the  ratio  of  detectability 
index,  tr(P  Q“^)  ,  of  modulated  to  unmodulated  signal  is 

r  «  1  +  Tr  (30) 


Since  b  cannot  be  greater  than  one,  it  is  sden  that. the  improvement  in 
detectability  is  rather  minor,, 

This  result  would  appear  to  violate  intuitive  concepts, since  practical 
experience  with  signals  of  the  form  shown  in  Fig,  1  indicates  them  to  be 
much  easier  to  detect  in  a  noise  background  than  steady  noise  signals. 
However,  tnis  is  true  only  if  the  signal-to-noise  ratio  is  quite  large. 

For  the  threshold  situation  considered  here,  where  the  signal-to-noise 
ratio  is  much  less  than  unity,  the  result  given  in  Eq.  (30 )  will  also 
properly  describe  the  practical  detectability. 

The  Optimum  Detector 

The  optimum  detector  is  a  likelihood-ratio  detector ;  i.e.,  it  performs 

the  operations  described  by  Eq,  (13 )«  The  only  term  of  that  equation 

1  1  -1  -1 

containing  the  signal  is  the  first  term  ^  x  Q  P  Q  x  ;  the  ether  terms 
are  constants.  For  the  white-noise,  white-signal -spectrum,  single- 
hydrephone  case  considered  hero 


n 


i  xV‘1P  Q_1x  -  ~  ~~  ^  (1  +  b  ccs  cd^)  x(t±) 

11  . 


,2 


ijl 
T 

1  (l  +  b  cos  at )  x"  (t )  dt 


2N  A 


(31) 


Thus  the  optimum  detector  squares  the  signal,  multiplies  it  by  the  modulating 
function  and  than  integrates  the  result.  This  is  shown  diagrammatically 
in  Fig.  2, 
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Fig,  2  Optimum  Detector 

The  part  of  the  detector  shown  inside  the  dashed  box  can  be 
implemented  by  a  linear  filter  matched  to  the  modulation  envelope.  Thus 
the  optimum  filter  has  essentially  the  form  that  one  would  use  intuitively 
on  a  signal  of  the  form  shown  in  Fig,  1, 

Case  2;  Uniform  Signal  and  Noise  Spectra,  Hydrophone  Array 
Here  we  assume  that  there  are  K  hydrophones  in  a  three-dimensional 

7 

isotropic  noise  background.  For  such  a  background  Faran  and  Hills  have 

shown  that  the  cross-correlation  function  for  noise  picked  up  at  elements 

separated  a  distance  d  ■  ex  ,  where  c  is  the  velocity  of  sound,  is 

s 

VT 

R(t,t  )  a  -±_  I  rV)  dT*  (32) 

8  2t 

s  J 

-T  +T 
S 

where  R*(t)  is  the  autocorrelation  function  of  the  noise  at  any  one 
hydrophone , 
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For  white  noise 


R  (t) —  6(t) 
2W 


(33) 


where  N  is  the  noise  power  and  W  the  band  width. 
Hence 


T  +T 

■s 


r(t^8)  - 


N 


UVfr, 


.  1  .  I 

6(t  )  dT 


-T  +T 
3 


N 


h&. 


N 


8Wt 


for 


for 


-T  <  T  <  T 
8  S 


T*T  or  To  -  T 
S  8 


otherwise 


If  white  noise  of  infinite  bandwidth  but  finite  power  is  assumed,  W->-oo 


(3U) 


so  that  R(t,t  )  =  0  for  all  t  /  0.  However 
S  8 


R(0,0)  =  N.  (35) 

Thus  Q  is  a  diagonal  matrix  of  order  m  »  K. 

The  elements  of  the  P  matrix  are  computed  under  the  assumption  that 
the  array  is  steered  "on  target,"  since  this  gives  maximum  detectability. 

Under  these  conditions  the  general  element  is 

2 

P^.  =  S  f  (t^)  for  ji  -  j  j  -  m,2m,3m...Km 

=■  0  otherwise  (36) 

where  m  is  the  number  of  observations  and  K  the  number  of  hydrophone  elements. 

_!  2 

Hence  the  diagonal  elements  of  the  matrix  (P  Q  )  are 
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for  1  <  i  <  m 


(P  a"1/] 

li 


S2  k 

K  %  f  4(t,  ) 

N 


K  %  fU(t.  ) 
N2  1“m 

K  $  fU( w 


for  itH-1  <  i  <  2m 


for  2m*  1  <  i  <  2m 


Therefore 


S2  k 

K~f 4(t.  ,v  )  for  (K-l }m+l  <  i  5  Km 
jj*.  l-lft-ljm 


m 


t.rCPQ”1)  aK2%£fU(ti) 


(37) 


n  r , 

i^i 


If  we  again  let  f  (t^)  =  1  +  b  cos  »t^  and  replace  th9  summation  by  an 
integration,  we  obtain 

T 

2  2  f 

tr  (P  Q”1 )  =  K2  ~  (1  +  b  cos  dt )  dt 

n  I 


2  T  S‘ 

K  1  7 


b2 1  2  S2  /  n  .  b2 

1  +  —  =  K  m  ~~x  1  +  — 

\  2 '  r  \  2i 


Thus  the  ratio  of  detectability  index,  tr  (P  Q"^)  ,  of  modulated  to 
unmodulated  signal  is  as  before 

■U2  ' 


(38) 


(39) 


Also  it  is  clear  that  the  optimum  detector  would  consist  of  a  filter  of  the 
type  shown  in  Fig.  2  attached  to  each  hydrophone,  with  the  outputs  of  all 
the  detectors  summed  to  give  the  total  output. 
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Case  3:  Harkov-Type  Signal  Spectrum  in  White  Noise 

In  this  case  the  power  spectrum  of  the  unmodulated  signal  is 


n  /  \  2aS 

"  ~2 - 5 

a  +  co 


•where  S  is  the  signal  power.  The  corresponding  autocorrelation  function  is 


R  (t)  =>  S  e 


For  the  moment  consider  only  a  single  hydrophone.  Then,  with  modulation  on 
the  signal  the  general  element  of  the  P  matrix  is  given  by: 


°  S  f(ti)  f(t^)  P 
-a(ti+1-t.)  -aA 

where  p  =  e  =  e  and  where  A  is  the  time  interval  between 


samples  as  before.  Also 


Q  =»  N  I 


as  before.  Thus  the  matrix  P  Q~  is  equal  to 


pf2fl 

p%fi 


pV2 


2 


PV3 

pf2f3 


p3fifu 

p2f2fU 


n-1.  „  1 
*  ’  •  p  flfn 


„n-l„  „ 
P  f  i, 
r  n  1 


•  •  f* 


where  f^  is  shorthand  for  f (t^ ). 


A5-15 


Then 


2  2 

(P  Q'1)  -  % 

N 


n 


2  2(i"l )  X 


X 


X 


i-1 


X 


X 


r  2 


n 


5]fi2p2(2-i,^f12p2(i"2) 

i-3 


i-1 


X  f- 


X  X 


3 

) 

_ I 

i=l 


n 


£  r.V^-i)  t2pn  !-3) 


i-U 


where  only  the  diagonal  terms  have  been  worked  out  in  detail.  Thus  we 
obtain 


2  p  n  r  3 

tr  (P  Cf 1 )  -  2*  y 

-  IT  ^ 1 

J-1L  i-l 


n 


f2Vf2p2(3~).t/V  f  2p2<i-3) 


i-j+l 


S  *  f2(s)  dsl 


2  2 
N  A 


■■■a  "5  /  f2(s)  ds 
N  A 


3  T 

f2(t)  p2^3-t^dt  +  f f2(t)  p2^"S^dt 


s  T 

„  -2a(s-t)  f  9  -2a(t~s) 

f  (t)  e  dt  +/  r(t)  e  dt 


L°  3  (42 ) 

2 

We  again  assume  that  f  (t)  ■  1  +  b  cos  cot.  The  integrations  involve 
considerable  algebra,  but  are  otherwise  straightforward.  If  all  high- 
frequency  terms  are  ignored  -  i.e„,  if  we  assume  as  before  that  cJT  >  >  1 
the  result  is 


2  2  1  2  2 
tr  (P  2-1)  *  4f  1  ♦  -pLj 
N  A  a  \  co  +  ha 


(U3) 
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and  therefore,  the  detectability  index  ratio  of  modulated  to  unmodulated 
signal  becomes  in  this  case 


r 


2  2 

,  2aV 

1  +  -5 - 

co  +  Ua 


5 


(WO 


This  result  is  actually  less  favorable  than  the  one  for  a  white  signal 
spectrum,  which  it  approaches  for  a-^  or  for  co-^0  j  in  other  words, 
if  the  modulating  frequency  is  much  less  than  the  width  of  the  signal 
spectrum. 

This  result  extends  directly  to  the  multiple  array  case,  A3  was  shown 
under  Case  2,  the  Q  matrix  is  still  diagonal,  and  if  the  array  is  steered  on 
target  the  P  matrix  is  a  K  x  K  array  of  m  x  m  identical  submatrices  of  the 
type  considered  above  for  one  hydrophone.  Thus  we  obtain  directly  for  an 
array  of  K  elements  that 


tr  (P  Q”1)  *  K2  44- 

N  A  a 


1  + 


2a2b2  ' 


2.2 

00  +  ua 


(US) 


The  detectability  index  ratio  is  then  again  given  by  Eq,  (Wi), 

Case  lit  Uniform  Signal  Spectrum,  Markov-Type  Noise  Spectrum,  Single 
Hydrophone 

This  is  the  only  remaining  case  for  which  results  can  be  easily 
obtained.  The  P  matrix  is  diagonal  and  has  the  form 

1"  2 
'  f 

1 


n 


(U6) 
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The  Q  matrix  is 


2  3  U  n-1 

1  p  p  p-'  p  P 

_  2  3  n~2 

P  1  p  p  P  ...  P 

2  2  n-3 

P  P  1  P  P  ...  P 


n-ln-2 
P  P 


...  1 


where  p  =  e”  as  before# 

The  inversion  of  the  Q  matrix  is  fairly  easily  accomplished}  the 


result  is 


1  -P 


-P  1  +  P  -P 


-P  1  +  p  ~p 


N(1  -  P  ) 


■P  1 


Hence 


,  1'  -  2,  2  2 

fl  +  fl  f2  p 


[pQ-1] 


2  S2 


N2(1  -  p2) 


r  2,  2  2  ,  ,  lt„  ,  2  ,  2.  2  2 

"f2  fl  P  +  f2  ^  p  '  "  ‘ 2  f3  p 


,  2-  2  2  2,.  2  ,  ,  3-  2  2 

^  ^  ^  *  P  /  *“  ^  P 


2  2  2  2 

-f  f  ,V  +  f 

n  n-l  n 
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where  the  off-diagonal  terms  are  not  shown.  Finally,  we  obtain 


tr(P  Q"1)2-  ^ 


-N2(l  -  p2) 


flUtrnUtZ  p2)^  p2  Z  A'  fi+l! 

i=2 


n-1 


i-1 


n 


i-2 


This  result  may  be  simplified  by  using 


-  2  2  a  /.  2U  .  196  f  2  v  .2  . 

i+1  '  1  at  1  *  at2  1 


«  2  -  2  d  /«  2 \  1  3 _  2 \  .2 

D  A  +  *2  2  ^  ""  °  •  ® 

3t  dt 


(50) 


Then 


tr(P  Q-1)2  ^ 


N2(l  -  P2) 


n 


n 


i=l 


i=»l 


a  ♦  t*  *  Pv  £  r,2  ^  (f,2) 


(f^  +  fnu)p2(i  +  P2)  +  p2a 


(51) 


For  large  n,  and  as  A  becomes  very  small,  the  only  significant  term  in  this 
expression  is  the  first  summation:  thus 

2  „2 


tr  (P  Q"1 )  ^  -S-J l-t.P1^)  v, 


(52) 


N2(l  -  p2)  i=l 
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p 

For  f^  =  1  +  b  cos  this  result  is  identical  with  that  obtained 
for  case  1,  Eq.  (27),  except  for  the  coefficient  multiplying  the 
summation.  Thus  for  T  ■  nA  very  large  so  that  cdT  >  >  1  ,  we  find  here 


that 


tr  (P  Q'1) 


S2(l 


Ju 


N2(l  -  p2) 


n 


1  + 


bc 


(53) 


The  ratio  of  detectability  indices  of  modulated  to  unmodulated  signal  is, 
as  for  case  1, 


r 


(5U) 


IV.  General  Conclusions 

The  four  examples  of  the  effect  of  modulation  on  the  detectability  of 
a  noise-like  signal  represent  the  simplest  cases  for  which  results  could 
easily  be  obtained.  Other  examples,  such  as  a  Markov  signal  in  a  Markov 
noise,  etc,,  could  also  be  considered.  The  computational  complexity  begins, 
however,  to  grow  very  rapidly,  especially  when  hydrophone  arrays  are 
considered,  and  it  is  unlikely  that  the  results  would  turn  out  to  be 
qualitatively  much  different  from  those  that  have  been  obtained  in  the 
simpler  examples.  All  these  examples  show  that  the  improvement  in 
detectability  resulting  from  modulation  of  the  signal  is  quite  small,  with 
the  maximum  increase  of  detectability  index  being  no  more  than  $0  %,  This 
value  is  obtained  if  the  modulation  index  b  is  unity  and  if  the  modulating 
frequency  is  very  much  lower  than  the  bandwidth  of  the  signal  spectrum. 

Also  implied  in  all  the  computations  is  that  the  form  of  the  modulation 
envelope  be  known  precisely;  this  would  not  normally  be  true.  Thus  the 
modest  improvement  in  detectability  that  has  been  computed  is  actually  a  very 
optimistic  estimate,  and  the  improvement  that  could  be  obtained  in  practice 
would  be  considerably  smaller. 
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If  the  background  noise  is  white,  the  results  extend  directly  to 
multi-element  arrays.  The  case  for  non-white  background  noise  in  a 
multi-element  array  has  not  been  completely  analyzed,  but  it  is  doubtful 
that  it  would  turn  out  to  be  radically  different  from  the  cases  that  have 
been  considered  in  detail. 

Qualitatively,  the  reason  for  the  poor  results  that  are  obtainable 
is  that  in  spite  of  the  coherent  modulation  the  detection  is  basically 
incoherent.  This  means  that  the  detector  is  sensitive  only  to  the  power 
in  the  signal  and  cannot  take  into  account  detailed  signal  structure.  The 
effect  of  the  amplitude  modulation  is  simply  to  increase  the  signal-to- 
noise  ratio  above  the  average  at  certain  times  and  to  decrease  it  at  other 
times.  The  detector  can  take  advantage  of  this  by  increasing  its  gain  when 
the  signal -to-noise  ratio  is  high,  and  by  decreasing  it  when  it  is  low0 
However,  even  if  full  advantage  is  taken  of  the  power  fluctuation,  the 
detection  cannot  actually  be  as  good  as  that  which  would  be  obtainable  on 
an  unmodulated  signal  having  a  constant  signal -to-noise  ratio  equal  to  the 
maximum  occurring  during  the  modulation  cycle.  Thus,  since  the  modulation 
peak  cannot  exceed  the  average  value,  the  maximum  improvement  in 
detectability  cannot  exceed  a  factor  of  2;  for  sinusoidal  modulation  it 
turns  out  to  be  a  factor  of  l„5o 

In  view  of  the  foregoing  discussion  it  seems  clear  that  major 
improvements  in  signal  detectability  can  only  be  obtained  by  some  form 
of  coherent  detection.  This  would  require  that  the  detailed  structure 
of  the  signal  be  taken  into  account.  There  is  some  indication  that 
coherent  wave  patterns  resulting  from  machinery  noises,  banging  of 
snorkel  tubes,  etc.,  might  be  found  at  the  low  end  of  the  signal  spectrum, 
A  suggested  avenue  for  further  study  is,  therefore,  the  investigation  of 
these  lower-frequency  coherent  signals. 


A5-21 


List  of  Symbols 


A  ■  any  matrix 

aij  °  general  matrix  element 

a  *  break  frequency  of  first-order  Markov  spectrum 
b  -  modulation  index 
c  =»  velocity  of  sound  in  water 
f(t^)  =  f^  =  modulating  function 

f(^/y)  -  conditional  probability  density  of  x  given  £ 
Gz(o)  =  power  spectrum  of  unmodulated  signal 
K  =  number  of  hydrophones 
K  »  likelihood-ratio  threshold 
=  likelihood  ratio 
m  ■  number  of  time  samples 
n  =  total  number  of  samples 
n(t)  =  noise  signal 
n  =  noise  sample  vector 
N  ■  noise  power 
P  =  signal  covariance  matrix 
=  general  element  of  P 
p(y)  =  probability  density  of  y 
Q  =  noise  covariance  matrix 
q(u)  =  noise  probability  density 

Rz(t)  =  autocorrelation  function  of  unmodulated  signal 
r  =  ratio  of  detectability  indices 
S  =  signal  power 
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T  =*  length  of  observation  time 

W  =*  noise  bandwidth 

x(t )  ■=  received  signal 

y(t )  =  transmitted  signal 

z(t)  =  unmodulated  transmitted  signal 

x  =  received  signal  sample  vector 

£  =  transmitted  signal  sample  vector 

A  =■  time  interval  between  successive  samples 

Tg  =  time  difference  between  hydrophones 

co  =■  frequency 

det  A  =  determinant  of  A 

/  \  =  conditional  averaging  operation  -  noise  only  present 

N 

/  ^  =  conditional  averaging  operation  -  signal  and  noise  present 

'S+N 

tr  A  =  trace  of  A  =  sum  of  diagonal  elements  of  A 

i 

x 


transpose  of  x 
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I .  Introduction 


The  object  of  this  investigation  is  to  obtain  a  realistic  estimate  of 
the  cost  in  target  detectability  of  clipping  the  outputs  from  an  array  of 
hydrophones  prior  to  data  processing *  The  problem  is  approached  by  considering 
the  performance  of  likelihood  ratio  detectors  operating  on  clipped  and  on 
unclipped  data;,  Likelihood  ratio  detectors  are  chosen  because  they  are  known 
to  be  optimal  under  reasonably  general  conditions „  Data  processing  is 
presumably  accomplished  by  digital  means  so  that  the  hydrophone  outputs  must 
be  sampledo 

Results  have  been  obtained  only  under  the  following  assumptions: 

(1)  Signal  and  noise  are  both  stationary  Gaussian  processes, 

(2 )  Signal  and  noise  spectra  are  flat  over  a  frequency  band  0  <■  £  <  W 
cps  and  vanish  for  f  >  VJ. 

(3)  The  noise  disturbances  received  at  different  hydrophones  are 
uncorrelated  (hence,  because  of  (l),  statistically  independent). 

(U)  The  signal  to  noise  ratio  at  each  hydrophone  is  small, 

(5)  The  array  of  hydrophones  is  steered  on  target, 

(6)  The  sampling  period  is  ^  seconds  for  the  clipped  as  well  as  the 
unclipped  system. 

The  second  of  these  assumptions  is  probably  the  least  realistic.  It 
could,  in  principle,  be  relaxed,  but  only  at  the  expense  of  greatly  increased 
computational  complexity.  At  the  present  time  it  does  not  appear  that  the 
additional  insight  thus  gained  would  justify  the  effort. 

If  assumption  (2)  is  satisfied,  the  disturbance  (signal  plus  noise) 
received  by  each  hydrophone  is  completely  characterized  by  a  series  of 
samples  taken  at  intervals  of  ^  seconds,.  The  continuous  signal  can  be 
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reconstructed  without  error  from  the  samples  so  that  use  of  the  samples  in 
place  of  the  continuous  signal  involves  no  loss  of  information.  Furthermore, 
assumption  (l)  in  conjunction  with  (2)  insures  that  the  samples  are 
statistically  independent* 

If  the  samples  are  clipped  they  remain  statistically  independent,  for 
functions  of  independent  random  variables  are  independent  random  variables. 

It  is  therefore  only  necessary  to  consider  the  detection  problem  for  a  single 
sample  from  each  hydrophone.  Extension  to  the  2TW  samples  that  will  be  taken 
in  T  seconds  is  a  trivial  matter. 

In  evaluating  the  results  of  the  present  analysis  it  is  necessary  to 
keep  in  mind  that  the  optimal  processor  in  each  case  is  presented  with 
sampled  data.  In  the  unclipped  system  this  represents  no  loss  of  information 
for  the  reasons  indicated  above.  The  clipping  operation,  on  the  other  hand, 
changes  the  original  spectrum  defined  in  (2 )  to  one  that  does  not  vanish 
beyond  any  finite  frequency.  Sampling  at  intervals  of  7^  seconds  (or  at  any 
other  finite  rate)  therefore  may  involve  seme  loss  of  information.  Hence, 
if  one  does  not  regard  the  sampling  at  intervals  of  7^  seconds  as  a  fixed 
constraint,  the  performance  of  the  clipped  system  might  be  improved.  The 
price  for  such  improvement  is  basically  the  requirement  of  greater  bandwidth 
(or  sampling  frequency)  in  the  data  processing  system. 

II.  The  Likelihood  Ratio  Detector 

The  instantaneous  value  of  the  likelihood  ratio  is  the  quantity  A. 
defined  by  <P(v/S)> 

/V-  - 5-  (1) 

P(V/0) 
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P(V/S)  is  the  conditional.  probability  that  the  sample  set  received  by  the 

M  hydrophones  assumes  the  value  V  when  the  signal  component  is  known  to 

have  the  value  S,  V  and  S  are  M~dimensional  random  variables  (V^,Vg, o».V.v.) 

and  (S1,,S2S0..SM)  respectively.  The  notation  /  indicates  an  averaging 

S 

operation  with  respect  to  the  probability  distribution  of  the  signal. 

P(V/0)  is  the  conditional  probability  that  the  received  sample  set  assumes 
the  value  V  when  signal  is  known  to  be  absent 0 

The  likelihood  ratio  detector  operates  by  deciding  "signal  present" 
if  /\  exceeds  a  certain  threshold  value  and  K signal  not  present"  if  A 
falls  below  the  threshold.  A  measure  of  performance  can  therefore  be 
obtained  by  comparing  the  change  in  the  average  value  of  A  brought  about 
by  the  appearance  of  a  signal  with  the  typical  fluctuation  about  the  average 
value.  Under  assumption  (U)  the  fluctuation  is  almost  entirely  due  to  the 
noise  and  one  therefore  arrives  at  the  following  figure  of  merit, 

R  (2) 

D(A) 

A  E(A.)  is  the  change  in  the  expectation  (mean)  of  A_  du®  to  the  signal 
and  D(/\)  is  the  standard  deviation  of  A.  in  the  absence  of  signals, 

The  quantity  R  will  now  be  computed  for  both  the  clipped  and  the 
unclipped  instrumentation. 

III.  Clipped  Data 

If  the  hydrophone  outputs  are  clipped  and  there  is  no  signal,  all 
possible  combinations  of  (V,, have  the  same  probability.  Hence 

P(V/0)  -  2’M  (3) 
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Thus  the  denominator  of  A  is  a  constant  independent  of  V .  Now 


where  N  is  the  average  noise  power  at  each  hydrophone  and  use  has  been  made 
of  assumption  (5),  which  implies  that  =  Sg  “  o«»  aS^=s,  The  clipping 
levels  have  been  written  a3  -  1,  obviously  without  loss  of  generality* 

Using  the  Gaussian  property  of  the  signal,  the  numerator  of  A  can  be 
written  in  the  following  form: 

<P<T1  -  *  1*  V2  '  *  VM  *  1  - 

Z  (Xj+S  )2  (Xjj+S  )2 

~  2X  2N  SH 

©  ^  G  c  •  •  6 

(5) 


00 


00 


(2nN  f^^lznjs 


*1  —  ^ 


ds 


£  is  the  average  power  of  the  signal  whose  mean  value  is  taken  as  zero. 
The  integrations  in  the  various  are  all  identical  and  can  be  carried  out 
separately.  With  the  change  of  variable 


,  -  *2 
(x±  +  s) 

2N 


(6) 
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one  obtains 


9  (xL+s ) 

1  /  2W  .  1 

V  2nN  J  V  n 


e“y  dy 


Thus  Eq„  (5)  becomes 


,  ,  fV2N 

1  +  1  -y  , 

5  'vri  e  13 


1+1  _S_  _  1/  __3_  |3 

2  VT  VzN  ^V2N I 


(^1  -  i  l»  v2  -  i  1,-..  v„  -  1 1)} 


00  2 
r  s 


/  d3  e~  ^  +  — L  .  i/_sj  ...  ^  _  JL.  +  i/_2_ 

V  2n$  J  2  V2t?  31V2n1  2  V2N  3ll/2N 


K  is  the  number  of  hydrophones  having  positive  outputs*  Under  the  stated 
assumption  it  is  clearly  immaterial  which  of  the  hydrophones  have  positive 
outputs . 

Expanding  the  two  bracketed  expressions  in  the  integrand  of  Eqe  (8) 


according  to  the  Binomial  law  and  retaining  only  terms  up  to  the  power 

i  \y2lf 


one  obtains 


<r(*i  -  i  1,  v2  -  i  1,...  v„  -  i  1  )> 


00  2 

r  s 


wr  ♦ 


-j  M-2 

O 

IYn\ 

c 

s 

T 

2N 

A/2n5 


^Terms  in  higher  powers  of  -~r,  are  ignored  because  of  assumption  (U), 

V2N 
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f 

\ 

2-M< 

i.i 

n 

K(K-l)  *  (M-K) (M-K-l)  _  K(|(_K) 

2 

4  > 

1 

/ 

Thus  from  Eqs.  (3)  and  (9) 


(9) 


K(K-l)  .  (H-K)(M-K-1)  _  K(M) 
2 


1 

N 


(10) 


In  the  absence  of  signal  the  probability  that  K  hydrophones  have  positive 

-M  MI 

outputs  is  simply  2  '(h  k)i'"k'j  0  In  Presence  signal  the  probability 


is 


,-M  MI 


_  2 

(m~K)I  K] ^  1  +  n 

l 


K(K-l)  +  (M-K ) (M-K-l )  _  K(M_K) 
2 


N 


Hence  the  quantity  A  E(A)  introduced  in  Eq.  (2)  is  given  by 
r 

1 

N 


M 


*<A.)  -Y 


l.l 

n 


K=0 


.  K  (m-k  ) 

2 


>  x 


,~M  MI 


(M-K);  K I  n 


K(K-l)  +  (M-K ) (M-K-l )  _ 
2 


£ 

N 


M 


■Mz 


Mi 


(M-K) I  Ki 


<S 


..  k(M-K) 

2 


£ 

N 


+  4 

71 


[k(K-1)  +  (M-K)  (M-K-l)  „miMl 

2 

'll 

2 

v 

2 

l»j 

J 

(n) 


The  K  sum  can  be  evaluated  without  great  difficulty  if  one  recognizes  j^by 
inspection  of  the  binomial  expansion  for  (l  +  l)^"j  that 


M 

I 

KnO 


MI 


•»M 


(M-K)!  Ki 


(12) 
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e 

Consider  first  the  coefficient  of  in  Eq.  (11) 


The  three  sums  will  be  evaluated  separately.  In  sum  (l)  the  change 


of  variable  K-2  =  j  leads  to 


f-T  M-2 

®  =  Y  - - -  =  — -  byEq.  (12)  (lU) 

(M-2- j  )  j  j«  (M-2  ) j 


Again  using  Eq.  (12) 


_  Vn  ,  9 M-2 

©  •  S  -  “  — -  (15) 

(M-2-K  )j  K!  (M-2  ),? 


Finally  the  change  of  variable  K-l  =  j  and  use  of  Eq,  (12  )  converts 
sum  (3)  to  the  form 


2  M-2 
(M-2 ) J 


Thus  expression  (13)  vanishes. 


(16) 
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2 


Next  consider  the  coefficient  of 
M 

2-H£ 


% 

T 


Ka0 


K(K-l)  +  (M.-K)(M~K-l)  _  K(M_K) 
2 


M 


2”M^  ^  (M2  -  I4MK  +  liK2  -  M) 


K=0 


M 


in  Eq,  (11) 

2 

M J 

(M-K)J  K 1 

(M-K  )i  u 

m 

(M-K)i  KJ 


2-(Mf2)^  _  I4K (K-K  )J 

K=0 

M  M-l 

2“(M+2'MjJm2(M-1  )  y  - - - 8M(M~1)  y 

I  JK  (M-K )'  KJ  .H 


K=0 


M-l 

+  i6Z 

.K-l 


(M-K)J  KJ  ^  (M-K-l)J  (K-l)J 


© 

K(M-K) 


(M-K-l)J  (K-l)J 


(17) 


Sums  (§)  and  (5)  can  be  written  down  Immediately  from  inspection  of  Eqs0  (12) 


and  (16)  respectively: 


Mi 


© 


,M-2 


(M-2  )i 


Sum  (6)  requires  some  further  manipulation: 

M-l 

(6)  =,  V  KK-1)  +  ll  [(M-K-l )  +  ll 
^  (M-K-l  )J  (K-l)J 

M-2  M-l 

-  y  — - —  +  (M-i 

j^2  (M-K-2 ) J  (K-2 )J 
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>£ 

K=1 


(M-K-l )J  (K-l)J 


(18) 

(19) 


I 


1 


(20) 


M-U  M-l 

? - i -  ♦  (m-1)7 

£<>  fci 


(M-K-l)J  (K-l)J 


where  j  2  K-20 

Now  from  Eqs.  (12)  and  (16) 

©  .  *  mx 

(M-U)J  (M-2 ).' 


(21) 


Finally,  combining  Eqs.  (17),  (18),  (19)  and  (21) 


H  r 

•My  |K(K-l)  + 
K=0  L 


-.2 


iMiiMzl)  .  k(k-i) 


Mi 


(M-K)J  Ki 


^  H2  (M-l  )2  -  |  M2  (M-l  )2  +  ^  M  (M-l )  (M-2  )(M-3  )  +  M(M-1  f 


M(M-1 ) 
1 


-  i  M(M-l)  +  ^  (M-2 ) (M-3 )  +  (M-l) 


M(M-1 ) 


(22) 


Hence,  substituting  into  Eq0  (11 ): 


A  E(/Y)  -  \  M(M-l) 


5 

N  / 


(23) 


<r 


The  reader  will  note  that  Eq0  (23)  does  not  contain  terms  of  the  order  -jjj- 

The  question  therefore  arises  whether  an  error  was  not  introduced  when 

.  _  lU 

terms  of  the  order 


A/2N 


were  discarded  in  Eqc  (9),  for  terms  of  this  form 


would  yield  contributions  of  order 


5 


— *  after  integration  over  the  signal 


distribution.  It  is  a  straightforward  though  tedious  matter  to  demonstrate 

J12 


that  the  added  terms  of  order 


cancel  so  that  no  error  arises.  The 


details  of  this  demonstration  are  carried  out  in  the  Appendix, 
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The  next  step  is  the  computation  of  D(/\) 

d(A)  -  ]j e(A2)  *•  [e(A)T 


(2U) 


both  means  being  computed  in  the  absence  of  signals 

From  Eqe  (10)  and  the  vanishing  of  (13)  it  is  immediately  apparent  that 

E(A)  "  1  (25) 

Thus  only  E(/Y  )  remains  to  be  computed., 

M  (  1  /]2 

e(/y2)  -  2-m  f  ji  + 1  „  K0WC)  i  — 

"  L  2  J  w  J  (M-K)J  U 

(26) 

/  <  \°  C 

By  Eq.  (12)  the  term  in  |  is  unity0  The  term  in  -rr  is  identical  in 

'  /  <  '2 

form  with  Eq.  (13)  and  therefore  vanishes t.  The  term  in  1 


N 


is  identical 


with  the  corresponding  term  in  Eqa  (11)  so  that  Eq„  (26)  can  be  evaluated 
by  inspection  of  previous  results, 

S'2 


Hence 


e(A.2)  -i  +  4  m(m-d  (4 

n  ' 

D(A)  -  ~Vm(M~1) 


-S 

N 


(27) 

(28) 


S  W 

in  Eq.  (9).  This  demonstration  is  carried  out 

V  2N/ 


Again  it  is  necessary  to  show  that  no  error  is  introduced  by  omission  of 

_  |U 

the  terms  of  order 
in  the  Appendix0 

The  figure  of  merit  for  a  single  sample  from  each  of  the  M  hydrophones 
can  now  be  written  down  by  use  of  Eqs.  (2),  (23)  and  (28): 

“  ~  Vm(m-i  )1  £• 


R 


Cl  15 


(29) 


Pc^  stands  for  the  single  sample  figure  of  merit  in  the  clipped  case. 
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The  generalization  to  the  group  of  2TWK  samples  observed  in  time  T 

follows  immediately  because  of  the  independence  of  samples  resulting  from 

assumptions  (l)  and  (2).  If  the  figure  of  merit  obtainable  by  processing 

clipped  samples  for  T  seconds  is  designated  as  R  ,  one  obtains 

c 


■Vtw1  I  Vm(m-D1  4 


(30) 


IV  c  Unclipped  Data 


The  likelihood  ratio  detector  for  operation  on  unclipped  data  has 
been  analyzed  in  Report  No.  3.  As  pointed  out  previously,  no  information 
is  lost  through  sampling  under  the  assumptions  of  the  present  analysis. 
Hence  the  results  of  Report  No.  3  can  be  used  without  modification. 
Because  of  assumption  (3)  the  array  gain  is  simply  M,  Therefore  Eq.  (30) 
of  Report  Noc  3  reduces  to  the  simple  form 


A  (DC  output) 
D (output ) 


R 

u 


T  $ 

^  T 


_ (  O' 

=  Vtw  m  ~ 


(31) 


R^  is  the  figure  of  merit  for  the  likelihood  ratio  detector  designed  to 
operate  on  unclipped  data. 

By  comparing  Eqs.  (30 )  and  (31),  one  obtains  the  desired  value  for 
"decrease  in  detectability"  or  "cost  of  clipping" 


5=  .  2  (,; 

Ru  n  M 

2 

The  performance  ratio  given  by  Eq,  (32)  is  evidently  close  to  —  for  all 
but  very  small  values  of  M.  The  fact  that  R  /R  =0  for  M  =  1  is  due 
to  the  definition  of  "useful  output"  adopted  in  the  analysis.  The  quantity 


A6-11 


A  E  (/\_ )  is  the  change  in  the  average  value  o.f  /\  when  a  target  appears 
that  was  not  previously  present 9  It  is  not  the  difference  between  the 
average  values  of  A  on  and  off  target®  Thus  a  single  hydrophone  operating 
on  unclipped  data  is  capable  of  detecting  the  appearance  of  a  target  by 
noting  an  increase  in  received  power,,  A  single  hydrophone  operating  on 
clipped  data  obviously  has  no  such  capability, 

V®  Conclusion 

2 

The  figure  of  approxir, lately  —  for  the  "cost  of  clipping"  defined  by 

Eq.  (32)  is  not  unexpected.  It  was  shown  in  Report  No,  3  that  the 

likelihood  ratio  detector  designed  for  operation  on  unclipped  data  is 

identical  in  performance  with  a  standard  power  detector  if  signal  and 

noise  have  the  properties  assumed  in  the  present  analysis.  It  can  be 

shown  also  (see  e.g.  Report  No®  l;)  that  a  standard  power  detector  operating 

on  clipped  data  derived  from  an  input  whose  spectral  properties  satisfy 

assumption  (2)  on  page  1,  has  a  figure  of  merit  lower  than  that  of  a 

2 

detector  operating  on  unclipped  data  by  a  factor  of  approximately  — . 

What  has  been  demonstrated  here  is  that  this  degradation  must  be  regarded 
as  an  inherent  cost  of  sampling  and  clipping.  No  detector  operating  on 
clipped  samples  could  exhibit  less  degradation  in  performance®  On  the 
other  hand,  a  detector  operating  on  a  clipped  continuous  signal  (or  a 
signal  sampled  at  a  rate  above  2V7  samples  per  second)  might  perform  somewhat 
better  (see  comment  in  section  I).  Since  any  possible  improvement  is 
certainly  limited  to  a  factor  of  about  ^  and.  the  analysis  of  the  continuous 
case  is  not  simple,  no  attempt  has  been  made  to  date  to  carry  it  out. 
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Appendix 


In  Eq„  (10)  A  is  represented  by  the  first  two  terms  of  a  power  corVec 
in  .  Since  A  E(A)  j^Eqs,  (11),  (23 )j  has  a  lowest  order  term 


proportional  to 


S\2 


it  becomes  necessary  to  investigate  whether  the 
*\2 


vu  vi.  jnou.  v  u  T *  j.  v  vt/ooai.  jr  vw  vu&i  u  no 

/  ‘A2 

omission  of  terms  of  order  in  Eq,  (10)  causes  any  error  in  Eqc  (23) 

'I  ,  _  >ii 


Returning  to  Eq„  (8)  and  carrying  terms  up  to  order 
an  expression  of  the  form 


-M 

lV2N 


one  obtains 


V  2n^' 


ds  e 


s 


1  +  (a^s  +  a^s 


K  , 


1  -  (ojS  +  a 


M-K 


_00 


where  a^,  a^  are  numerical  constants. 

Expanding  the  bracketed  expressions  as  binomials: 

(p<Ti  ’ 1  l*  7a  '  1 1>-  vn  ■  *-  4  ‘ 


00  2 
s 


V  2.5 


ds  e 


25 


JOO 


1  +  K(a^s  +  a,s^ )  + 


k: 


(K-2 )2  22 


(a^s  +  a3sj) 


+  - — -  (a  s  +  as3)  +  - — -  (a^s  +  as3) 

(K— 3  )2  3.'  3  (K-U)i  U.'  3 


(A-l) 


1  -  (M-K)(a1s  +  a3s3) 


*  (a  a  *  a,*3/  -  (a  ,  *  a  a3)3 

(M-K-2 )2  22  3  (M-K-3 )J  32  1  3 


+  (M-K  )l - (  +  3} 

(M-K-U)2  i*J  3 


(A-2) 
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Only  terms  of  order  s4  in  the  integrand  are  of  interest  for  the  present 
purposes.  These  are 


-CO 


a  a  +  - QkVl —  gb  -  2K(M-K)  a.a. 

(M-K-2  )J  1  i  (M-K-U)l  Ui  1  J 


K(M-K)J  U  +  K l  (M-K)J  a  h  +  Kj  _ 


(M-K-3 )•  3J  1  (K-2)J  (M-K-2  )£  h  1  (K-2)J  1  3 


K|  (M-K)  U  KJ  a  U 

—  a_  +  a_ 

(K-3)J  3 i  (K-U)i  UJ 


2“M0  /S'2  KJ  (M-K)J  <! 


fr 


[(M-K-2)*  KJ  (K-l)J  (M-K-l)J  (M-K)J  (K-2)J 
1  .  1 


V3 


[  (M-K-U)l  KJ  UJ  (K-l)J  (M-K-3 )J  3J  (K-2  )J  (M-K-2)?  U 
1  .  1 


(M-K-l)J  (K-3  )J  3J  (K-U)J  (M-K)J  42 


ik» 


(A-3) 


It  is  apparent  from  Eq„  (11 )  that  the  contribution  of  the  above  term  to 
AE(/V)  is  simply 


M 


2-M-3  5 2  MJ  Y,  l 


K=0 

1 


[  (M-K-2  )J  KJ  (K-l)J  (M-K-l)J  (M-K)J  (K-2)JJ 

1  .  1 


ala3 


[  (M-K-U)J  KJ  4J  (K-l)J  (M-K-3 )J  3J  (K-2)J  (M-K-2)?  4 


(M-K-l)i  (K-3  )l  32  (K-U)J  (M-K )  J  hi 


(A-4) 


Where  factorials  of  negative  numbers  appear  in  this  expression,  the  term  in 
question  is  to  be  interpreted  as  zero* 


;/>-H 


It  is  immediately  apparent  from  Eqs.  (13)  to  (16)  that  the  coefficient 


of  a^a^  vanishes.  The 

coefficient 

u 

of  a^  can  be  rewritten  in  the  form 

M-U 

M-U 

M-U 

o  r _ i 

__2  V 

A  +  V  A 

'  j£0  (M-K-U)J  Ki 

Uj  L 

UJ  K=0 

(M-K-U)J  KJ  3J  (M-K-U)i  KJ  U 

i  -i  + 

12  3 


1 


l 


(M-U)J 


0 


(A-5) 


Thus  A  E(/\_)  is  unchanged  by  inclusion  of  the  higher  order  terms. 
Furthermore,  it  is  clear  from  Eq.  (26)  that  the  additional  term  of  order 

A  o 

in  E(/Y“)  is  proportional  to  Eq0  (A-U).  Hence  there  is  no  change 
in  E (A2). 
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I .  Introduction 


An  important  class  of  sonar  targets  has  the  following  two 
characteristics  s 

1)  The  signal  emitted  by  the  target  contains  periodic  components 

in  a  low  frequency  band  0  <  f  <  ,  The  frequency  of  these  periodic 

components  is  not  known  a  priori. 

2)  The  high  frequency  portion  of  the  signal  emitted  by  the 

target  (W^  <  f  <  )  has  essentially  the  properties  of  white 

Gaussian  noise  amplitude  modulated  by  the  same  low  frequency  periodic 
component. 

Since  the  ratio  of  total  signal  power  to  total  noise  power  at  the 
receiving  array  is  likely  to  be  low,  the  question  arises  whether  these 
special  properties  of  the  signal  can  be  used  to  improve  target 
detectability. 

Utilization  of  the  periodic  property  of  the  low  frequency  signal 
is  theoretically  attractive  but  encounters  the  practical  difficulty  that 
self -noise  at  the  receiving  array  also  has  periodic  components  in  the 
same  general  frequency  range.  In  the  absence  of  a  priori  knowledge 
concerning  the  frequency  of  the  periodic  target  signal  it  would  be 
necessary  to  monitor  continuously  the  periodic  components  of  self  noise. 

A  second  possible  procedure  deals  only  with  the  high  frequency  band 
W1  <  f  <  and  attempts  to  improve  detectability  by  using  the  fact 

that  the  signal  is  amplitude  modulated  by  the  same  low  frequency  periodic 
component.  This  possibility  was  investigated  in  Progress  Report  No.  5» 
which  showed  that  the  best  possible  detector  of  this  type  could  in 
general  provide  only  small  improvements  over  a  standard  power  detector. 
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The  present  report  investigates  a  scheme  that  attempts  to  use  the 
periodic  low  frequency  signal  component  without  requiring  knowledge  of  its 
frequency  or  monitoring  of  the  self  noise.  Thi3  is  accomplished  by 
demodulating  the  signal  in  the  high  frequency  range  and  employing  the 
detected  modulation  envelope  as  a  "local  oscillator"  which  can  then  be 
cross-correlated  with  the  low  frequency  component  of  the  signal.  The 
difficulty  with  this  procedure  is,  of  course,  that  the  "local  oscillator" 
signal  is  extremely  noisy, 

II.  Problem  Statement  and  Assumptions 

A  block  diagram  of  the  proposed  instrumentation  i3  shown  in  Fig.  1. 


block  D.C. 


Figure  1  Cross-Correlation  Detector 


The  output  from  the  transducer  array  is  split  into  two  parts.  The  low 
frequency  portion  (0  <  f  <  Wn  )  is  given  by 

W(t)  =*  n^(t)  +  A^cos  2nf^t  +  A2C0S  2nf2t 


n^(t)  has  the  properties  of  Gaussian  white  noise  with  average  power 
limited  to  the  band  0  <  f  <  A^cos  2nf^t  is  the  periodic  signal 
component.  Agcos  2nf2t  is  a  periodic  component  of  self  noise. 


(1) 
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The  high  frequency  portion  (W^  <  f  <  V,'2 )  of  the  array  output  is  given 
by 

x(t)  »  n2(t)  +  8(t)|l  +  a  cos  2nf.jtJ  (2) 

ng(t)  is  the  noise  component  and  s(t)Jl  +  a  cos  2nf^tj  the  signal  component. 
Both  ng(t)  and  s(t)  have  the  properties  of  Gaussian  white  noise  bandlimited 
to  <  f  <  Wg  ,  with  average  powers  N2  and  S  respectively,  a  is  the 
modulation  index  of  the  sinusoidal  amplitude  modulation.  In  most  cases  of 
practical  interest  a  <  <  1,  The  implicit  assumption  that  the  modulation 
envelope  and  the  low  frequency  sinusoid  are  in  phase  is  probably  not 
entirely  realistic.  However,  it  certainly  represents  the  mos*  favorable 
situation  and  will  therefore  serve  the  purposos  of  this  study. 

The  performance  of  the  cross-correlation  detector  of  Fig.  1  will  be 
evaluated  by  comparison  with  that  of  a  simple  power  detector  operating  on 
the  band  W.  <  f  A  block  diagram  of  the  latter  is  shown  in  Fig,  2, 


x(t)  -  n2(t)  +  s(t)[l  +  a  cos  2nf^t] 

11 

\j 

y(t )  =  x2(t) 

Av. 

filter 

squi 

arer 

1 

1 

(o,wo) 

Figure  2  Power  Detector 


Up  to  the  low  pass  filter  it  is  clearly  identical  with  the  high  frequency 
channel  of  the  cross-correlation  detector.  The  power  detector  uses  the 
change  in  average  output  due  to  the  appearance  of  a  signal  for  detection 
purposes.  The  performance  of  each  instrumentation  will  be  characterized  by 
a  noise  to  signal  ratio  defined  as  follows: 

bf  .  Mean  square  value  of  fluctuating  component  of  detector  output 

5  Square  of  average  signal  component  of  detector  output 

_  (3) 

^In  practice  would  be  a  relatively  low  frequency  compared  to  so 

that  most  of  the  total  signal  power  would  be  concentrated  in  the  range 
Wx  5  f  $  W2. 
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III.  Derivation 


The  mean  square  fluctuation  at  the  detector  output  is  calculated  most 
readily  by  integrating  the  output  power  spectrum.  Since  the  quantity  y(t) 
is  of  central  importance  in  each  instrumentation,  the  first  task  will  be 
the  computation  of  its  spectral  density  G  (f  ).^ 

y 


y(t)  «  x2(t)  “j^Ct)  +  s(t)Jl  +  a  cos  2nf^tj| 

-  n22(t)  +  2n2(t)  s(t)jl  +  a  cos  2nf^tj  +  s2(t) 


1  +  a  cos  2nf. 


I2 

i*]  <u) 


Let 


n2  (t)  =■  n2  (t)  -  N2 
s  (t)  =  sc(t)  -  S 


(5) 

(6) 


2  2 

where  N2  and  S  are  average  values  of  n0  (t)  and  s  (t)  respectively.  The 

!  t  2 

spectra  of  a,  (t)  and  of  s  (t)  can  then  be  written  down  immediately. 

QO 


G  i  (f)  =  2  /  G  (f ' )  Gn  (f  -  f')  df’ 
"2  )  n2  n2 


(7) 


-OO 


Ga*(f)  -  2  I  Ga(f')  Ga(f  -  f')  df  * 


-.0 O 


Since  n2(t)  and  s(t)  are  statistically  independent,  the  spectrum  of 

n2(t)s(t)  can  also  be  written  down  without  difficulty, 

00 

-JOO 


(8) 


(9) 


^The  notation  G  <f)  =  spectral  density  of  z(t)  will  be  used  throughout. 
Spectral  densities  zare  defined  in  such  a  manner  that 

f  G  (f)  df  =  e[z2]„ 

Jcc  " 

2W.  B.  Davenport  and  VI .  L,  Root,  Random  Signals  and  Noise,  section  12,2, 
McGraw-Hill  Book  Co.,  1958. 
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Before  proceeding  further  it  is  convenient  to  rearrange  Eq.  (U)  as 


follows: 


y(t)  =  Ng  +  S 


^  v. 


(t)  +  o' (t) 

M 

+ 

_ i 

L  J 

+  2a 


n„  (t )  s  (t )  +  s  (t ) 


a2  i 

cos  2rtf^t  +  -g-  s  (t )  cos  Unf^t 


yi. 


+  2a  S  cos  2nf^t  +  -jy-  S  cos  Urtf^t 


(10) 


y6 


Note  that  yn  is  the  only  time -dependent  term  in  Eqc,  (10)  that,  does  not 
vanish  as  a->0.  If  a  <  <  1  y0  will  therefore  be  the  dominant  component 
of  fluctuation  noise,.  Its  three  terms  are  uncorrelated  so  that  its  spectrum 
can  be  calculated  easily  from  Eqs.  (7),  (8)  and  (9),  The  result  is  shown 


in  Fig.  3c 


Figure  3  Power  Spectral  Density  for  ^2  lsee  Ecl*  0-0' 
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Results  for  the  Power  Pete c tor 

The  change  in  average  output  due  to  signal  for  the  power  detector  is 
5 imply  the  signal  component  of  y^.  Hence 


Under  the  very  reasonable  assumption  Wq  <  f^  the  terms  y^  and  y^ 

do  not  appear  in  the  power  detector  output  at  all*  Furthermore,  for 

a  <  <  1  terms  y^  and  y^  are  clearly  small  compared  to  yg*  In  practice 

the  final  averaging  filter  undoubtedly  has  a  bandwidth  VJq  satisfying 

WQ  <  <W  .  In  that  ease  the  spectrum  of  y2  may  be  regarded  as  essentially 

flat  over  the  band  G  <  f  <  W  and  one  obtains 

o 

N  2  +  2  N„S  *•'  S? 

]\[  =  2  G  (0)  W  ~  2  W  — - - -  (12) 

y2  0  °  V.f2  - 

Equations  (11 )  and  (12)  immediately  yield  the  desired  noise  to  signal 
ratio  of  the  power  detector 


Results  for  the  Cross-Correlation  Detector 

The  blocking  capacitor  in  the  high  frequency  channel  of  Fig*  1 
eliminates  the  term  y^  in  Eq.  (10).  Its  effect  on  the  other  components 
of  y(t)  will  be  regarded  as  negligible*  The  output  z(t)  of  the  blocking 
capacitor  is  therefore  given  by 

z(t)  =  y2(t)  +  y.Ct)  +  y^(t)  +  y^(t )  +  y6(t)  (lU) 
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z  (t )  is  one  input  into  the  multiplier*  The  second  input  W (t )  contains 
three  components.  Hence  the  multiplier  output  u(t)  is 


15 

u(t)  -  ^  Uj  (t ) 
i=l 

where 

u^(t)  =  ^(t  )y2(t) 

u^(t )  =  y^tjA^cos  Sn^t 

W 

O 

O 

CM 

+> 

CM 

1! 

u2(t)  -  TXj  (t  >y3  (t ) 

•Uy  (t )  =  y3(t)A1cos  2nf1t 

u12(t)  =  y3(t)A2cos 

u^(t)  =  ^(t  )y^(t) 

«8(t)  =  y^(t )A1cos  2nfxt 

^13  (t )  -  y^(t)A2cos 

u^(t)  -  n^tjy^t) 

u9(t)  -  y^(t  V^cos  2nf1t 

(t )  =  y^(t  )A2cos 

u^(t)  =  n^tJy^Ct) 

ttj_0(t)  ■  7^)^oa  2nf1t 

(t)  <=■  y6(t)A2cos 

(15) 


2nf2t 

2nf2t 

2nf2t 

2nf2t 

2nf2t 


(16) 


Inspection  of  Eqs.  (10)  and  (l6)  makes  it  clear  that  only  u^(t)  provides  an 
average  component  of  output c  Hence  immediately 


5  -  (a  \  S)2  (17) 

The  computation  of  is  simplified  by  observing  that  all  terms  of  Eq.  (l6) 
except  u^(t),  u^(t)  and  u^(t)  contain  factors  of  a  or  a  Therefore  to  a 
first  approximation  for  small  a  only  u^(t),  u^(t)  and  u^(t)  need  be 
considered.  As  in  the  case  of  the  power  detector,  is  assumed  to  be 
sufficiently  small  so  that  only  the  zero  frequency  value  of  G^(f )  is  of 
interest.  The  two  factors  of  each  of  the  significant  u.^(t)  are  statistically 
independent  so  that  the  corresponding  spectra  can  be  written  down  without 


difficulty. 


G  (0) 

V 


G  (f)  df 
y2 


(18) 
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or  using  0  (f )  as  given  by  Fig.  3 
y2 


(N2  +  3  f  N  /  W 

0,  (0)--2 - i 


w2  -  wi  M  w2  -  W1 


Similarly 


5U  (0)  '  —  K  <V  +  °7  *  —  (l - S>- 

U6  h  L  y2  1  y2  1  J  2  \  w2  "  Wi/  w2  "  W1 


0  (0) 


\  (n0  +  sr 


To  (f  )  ♦  0  (-f2)l  *-2-  1 - M  - - 

U  L  y2  ^  y2  *  J  2  (  W2  "  W!'  W2  “  \ 


u^(t),  u^(t)  and  u^(t)  are  uncorrelated  so  tnat  the  approximate  spectrum 
of  u(t)  is  simply  the  sum  of  the  three  spectral  components.  Recalling  the 
previous  assumption  that  f^,f 2  <  and  <  <  W2  one  can  neglect 


W2  -w1  w2  -  W;L  W2  -  wx 


compared  to  unity.  Then 


(N  +  S)2  /  A  2  A  2\ 

]\T^  2  W  — - -  N.  +  +  — 

x  °W-W  12  2 

2  1  \  ! 

Hence  the  figure  of  merit  for  the  cross-correlation  detector  is 


M  Ii  ,ia22 

%d.  a  w2-V  s'  V  2 -4. 


Conclusions 


The  ratio  of  the  two  figures  of  merit  j^Eqs.  (13)  and  (23  )~J  is 

w!c,. .  i  L  Ni  .*»*] 

Nfti  2a  A12/2  A^2 
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This  ratio  can  be  smaller  than  unity  only  if  a  >  0.707  and  then  only 
if  the  power  — — -  of  the  sinusoidal  signal  component  is  large  compared 

V 

to  N,  and  to  —  .  It  should  be  noted  that  the  condition  a  >  0.707 
1  2 

is  contrary  to  the  assumption  a  <  <  1  repeatedly  used  in  the  analysis. 
0n9  is  therefore  forced  to  the  conclusion  that  the  cross-correlation 
detector  is  decidedly  inferior  to  the  simple  power  detector  under  most 
operating  conditions  of  practical  interest. 

One  could,  of  course,  use  the  cross-correlation  detector  to 
supplement  rather  than  to  replace  the  power  detector.  However,  the 
results  of  the  present  analysis  indicate  that  the  gains  to  be  made  by 
such  a  procedure  would  be  slight. 
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